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XVI. Scattering of Electrons by Charged Dislocations in Semiconductors 
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ABSTRACT 


Experiments suggest that dislocations in germanium act as rows of 
closely spaced acceptor centres. In n-type material a dislocation accepts. 
electrons and becomes a line of negative charge surrounded by a eylin- 
drical region of fixed positive space charge. An array of parallel disloca- 
tions will have relatively little effect on electron mobility parallel to: 
the dislocations. The mobility normal to the dislocations will be reduced 
because: (1) The space charge cylinders deflect electrons and scatter 
the electron momentum normal to the dislocations; this reduces mean 
free time. Dislocation scattering is characterized by a mean free path 
that depends on direction. (2) Even when the dislocations are far 
enough apart that the mean free time is not affected, the macroscopic 
mobility may be reduced because the electrons cannot drift parallel 
to the applied field ; instead they follow curved paths winding between 
the space charge cylinders (which in typical cases are the order of a 
micron in diameter). 

This paper deals with both drift and Hall mobility when the current 
flows normal to the dislocations and a magnetic field acts (1) parallel to 
the dislocations and (2) normal to the dislocations. 


§ 1. INTRODUCTION 


THE experiments of Pearson, Read and Morin (1954) indicate that 
dislocations in germanium act like rows of closely spaced acceptor centres. 
In n-type material these acceptors accept electrons from the conduction 
band and the dislocation becomes a line of negative charge. Surrounding 
the dislocation is a cylinder of fixed positive space charge. The purpose 
of this paper is to study the effect of these charged dislocations on the 
motion of conduction electrons. . 

In a previous paper, the author (1954) discussed the geometry of dis- 
locations in the diamond cubic structure and described the acceptor type 
behaviour of dislocations by a simple model (due to Shockley 1953) 
aie ee ee ee 
* Communicated by the Author. 
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in which the dislocation acceptors are identified with the dangling bonds 


on the edge of the extra plane of a dislocation having an edge component 
(that is, a dislocation which is not parallel to its slip, or Burgers vector). 
The electrostatic energy in the space charge region was calculated as 
a function of the fraction f of the dislocation acceptors which accept 
electrons and acquire a charge g, where q is the electronic charge. 
Approximate formulas were derived giving f as a function of temperature 
for given values of the chemical donor and acceptor concentrations and 
the spacing c between the acceptor centres on the dislocation. Several 
experiments were suggested to test the theory. The simplest array of 
dislocations on which to make experiments is the array produced by 
plastic bending. If the bending is done at a high enough temperature 
that there are no residual macroscopic stresses, then the dislocations are 
all parallel (and parallel to the axis of bending) and the density of 
dislocations N per cm? is proportional to the curvature. Such known 
and simple arrays were introduced by Gallagher (1952) and Pearson, 
Read and Morin (1954) and observed by the etch pit technique of Vogel, 
Pfann, Corey and Thomas (1952). 

As discussed by the author (1954), the easiest way to test the predicted 
} vs. T relation is to measure the conductivity parallel to the dislocations, 
so that scattering by the dislocations is a minimum. In this paper, 
however, we are interested in scattering by the dislocations and we 
consider the case where the current is flowing normal to the dislocations. 
We also deal with the Hall effect when the magnetic field is (1) parallel to 
the dislocations and (2) normal to the dislocations. 

The dislocations affect conductivity and Hall angle because : 

(1) The average concentration of conduction electrons is reduced by 
the number of electrons accepted by dislocations. The electrons in 
dislocation acceptors cannot carry current in a direction normal to the 
dislocations (and perhaps not in any direction). 

(2) Conduction electrons are scattered by the space charge region 
surrounding the dislocation. Thus the mean free time is reduced. 

(3) The current streamlines are distorted by the space charge cylinders ; 
the carriers cannot drift in straight lines parallel to the applied field ; 
they have to follow curve paths that wind between the cylinders. Even 
when the distance between dislocations is large compared to a mean free 
path (so that the dislocations have a negligible effect on mean free time), 
the dislocations may have a pronounced effect on apparent, or measured, 
mobility because of the distorted streamlines. 

Sections 3 through 5 deal with the case where (2) above can be neglected. 
This will usually be the case in germanium at room temperature— 
see Pearson, Read and Morin (1954). Here the problem is one of flow 
around obstacles and a continuum rather than a particle approach is 
appropriate. Sections 6 through 12 deal with the effect of dislocations 
on the mean free times for both drift and Hall mobility ; the theory is 
worked out by extending Chapter 11 of Shockley (1950). 


Eee 
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§ 2. DEFLECTION oF ELECTRONS BY THE SPACE CHARGE CYLINDER 


The electrostatic potential around a dislocation was discussed in 
the previous paper (Read 1954). The electrostatic energy —qi of an 
electron increases in going toward the dislocation. Consider an electron 
that is moving in a straight line toward the dislocation. When it enters 
the space charge region it experiences a field opposing its motion. It will 
penetrate the space charge cylinder to a distance where —gq7/ is equal to the 
initial kinetic energy ~kT' of the electron. As in the previous paper 
and in the theory of the abrupt pn junction, we shall assume that the 
space charge region has a sharp boundary. This is a good approximation 
when the potential % rises by k7'/q in a distance that is small compared 
to the radius R of the space charge cylinder, as will be the case at low 
enough temperature. Thus an electron that strikes the cylinder head-on 
will penetrate only a negligible distance (compared to R). Electrons of 
glancing incidence will penetrate even shorter distances. We can there- 
fore regard the electrons as making elastic collisions with rigid cylinders. 

As discussed in the previous paper (Read 1954) the radius A of the 
cylinders is related to the spacing c between neighbouring acceptors on 
the dislocation, the spacing a between accepted electrons, and the excess 
concentration of chemical donors V,—N, by the condition of electrical 
neutrality. The number of electrons per unit length along the dislocation 
is lla=f/c. This must be equal to the number 7h?(N;—N,,) of excess 
donors per unit length of the space charge cylinder. Thus 


Pisin Bec(NGae Na) e Gorge | 2 tee a 


In the previous paper we discussed a particular type of dislocation in 
germanium (shown in fig. 1(b) of that paper), for which c=4A. In 1-7 
ohm em n-type germanium (V,—N,= 1015 cem-3), f was found to be of 
the order of 1/10 in the temperature range 0<7'<150°. This was based 
on a dislocation acceptor energy level 0-225 ev below the conduction 
band—a value which is consistent with the experiments. For these 
values R=28304, which is large compared to the spacing a=cjf=404 
between electrons on the dislocation. Tt is, however, only about 3 times 
the spacing (V,—N,)* between excess chemical donors. Thus the 
boundary of the space charge cylinder will not be smooth. The roughness 
of the boundary will give some scattering of electron momentum parallel 
to the dislocations. As discussed in the previous paper such scattering 
can be measured by measuring conductivity and Hall angle with the 
current parallel to the dislocations. The scattering (by dislocations) 
of momentum normal to the dislocations will be much greater. In what 
follows we shall assume that it 1s enough greater than dislocation scattering 
of momentum parallel to the dislocations can. be neglected. 

Let the dislocations run in the y direction, and let N be the number 
per em? cutting the xz plane. If the specimen 1s beat, into a us of 
5 em. as in the experiments of Pearson, Read and Morin (1954), the 


: : i der of 5X 10cm. The exact value 
density N of dislocations 1s of the order of 5 x 10° em™. 12 EXAC 
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will depend on the orientation of the slip systems—see, for example 
Nye (1953) or Read (1953, p. 38-39). 


§ 3. Errrcr or DistocaTIons ON CURRENT STREAMLINES 


In this section we consider the case where the mean free path is small 
compared to R. The reduction in conductivity then comes from the 
effects (1) and (3) discussed in §1. We consider these in order. 

Carrier Concentration. The electrons in dislocation acceptors cannot 
contribute to conductivity normal to the dislocations. Let n be the 
electron concentration in the normal (electrically neutral) n-type material 
and let «=7R?N be the fraction of the volume occupied by the space 
charge cylinders. Then the average electron concentration <n) is given 
by 


* <«n)=n(1—e). fs eee a es 
From the definition «=7h?N and the neutrality condition 


TR (N Ne ) =f/c 
we have 


Nf N, oe 
mgy-Nj (wf OF 


where V.—WN/c is the number of dislocation acceptors per cm®. Thus, 
in a sample of known impurity concentration and with a known array of 
dislocations (which can be introduced by controlled bending and observed 
by the etch pit technique of Vogel, Pfann, Corey and Thomas), the fraction 
f of acceptors that are full can be found from «. We now consider the 
measurement of «. 

Distortion of Streamlines. The <n) conduction electrons per cm? 
cannot drift in the direction of the applied field; they have to follow 
curved paths that wind around the space charge cylinders. Such 
distortion of the streamlines reduces the apparent mobility ; that is, 
it reduces the measured conductivity by more than the ratio (n/n. 

Let the applied field be in the « direction. Then the average field and 
current will also be in the x direction. The local field and current will 
vary throughout the specimen ; the current vanishing in the space charge 
cylinders where n=0. The local field E and current density I are related 
by I=qunE where is the drift mobility of electrons in the n-type material, 
(We assume that hole concentration is negligible compared to Nn.) 

Let ¢) refer to an average over the entire specimen. Then 
(1, )=qu(nE ,), where n=0 in the space charge cylinders (which occupy 
a fraction ¢« of the total volume). Thus 


(L,)=quen(l—e)(E,) ,=qu(n) (EY » 


where ¢),, refers to an average just over the normal n-type material. 
The measured current is (7,,) and the applied field is (2 2) Thus the 
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measured conductivity o*=(J,,)/{E,,) is related to the conductivity 
o=qun of the n-type material by 


Att) A rs eee TEES) 


where g(e) is a factor that measures the relative contribution of the field 
in the neutral n-type material to the average field. For example, if 
the space charge cylinders were touching (so as to form a close packed, 
hexagonal array) then no current could flow and all of the voltage drop 
would occur in the space charge cylinders and none in the n-type material. 
By definition 

=. (E,) rn 

(Hy) 

As «+0, g(«) approaches unity. When the space charge cylinders get 
closer together, so that « increases, then g(¢) decreases. Finally when the 
space charge cylinders are touching g(<)=0. 

S. P. Morgan (private communication) has shown that for small «, 
g(«)=1—e. Thus the apparent mobility o*/¢(n) and the average electron 
concentration (7) decreases with « at the same rate. For larger values 
of «, Morgan has found g(¢) by cutting circular holes in strips of resistance 
paper and measuring the conductance. He has shown that, if the holes 
form a hexagonal array, the conductivity is isotropic. Figure 1 is a plot 
of his experimental results for the hexagonal array.t 

If we can measure (1) the conductivity o* of a sample into which 
dislocations have been introduced by controlled plastic deformation, 
and (2) the conductivity o of the sample before deformation (or the 
conductivity of a control), then from the ratio o*/c, eqn. (3.3), and the 
curve of fig. 1 we have «. From « and eqn. (3.2) we find f. Thus the 
f vs. T' curve can be plotted and compared with the theoretical curves 
to test the theory and determine accurately the energy level of the 
dislocation acceptors, as discussed in the previous paper. 


q(«) (3.4) 


§4. Magnetic Frenp PARALLEL TO THE DISLOCATIONS 


In this section we consider the same case as in the preceding section 
except that, in addition to the electric field H,, we have a magnetic 
field H, parallel to the dislocations and a Hall field #, normal to the 
dislocations and of such magnitude that no current flows in the z direction. 
The current density I at any point is then related to E and H by 


I=qgnE—“2IxH . 2. 1... (41) 


oe 

+ The polar plot of conductivity in the xz plane is an ellipse. Therefore the 
conductivity is isotropic in the xz plane if the parallel dislocations form a regular 
array having two or more fold symmetry. The g(«) plot however would be 
different for different symmetrical arrays. 
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where jy is the Hall mobility and c is the velocity of light. (Not to be 
confused with the c which is the spacing between acceptors along the 
dislocation ; the two c’s never occur in the same formula and the velocity 
of light occurs only in the quotient H/c.) When the mean free time is 
determined by thermal scattering, then the drift mobility «4 and the Hall 
mobility yz are simply related by 


oe 8 
L= Pee . . . ° . . e 


See, for example, Shockley (1950), p. 278. Shockley (p. 287) also treats 
the more general case where both thermal and impurity scattering affect 
the mean free time. 


Fig. 1 


(0) 0.2 Oo4 0,6 0.8 1.0 


Variation of the dimensionless factor g and (1—e)g with the fraction ¢ of space 
occupied by the space charge cylinders when the cylinders form a 
hexagonal array. 


We now average eqn. (4.1) over the entire volume of the specimen. 
From the a components (/,)=qu<(nH,) as in the preceding section 
and we again have eqn. (3.3). The z components in eqn. (4.1) give 
Py ,/c(1,,.=qe<nk .). The x and z directions are both normal to the 
dislocations. Unless there is some marked anisotropy in the arrangement 
of dislocations we can take (nE ,)/(E.)=(nE,,)/(E,)=(n)g(e). This 
will be a good approximation if the (parallel) dislocations are arranged 
either randomly or in some regular array having two or more fold sym- 
metry (so that conductivity is isotropic in the xz plane). We thus have 
MyHl ,/c¢1,.)=quin)(#,)g(e). When we divide this by 


Lo =qe Kn (Ey y¢) 


EEO 
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to obtain the measured Hall angle ¢= (H,)/(E,), the factor g(e) drops 
out and we have the usual formula | 


abl As 


: (4.3) 


Thus we reach the conclusion that an array of parallel dislocations 
should have no effect on the measured Hall angle when the electric fields 
are normal to the dislocations and the magnetic field is parallel to the 
dislocations. This provides a sensitive test to determine whether the 
dislocations in a bent sample are really parallel. 


§5. Macnetic Frntp NorMAL TO THE DisLOCATIONS 


Consider the situation discussed in §4 but let the magnetic field be 
in the z direction (normal to the dislocations) and the Hall field, in the 
direction (parallel to the dislocations). In this section we shall assume 
that the electrons in the dislocation acceptors give a negligible contribution 
to the conductivity in the y direction. As discussed in the previous 
paper, this question remains to be investigated. We would expect 
that an electron moving from acceptor site to acceptor site along the 
dislocation would have a considerably lower mobility than a conduction 
electron in normal n-type material. Thus, if the density {NV , of electrons 
in dislocation acceptors is no larger than the density of conduction 
electrons, we would expect the conductivity parallel to the dislocations 
to be due mainly to conduction electrons in the n-type material. In the 
bound range of temperature, where n<NV,—N,, the dislocations might 
make a significant contribution to conductivity. This remains to be 
studied experimentally by measurements of conductivity and Hall angle 
vs. 7 at low temperature. 

Consider the variation of the electric field E throughout the sample : 
E is constant in the y direction ; therefore the component #, must be 
constant in the x and z directions as well. Then taking averages in (4.1) 
we have, for the y components, 


qexn Lf =H cle ) ce 


while the x components give qu(n){E,)g(«)=(/,,) as before. Thus the 
Hall angle @ is 


pe: oe tg ERD 


The dislocations therefore reduce the Hall angle. Measurements of 
the Hall angle before and after dislocations are introduced should therefore 
give g(e) from which f can be found as described in § 3. 

The measured mobility is o*=qu<n )g(€) independently of the magnetic 
field. If we divide o* by 0 to get (), the factor g/e) dreps out and we 
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have <n) directly in terms of quantities that are known independently 
of the distortion of the streamlines. 


H. ok 
(ny a AE 8 dy ae Sean 


This is the usual formula for average carrier concentration. 

Pearson, Read and Morin (1954) used the standard formulas, eqn. (4.3) 
and uy—cb/H ., to calculate py and <n) with the current and magnetic 
field both normal to the dislocations. Thus the calculated values of 
<n) are correct even though the mobilities are wrong by the factor g(e). 


§ 6. Errect oF DisLocaTions ON MEAN FREE PATH 


In this section we find the mean free path associated with dislocation 
scattering. 

First consider an electron moving with velocity v. Let v, be its 
velocity in the xz plane, or normal to the dislocations. As before, we 
treat the dislocation and its space charge as a smooth cylinder impene- 
trable to electrons. The conduction electrons collide with the dislocation 
and are deflected; v, changes on collision, while the component v, 
remains the same. We now find (1) the collision rate 1/7, for a given 
velocity v and (2) the average change in velocity Av per collision. The 
product of these is the average rate of change of velocity dv/dt of 
electrons having incident velocity v. 

Collision Rate. An electron will collide with any dislocation within 
a distance R of its path. The probability of collision per unit time is 
thus 2Rv, times the number NV of dislocations per cm? cutting the xz 
plane. Thus the collision rate per electron is 


il 
ee Cis: toot eee See re 
Average Change in Velocity. Figure 2 illustrates the reflection of an 
incident electron having velocity v, in the xz plane. The electron comes 
in at a height / above a plane through the dislocation and parallel to v,. 
The change Av in velocity on collision will be equal to the change dv, 
in v,. Suppose that, for any given velocity v, all values of h are equally 
probable. (We discuss this point in more detail later.) Then the 
probability (per unit @) of an electron being deflected through an angle 6, 
is (d/d0)h/R=} cos d=} sin $6,. The average change in velocity { 4v) is 
—v,<(1—cos 0,), or 
ope ih ee Sema ey: a in 40,d6,= = , 
=a Nash ef 1) 58D 264 lestne ahby . (6.2) 
The time rate of change of velocity for electrons having incident 
velocity v is 
ON RE EAS. 
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Since the y component does not change we have, using eqns. (6.1) and 
(6.2) 


dt aia 
where 7, is the mean free time for dislocation scattering and is given by 
hl 0s 02) a, 


Ty, tin G 


(6.3) 


where 
») 


‘= 3RN 


(6.4) 


is the constant mean free path for scattering of momentum in the plane 
normal to the dislocations. 


Fig. 2 


v,+ Av 


Reflection of an electron from an impenetrabie space charge cylinder, 


§ 7. Basic ASSUMPTIONS 


In this section and the following one, we consider the combined effects 
of applied fields and scattering. The analysis will be an extension of 
Chapter 11 of Shockley (1950) ; the same terms are used and some of the 
equations are taken over directly. However, it is necessary to take a 
somewhat more general view from the beginning because, when disloca- 
tion scattering is included, the mean free time is not constant over an 
energy surface in momentum space. ‘a ; ar 

This section discusses the basic assumptions peculiar to dislocation 
scattering. Section 8 will treat the motion of electrons in aaa 
space. Section 9 returns to scattering and uses some pauls of $8 to 
treat the combined effects of dislocation and thermal scattering. The 
following sections consider combined thermal and dislocation scattering 
under various orientations of the electric and magnetic pode eee 
steady state nonequilibrium distribution of electrons in momentum space 


120 W. T. Read, Jr., on the Scattering of 


is found by equating the rates of change due to the fields and scattering. 
The approach is similar to that given by Shockley (1950) except that we 
consider each element of momentum space separately instead of finding 
the rate of change of current in an energy shell. 

We base the analysis on the following two approximations : 

(1) The fields E and H are constant throughout the material. 

(2) The distribution in velocity of the electrons incident on a reflecting 
cylinder is the same as the average distribution over the whole volume of 
material. 

Neither of these assumptions is exact but they become increasingly 
good approximations as the distance between dislocations increases. 
Assumption (1) is most nearly fulfilled when the mean free path for thermal 
scattering is longer than that for dislocation scattering. As thermal 
scattering begins to dominate, the errors in the calculation of dislocation 
scattering are greater; however, they have a smaller effect on the final 
result. The two assumptions are not necessarily exact even when there 
is only one dislocation ; as the following argument will show. 

Consider the distribution of electrons moving in the direction opposite 
to the applied electric field and net current. Let this be the +. direction. 
The electrons are incident on the —w side of the cylinder. Take the 
electron distribution at the point where the incidence is normal, for 
example point A in fig. 2. If the distribution of electrons in velocity 
space is the same at A as in the average for the whole specimen, then there 
will be an excess of electrons going in the +. direction. However, the 
number of incident and reflected electrons are equal at A ; hence there is 
also an excess of electrons going in the —a direction. The reverse is 
true on the -+-x side of the cylinder at the point B ; B will therefore have 
a net electron deficiency. The same general argument applies to electrons 
moving in other directions and leads to the conclusion that the electron 
density is greater on the —w side of the cylinder than on the +2 side 
(the applied field being in the x direction). This nonuniformity of 
electron density constitutes a line dipole, which will induce a compensating 
dipole by displacing the axis of the space charge cylinder in relation to 
the charged dislocation. The effect of the two opposite line dipoles is to 
produce some nonuniformity in the electric field near the cylinder. 


§ 8. Morton or ELectrrons IN MoMENTUM SPACE 


7 


Sections 6 and 7 dealt with the motion of electrons in the semi- 
conducting material. Now we deal with the motion of electrons in 
momentum space and the rate of change of electron distribution under 
electric and magnetic fields. The symbols used and some of the results 
are from Chapter 11 of Shockley (1950). However, the analysis is 
somewhat more general and is based on an approach suggested by 


C. Herring. As in Shockley’s treatment we assume the energy surfaces 
are spherical. 


——— 
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Consider a point P in momentum space. Let f(P) be the fraction of 
the energy states at P that are full at thermal equilibrium. Let f’(P) be 
the fraction that are full in the nonequilibrium conditions, that is, when 
a current is flowing. When there is an electric field E, then J =h 
where /; is proportional to H. Shockley treats the current (in an energy 
shell) associated with f,, then adds the magnetic field H and finds the 
Hall current (current normal to H and E) from the current due to E. 
We cannot follow Shockley’s procedure, however, since. it requires that 
the mean free time 7 be constant over an energy surface, where the 
magnitude v of the velocity is constant. For dislocation scattering 
7T1=7,(V;). Hence, the mean free time for combined dislocation and 
thermal scattering (which we shall derive in §9) is a function of both 
v and v,. It is therefore constant, not over a sphere, but rather over 
a circle around the P, axis. 

We begin with the equation of continuity in momentum space 

Tg (e+"X") Safle ea Meee mee 6!) 

It will be understood that all time derivatives in this section refer to 
changes due to the fields only. Later we add these rates and the corres- 
ponding rates due to scattering and set the sum equal to zero. The 
resulting equations give f’=/’(P). from which we find the currents for 
given electric and magnetic fields. 

The general approach suggested by Herring is the following: /’ is a 
function of EandH. We expand /’ in a series in powers of the magnitudes 
E and H. The terms in H alone vanish (H simply makes electrons flow 
around in circles and therefore leaves the distribution unchanged). 
Since we are interested in measuring conductivity and Hall angle. which 
are linear in the fields, we neglect second and higher powers of # and /. 


=> Tepe) eo has eg See nor 
where f, is proportional to Z, as in Shockley’s treatment, and /, is 
proportional to #H. 

The next step is to substitute (8.2) into (8.1). First however consider 
the meaning of the derivatives E.Vp and vXH.Vvp. The first is 
simply proportional to the derivative in the direction of the electric 
field; v xH.v> is the derivative in a direction tangent to the constant 
energy surface, which we shall call S, and in a plane normal to H. Thus 
it is the derivative along a circle on the sphere S and around the direction 
of H. The equilibrium distribution / (Py es/( é )=f(v) is constant over the 
spherical energy surface. Using this result in (8.1) and separating the 
terms in # from the terms in HH and neglecting terms in H? and HH” 


we have af 
~1— gE.V pf eee eee 8.3) 
and 
fe _1y xH Vines 5 : ; < P (8.4) 
C 
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In the following section we find f, from (8.3) together with the rate of 
change of /', due to scattering. When /, is known /, is found from (8.4) 
and the formulas for scattering. To conclude this section we simply 
state the following results which are discussed more fully in § 11.4 of 
Shockley (1950). 

1— . Vv 4 
Vel = 74 Vp - A= a Bere irs AS), 


where the approximation holds for non-degenerate semiconductors, 
for which f<1 for electrons in the conduction band. 
We now have for (8.3) 


of gE .v ql f 


Thus over an energy surface (where f is constant) the rate of change of 
J, is proportional to —v,. 

We shall show in the next section that /, is also proportional to v,. 
This result, however, does not necessarily extend to f,. We shall see 
in §12 that, when H is in the z direction, so that f, is associated with 
a flow in the y direction, then f, (on an energy surface) is not simply 
proportional to v,. 


§9. THERMAL AND DISLOCATION SCATTERING 


In this section we consider combined dislocation and thermal scattering 
with an electric field Z,, in the x direction, that is, normal to the disloca- 
tions. We shall show that 


qh 7h 
h=— eo, (9.1) 
where the mean free time 7 is given by 
1 1 1 
-=—+— ete ee ee 


T TY Lat 


where 7, and 7, are the mean free times for thermal and dislocation 
scattering respectively when each acts alone. That reciprocal mean 
free times are additive is neither obvious nor generally true. The proof 
for the present case follows. 

Again we use the notation and basic relations of Chapter 11 of Shockley 
(1950). W,; is the matrix element giving the probability that an electron 
at P; in momentum space will make a transition to another point P, 
on the same energy surface S when all the states at P; are empty. The 
rate of scattering from P, to P; is then proportional to Joli) Wee 
The net rate of scattering from P; to P; is 


(Aff AW = (fi Py p= (ffs) Wy 
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(since f is constant over §).* In the cases that Shockley treats W,; 
is a function only of the angle 6 between P; and P;. For thermal 
Scattering alone W’,; is a constant W, The mean free time T, for thermal 
scattering is then equal to the probability of a collision ; so 


cee W ,dS=W,S ee te 9:3) 
ie OS 
where S is the area of the spherical energy surface in momentum space. 
For dislocation scattering W,; vanishes except when P, and P, lie on 
the same circle C around the P,, axis (since the y component of momentum 
is not changed in a collision). When P, and P, do lie on the same circle 
C, then W,,, for dislocation scattering is equal to a function—call it Wi 
of the angle @, between the projections of P, and P; in the vz plane. This 
is the same 6, used in §6 and illustrated in fig. 2. By definition the 
collision time 7,, eqn. (6.1), is related to W,(@,) by 


; 
—= | W1(@,)aC. ft eo ee 
To CY 


W, is a probability per unit length in momentum space: W, was 
a probability per unit area. We found in §6 that W,(@) was proportional 
to 3 sin 36,. Kquation (9.4) determines the constant of proportionality. 
The result is W,=1/4R7, sin 30,. The mean free time 7, for dislocation 
scattering is defined in terms of W, by 

a = as (1==cos 07) W (0, dC =a a en (923) 
Ua Te O- 
which gives 7,= 37,. 

From these definitions we now go on to prove eqns. (9.1) and (9.2). 

The W’s for various scattering mechanisms are additive. In the 
present special case this will lead to the result that the reciprocal mean 
free times are additive. 

Consider the rate of change of f’(P). In $8 we used 0/d¢ for the time 
rate of change due to the electric and magnetic fields. We shall continue 
this notation. To distinguish acceleration and scattering terms we use 
a dot to indicate the time rate of change due to scattering. Then f’=/; is 


fi(P)= |OA(P)—f(P)I Was" _| EAP) ft (6,)0C" - . (9.6) 
S 


The integrals are taken over S and C with P’ as the variable of integration. 
Equation (8.6) gives 0f,/0t, which is due to the field # 2 setting the sum 
of the two rates equal to zero gives an integral equation for /,(P) in terms 


* The same result would have been obtained by neglecting the (1—/) factors. 
In other words, we make no error in the final result if we assume that - 
probability of a transition from one state to another is (a) proportional “ t ie 
probability of occupation of the state that the electron is spite and (0) 
independent of the probability of occupation of the state = ee estan 
is entering. (See Shockley 1950, p. 256.) We shall use this result in § 12. 
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of known quantities. We now show that (9.1) is a solution to the integral 
equation. 

The quantities 7 and f in (9.1) are constant over the circle C and over 
every circle around the y axis. Thus /, is an odd function of v,. and 
therefore of P,,. So the first integral in (9.6) becomes 

i(P) 
are 
The second integral involves only the w and z components of momentum 
and velocity. We can express everything in terms of the projection 
of the velocity in the zz plane. From simple trigonometry we have 
v,/=—v,(1—cos 6,)+v,sin@,. Thus substituting (9.1) into the second 
integral in (9.6) gives 


gH af, v, | 1-cos 0) W(@,)dC +0. [sin 6, @,)d0 J. 
[nf ag C- 


—W,Sf,(P)=—' 


Since W(6,)=W(—8,), the coefficient of v, vanishes. Then using (9.1) 
and (9.5) we are left with —/,(P)/7, for the second integral in (9.6) 


so (9.6) becomes 
: 1 i! : 
fi==h (=+=) ge ee aay 


Adding the scattering rate (9.7) to the acceleration rate (8.6) and setting 
the sum equal to zero gives (9.1). 

If there were impurity as well as dislocation and thermal scattering 
then the reciprocal mean free times would not be additive. 


3 


§10. Drirr Mosiniry 


As in the preceding section let there be an electric field 2, which gives 
rise to a current J,. Suppose J, is scattered by both dislocations and 
thermal vibrations. In this section we find the drift mobility « in terms 
of the mean free path /, for thermal scattering and the mean free path J, 
for dislocation scattering ;/, is given by (6.4). The mean free path for 
thermal scattering is derived by Shockley (1950). It can also be deter- 
mined experimentally as a function of temperature by measurements on 
samples containing a small density of dislocations and impurities. 

Consider an element of momentum space having a number ds of states. 
Let dl be the element of current density due to these states. The net 
current arises from f,=/’—f. Thus dl=qvf,ds. From (9.1) 

2 
di = a v,H Voids. 

Now dl is expressed in terms of the equilibrium distribution jf. The 
equilibrium number of electrons in the ds states is dn= fds. Integrating 
over all of momentum or velocity space gives the total current density 


Glew, 
L= re | atrdn. ae Oe er clien 
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For spherical energy surfaces the other components of current vanish. 
The drift mobility .=TJ ,/qnE , is found from (10.1) : 


é q Je at dn 
ay a ee ey aes « (LOS2, 
iT Tin oe 
It is convenient to use spherical coordinates in velocity space with 
V,;=0 sin 0, v,=v sin 0 cos ¢, | 
(10.3) 
V,= cos 0, V,=V sin 0 sin ¢. | 


The mean free times are related to the mean free paths by 7,=J,/v 
and 7,=/,/v,. The paths /, and /, are constant over velocity space at 
a given temperature. Thus 

i 


paw eh (10.4) 
‘ al ri sin 6 
nT 
The number dn of electrons in an element of velocity space is 
dn=v" sin 6 exp (—mv?/2kT)dé d¢ dv. 
It is convenient to introduce the average thermal velocity 
ly dn QE 
(oy= | =. ike (10.5) 
[ 46 7m 
Now (10.2) becomes 
_ @ &v) (7? sin? 0d0 ix 
a aby sme TRL oe) 


i’ |, 
When there is only thermal scattering /;=0o and the integral equals 
21,/3. For dislocation scattering /,= 00 and the integral is 7/41,. It is 
therefore convenient to define a mean free path /,, where 


[ae eas 
and set 

qd 1 9 
L= TF (vl, Fy pee Cm idl rere aL OLS) 

where F, is a function of the ratio /,//; and is given by 
F bs sige s+ 2) (ee. eee er. 9) 

AE Dat aa taal) 6 
ia + sin 
t 


The left hand side of fig. 3 shows 7’, vs. 1, /l, for by. The right hand 
side shows F’,vs. 1,/1, for 1;<J,. As insured by the definition, 
F,(0)=F,(«0)=1. 

0-99 and 1 it is sufficient for practical purposes 


Since F, is always between Lee 
Tn other words the reciprocal of the 


to use the approximation f,=1. 
Suhe7, VOL..40, NO. 373.—FEB. 1955 


K 
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resultant mean free path is simply a weighted average of the reciprocal 
mean free paths for dislocation and thermal scattering individually. 


§11. Macnetic Fretp PARALLEL TO THE DISLOCATIONS 


Now that we have f, we can proceed to find f,. Since f, is uniquely 
determined by the field #,, it is independent of the magnetic field H. 
In this section we take H=H_, in the direction parallel to the dislocations. 
The procedure will be as follows : first we find df,/0t due to the fields by 
substituting eqn. (9.1) for f, into eqn. (8.3). It will turn out that of,/ot 
due to the fields has the same form as df,/0t; so the same argument 
(§ 9) that gave /, can be used to find f,. We then evaluate the current J, 
associated with f,. From this we find the Hall angle and Hall mobility. 

Substituting f, from (9.1) into (8.3) gives 


Of» gi, Hi 
~~ =— —. on eek 
Ot kT VEX Cc V p(zfvz) ( ) 
Fig. 3 
100 [Kete) 
F 
98 93 
96 96 
94 94 
0.5 10 1.0 05 ia) 
PU cer = bf 0, 


Variation of the dimensionless factors F, and F, with the ratios i,/l, and L,/i,, 
where /, is the mean free time for thermal scattering and 1,, for 
dislocation scattering. 


As discussed in §8 in the paragraph preceding eqn. (8.3), the operator 
v XH .v> takes the derivative along the circumference of a circle on the 
energy surface and around the axis parallel to H. In this section H 
is parallel to the y axis ; so the circle in question is the circle C on which 
the mean free time 7 is constant. Then, since f=f(v) is also constant on 
CG; 
V XH .Vp(zfv,)=7f/v XH .vpv,. 


The vector Vpv, is in the x direction and has the magnitude 1/m where 
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m is the mass of an electron (which is constant for spherical energy 


surfaces). Hence v XH. Vpv,=—v,H,/m and (11.1) becomes 
Pe 
al ie Je. (11.2) 


where w=qH/me is the cyclotron frequency (frequency with which 
electrons move around in circles in velocity space when H acts alone). 

Equation (11.2) shows that df,/ét due to the fields is proportional to 
v,. The mean free time in the z direction is the same as that in the 
« direction (unless there is some marked anisotropy in the arrangement 
of dislocations). Thus the arguments applied to f, in §9 apply also to f, 
and we have due to scattering 


i=. Se ee es Sea: 


- Adding the two rates in (11.3) and (11.4) and setting the sum equal to 
ZeTO gives 


p= ort free. Ce Raa FP tabKey 
The current associated with /, is the integral of qvf, over all the states 
in velocity space. All the components of I, vanish except 


wP?E ; 
Hege 


I= rv 2 dn. eg Pn 115) 


Again using spherical coordinates (10.3), we have 


a n (7 sin? 0 a0 
=F 5 | (11.6) 


on AL sing \ =e 

(cay) 

For thermal scattering alone, the integral equals 3/,. For dislocation 
scattering alone it equals /,. 

We will usually be interested in the open circuit Hall voltage. Let F, 
be the applied Hall field required to give zero current in the z direction. 
We are assuming the x and z directions are equivalent as far as scattering 
is concerned ; so J.=qunH ., where the drift mobility u is given by (10.7) 
and I,, by (11.6). The measured Hall angle is 0=#./H,. The Hall 
mobility p,, is defined by 

cé cH, 
be ve Lik i H yh ' 

In the present case #,/E,=I ,/I,=1 ,/qmpk a Thus: {44 18 propor- 
tional to the integral in (11.6) divided by the integral in ( 10.6). The 
quotient of these integrals is 1, when there is only thermal scattering and 
(4/7)l, when there is only dislocation scattering. Thus it is convenient 
to define a mean free path /;, where 


(11.7) 


ee ae (11.8) 
liz L, 4l, 
Then rata ae, eee, «=| (L1,9) 
[dN vy 2 
K2 
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where F, is a dimensionless quantity equal to the integral in (11.6) 
divided by 1, times the integral in (10.6). #,—1 when either dislocation 
or thermal scattering acts alone. Figure 3 also shows the variation of 
F,, with the ratios 1,/1, and L,/l,. F, varies from a maximum of unity 
to a minimum of 0-956. 

The Hall mean free time 


g n ; 
Sees oN Ne, eee 
w q : 
18 NOW 
ly 
Ty = > ; tint We Meteor. | ds il.tt 
ee is ( ) 


§12. Magnetic FreELp NORMAL TO THE DISLOCATIONS . 

In this section we take the magnetic field in the z direction, that is, 
normal to the dislocations. The electric field is still in the x direction. 
Again we have equation (11.1). Now v xXH.vV> is the derivative along 
a circle about the z axis. The mean free time 7 is constant on a circle 
about the y axis. Thus the results will involve the gradient of + in 
momentum space. Equation ae becomes 


lame _ 
ot 


where v XH. Vpv,=v,H./m. Now consider the second term in the 
brackets : 


H 
TV x “Vpvz ees Vv ey vor | . . (12.1) 


Vv XH .V pTt=— TV xH ve(7 +7) 
t 


a ee xH v.= TH Use; 
uote a lmv, * 
Thus (12.1) becomes 
Of, qu, pase 
3a (*- FE) Mots ae (12.2a) 


where w=qH ./mc. Clearly f, gives a current J,=J, in the y direction. 
Since f, is not constant on a circle about the y axis, we cannot take 
fo=—fe/7, Instead let us consider the change in J, due to changing 
occupation of the states along two equal circles about the y axis in 
velocity space. The circles have equal values of v, v, and 7 but opposite 
values of v,. Call the two circles together C. The fields cause a rate of 
change of current proportional to (v,(0f,/0t)), where the average is taken 
over the two circles and simply involves inveprating with respect to ¢. 
From (12.2a), and (v,?),=v,?2/2 we have 


Of» qok , Ty 
€ y “Ot \ = kT (-3) Vy a (12.25) 


Scattering has the following effects on occupation of the states along 
the circles C’: 
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(1) Dislocation scattering causes a redistribution of electrons on each 
circle. 

(2) Electrons are thermally scattered out of the states along the circles. 

(3) Electrons are thermally scattered into the circles. 

Clearly (1) has no effect on J,. As discussed in §9 (footnote) we can 
assume that (a) the probability ane an electron will be scattered out of 
a group of states is proportional to the fraction of those states that are 
full and (6) the probability that an electron will be scattered into a group 
of states is independent of the occupation of those states. If electrons 
are equally likely to be thermally scattered into any point on the circles, 
then (3) gives a change in current proportional to (v),=0. The only 
non-vanishing contribution comes from (2) and is proportional to 
—lj7,{v,f2).._ Setting this equal and opposite to (12.26) we have 


defo) = FE (1 Rye ape P8183) 


Expressing this in spherical coordinates and integrating q(v,,f.), over all 
states in velocity space gives 


gH, om tile (12.4) 


hair kT 
where the length /, is defined by 
amt/2 20 om/2 Rf) 
ess cos?6 sin Od@ 1 | sin? dé (12.5) 


Peel) sain 00 9h) 212 a (leeeeoinne yon 
Oa (+) 

When there is only thermal scattering (/,= ©), then /,=1,/3. As 1,/1, 
approaches zero /, approaches zero. 

There are two ways of defining the Hall mobility uy=cb/H,. First 

we can measure the short circuit Hall current J, and define @ as the ratio 
I,|/I,.. Secondly we can apply a Hall field #, so that no current flows 
in the y direction and define 6 as the ratio # jE. In the present case, 
where the x and y directions are not equivalent, these two definitions of @ 
will not give the same result. We shall define 6=4,/H, because this is 
what is usually measured. The scattering in the y direction is entirely 
thermal. Therefore /,=qunl, where py=q(v)l/skT is the drift 
mobility for thermal scattering. Setting the J, due to #, equal to the 
Zin eqn. (12.4) we have 


ch — 3ql, 
BEB ROY | 
It is convenient to express the results by the ratio of uz, to the Hall 
mobility py =8ey for thermal scattering alone. Then, using (10.5) 


(12.6) 


we have 2 
MH MS e O ® Fy A : ° (12.7) 


ut Ly 
Figure 4 is a plot of wg/ea: Vs. 14/1; for 1, <1, and vs. 1,/l, for 1;<4,. 
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§ 13. CONCLUSIONS 


This section reviews the range of validity of the results derived in the 
preceding sections and considers possible experiments and further 
theoretical analysis. ot 

As discussed in § 7, the analysis of mean free time will be most accurate 
when the dislocations are farthest apart and the fraction « of space 
occupied by the space charge cylinders is small. As « increases, assump- 
tions (1) and (2) in §7 become increasingly poor approximations. Another 
effect that limits the validity of the results at large « is the following : 
according to the derivation in §9 a probability can be assigned to any 
given sequence of transitions, or deflections. For any such sequence we 
can find the positions of the deflecting cylinders. We have assumed, 
however, that the space charge cylinders do not overlap. Thus certain 
sequences of transitions are forbidden because they would require over- 
lapping cylinders. This effect is probably not important except when 
the average minimum spacing between cylinders becomes less than R. 


Fig. 4 
Lada by 
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Variation of the ratio of Hall mobility .,, for combined dislocation and thermal 
scattering to Hall mobility ,.,,, for thermal scattering alone. The electric 
and magnetic fields are both normal to the dislocations. 


To study the effect of dislocations on mean free time it is desirable 
not only to have « small but also to have 1/1, at least comparable to 1/I,. 
This requires R to be small since J, is proportional to R and « to R. 
However, to make & small it is necessary to have a large excess donor 
concentration N;—N,, which means that impurity scattering becomes 
important. Because the effect of dislocations on mean free time is most 
pronounced at low temperatures, it would be useful to work out the theory 
of combined dislocation and impurity scattering. However, in any case 
the formulas of §§8 to 12 with 1, will be valid when dislocation 
scattering is the dominant scattering mechanism. vt 


Electrons by Charged Dislocations in Semiconductors 131 


Over most of the temperature range where thermal scattering is more 
important than impurity scattering we have to consider both the distortion 
‘of the streamlines and the reduction in mean free time. As an approxi- 
mation we may use the formulas of §§3 to 5 with the drift and Hall 
mobilities given by the formulas in §§8 to 12. Then, to find the effect 
of dislocations on mean free time, it would be best to have the magnetic 
field along the dislocations so that the factor g(«) cancels. For the same 
reason to find (n) we would want the magnetic field normal to the — 
dislocations. 
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ABSTRACT 


Single phase ternary alloys of Mn, Zn and C have been prepared. 
These alloys have an ordered structure of the perovskite type and show a 
spontaneous magnetization at room temperature. The observed magnetic 
properties are exceptional and indicate a new type of ferrimagnetic 
behaviour postulated by Néel but hitherto unobserved in experiment. 


$1. INTRODUCTION 


Iy the course of work on binary alloys of manganese and aluminium 
G. R. Piercy noted the occurrence of magnetic properties in certain alloys. 
It was found that these alloys contained carbon as impurity. Sub- 
sequently further work identified the magnetic phase as a ternary alloy 
of manganese aluminium and carbon possessing a face centred cubic 
structure, and single phase ternary alloys were obtained over a wide range 
of composition. Further work by Morgan and Butters (Morgan 1954) 
showed that magnetic alloys of similar strueture could be produced when 
the aluminium was replaced by indium, tin or zine. 

The magnetic properties of these alloys are being investigated and this 
paper deals with the results obtained for the manganese zine carbon 
system, these alloys having particularly interesting magnetic behaviour. 


§ 2. PREPARATION OF THE ALLOYS 


The materials used in the preparation of the alloys for this work were 
manganese of 99-9°, purity donated by the Electromanganese Corporation 
of America, zine dust of 98° purity} and graphite of spectroscopic grade. 

The use of zinc dust was necessary because the alloys were produced 
by sintering the intimately mixed powdered components. In some cases a 
master alloy of manganese and carbon was first prepared by induction 
heating under an argon atmosphere. This alloy was then crushed and 
mixed with the appropriate amount of zinc dust to produce the mixture 
for sintering. The mixture of powders was pressed into a clean fused 
quartz tube and held in place by a quartz plug. The tube was evacuated 


and sealed, and the sintering achieved by heating to 900°c for 72 hours. 
ee eee 
* Communicated by Professor W. Sucksmith. 
+ The major impurity being oxygen. 
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For alloys containing less than about 15 atomic % zine this treatment 
produced hard compact rods of sintered material whose composition was 
close to the proportions of the materials used. 

Tn alloys with zine contents greater than 15 atomic °/ a separation of 
zine invariably occurred, the zinc condensing at points where there was 
free space; this separation of zinc was accompanied by the rejection of 
some graphite from the sinter. To avoid this loss of zinc by evaporation 
a lower sintering temperature of 600°C was also used. The zinc loss was 
thereby reduced but not completely eliminated, and it was not possible 
to produce an alloy of the desired optimum composition. The sinters 
were homogeneous but not so hard or compact as those prepared by heating 
at 900°c. 

The alloys decompose in moist air but not so quickly as the isomorphous 
aluminium alloy which at certain compositions will disintegrate completely 
in a few hours. They were therefore kept in a dessicator. 

Single phase alloys were produced containing varying amounts of zinc ; 
they all showed a spontaneous magnetization at room temperature and 
were thought to be ferromagnetic. A full investigation, however, revealed 
strong evidence for their being classed as ferrimagnetic substances. 


§ 3. STRUCTURE AND LATTICE PARAMETER 


The homogeneity of the alloys was determined by taking Debye— 
Scherrer diffraction photographs of powders selected from different 
portions of the sintered mass. Iron radiation and a manganese filter were 
used. Although the sinters were very hard the material was capable of 
plastic distortion and sharp diffraction lines were observed only with 
annealed powders. 

Within the range of composition studied, C-20 atomic %, Zn—10—20 
atomic %, Mn—70-60 atomic °%, the x-ray diffraction patterns showed the 
structure to be face centred cubic and the presence of extra diffraction 
lines indicated pronounced ordering. It was apparent that the ordering 
would be most complete for an optimum atomic composition Mn,ZnC, 
and the ordered structure has been determined for an alloy close to this 
composition by comparing the intensities of the diffraction lines measured 
by a geiger counter spectrometer, with the calculated intensities expected 
from an assumed structure. 

The diffraction lines typical of a face centred cubic structure were strong 
and this indicated that it was unlikely that carbon occupied a face 
centre or cube corner position. The structure assumed was the perovskite 
type shown in fig. 1. The calculated intensities of the x-ray reflections for 
this structure are given in table 1, together with the measured intensities. 
Assumptions were necessary in these calculations because Fek« radiation 
has a wavelength close to the absorption edge of manganese and the atomic 
scattering factor is reduced. The decrement in the atomic scattering 
factor was taken as —1-70 for Zn, —4-22 for Mn and 0 for carbon. The 


temperature factor was omitted. 
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In comparing the calculated intensities with those observed for a speci- 
men of composition Mngo.gZN 19.429 reasonably good agreement is found 
providing a quantitative justification for the assumed perovskite structure. 
Although diffraction lines of index » h2=1, 5, 9 and 13 were observed 
photographically it was not possible to estimate their intensity values with 
the geiger counter spectrometer. 


Table | 
Line index Intensity 
De Visual estimate Measured Calculated 
1 v. weak not measurable 0-8 
2 weak 11-5 7:8 
3 v. strong 100-0 100-0 
4 strong 64-9 52-5 
5 v. weak not measurable 0-86 
6 weak 3-0 2-63 
8 medium 35-4 33-6 
9 v. weak not measurable 0-5 
10 v. weak 1-1 0-74 
ll strong 45-2 40-5 
2 medium 20-0 17-0 
13 v. weak not measurable 0-5 
14 weak 7:4 d-1 
16 medium 39-0 39-6 


The measured and calculated values of the line intensities are based on a 
value 100 for the strongest line Xh?=3. 
Specimen composition Mng9.9Z49-;Cao- 
Lattice parameter 3-92394. 
Fig. 1 


Proposed structure of Mn,ZnC, 


The measured lattice parameter of an alloy containing 19-1 atomic % Zn 
was 3:923,4 and extrapolation over a composition range gives the value 
3-924, 4 for the optimum composition Mn,ZnC. 
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The ordered structure was extremely stable and occurred when the 
alloys were slowly cooled or quenched from the sintering temperature. 

The phase was stable only for carbon contents close to 20 atomic %, 
increasing the carbon content beyond this value caused the separation 
of graphite from the sinter whereas a mixture containing 17-5 atomic % C 
produced a two-phase alloy on sintering. In all alloys the carbon content 
was therefore maintained at 20 atomic °% and the composition varied by 
changing the ratio of the manganese to zinc contents. On account of the 
simple sintering technique the Mn : Zn ratio could not be made less than 
3: 1 due to the difficulties in preventing zinc loss by evaporation. (It is 
thought, however, that if an improved sintering process made it possible 
to reduce the Mn: Zn content below 3, useful alloys would not result 
because the isomorphous Mn, Al, C system shows two phase alloys when 
the ratio Mn : Al is less than 3, for 20 atomic °4 carbon.) 

A change in composition was produced by increasing the Mn : Zn ratio 
to values greater than 3: 1, and several alloys with varying zinc content 
were prepared but at 8 atomic °% Zn instability occurred, the sinter being 
composed of two coexisting face centred cubic structures with different 
lattice parameters. 

A preliminary determination of the structure of the alloy containing 
19-1°% zinc at low temperatures was kindly made by Dr. W. B. Pearson of 
the National Research Council, Ottawa. A Debye—Scherrer x-ray 
diffraction study was made using the method of Smith and Lonsdale in 
which liquid nitrogen is passed over the specimen. These experiments 
showed that at —196°c the original face centred cubic structure of the 
alloy observed at room temperature is slightly distorted becoming face 
centred tetragonal with a=3-9214, and c/a=0-9947. Superlattice lines 
were also present indicating that the tetragonal structure was also ordered. 


§ 4. MAGNETIC PROPERTIES 


All the alloys studied possessed a spontaneous magnetization at room 
temperature. The saturation intensities of alloy specimens were measured 
on a ring balance* similar to that designed by Sucksmith (1939). This 
apparatus permitted measurements to be made from —196°c to the Curie 
temperature, the high temperature measurements being made in vacuo or a 
reduced atmosphere of argon. Paramagnetic measurements were made on 
a similar but ‘more sensitive ring balance (Sucksmith 1929, Sucksmith 
and Pearce 1938). a 

Coarsely powdered. specimens were found to saturate easily in magnetic 
fields greater than 7000 oersteds and all ferromagnetic measurements 
were subsequently made in a field strength of 16 200 oersteds. 

Table 2 contains the measured Curie temperature and values of satura- 
tion intensity at 20, —45, and —196°c for alloys of different composition. 
The alloy closest to the optimum composition has the lowest Curie 


* We are indebted to the Beryllium Corporation of America for the gift of 
Beryllium copper strip which was used in this apparatus. 
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temperature, the latter parameter increasing progressively with increasing 
Mn: Zn ratio, which corresponds to a decreasing lattice spacing. The 
graph in fig. 2 shows the variation of Curie temperature with lattice 
parameter. 


Table 2 
Composition Mn: Zn Lattice oy =16 200 at Curie 
atom °% atom parameter —196°c —45°c 20°C point 
ratio A ao 
Miia tllere on 7:32 3-899, 61-6 60-4 58-6 488 
Mi, 2092 ay 4-94 3-912, 82-5 79-0 74-0 325 
Mnjs520; 64 Ose 3°86 3-920, 83-6 88-6 77-0 177 
Nitige li yrs ap 3-48 3°921, 771 39-0 76-0 143 
Mn go-570ie4 Can 3:18 3-923, oe 85-5 70-0 107 
MN gp ZNoqQ a9 3-0 3-924, — — — 80 


The values quoted for Mn,ZnC were determined by extrapolation. 
Curie temperatures were obtained using o?—T' graphs. 


Fig. 2 
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The variation of Curie temperature with lattice parameter. 


In fig. 3 the temperature variation of the saturation intensity is shown. 
The most interesting feature is the behaviour at low temperatures. 
For alloys with large zinc content a marked maximum in the magnetization 
occurs in the region of —40 to —50°c, Above this temperature the mag- 
netization decreases in the usual fashion becoming zero at the Curie point. 


i 


Structure and Properties of Ternary Alloys of Mn, Zn and © 137 


Below this temperature region the magnetization decreases but the 
decrease becomes less the lower the temperature; the graph is slightly 


‘convex towards the temperature axis. Unfortunately on account of the 


16 200 oersteds 


Saturation intensity per unit mass H 


lack of appropriate cryostats we were unable to make measurements at 
still lower temperatures and extrapolation to 0°K although possible is 


Fig. 3 
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somewhat uncertain. It is evident however, that the magnetic saturation 
intensity at 0°K will be only slightly lower than that at —196°o, although 
measurements at still lower temperatures are needed to establish accurately 
the manner the magnetization varies as the absolute zero of temperature is 
approached. 

This maximum in magnetization was observed in several alloys of 
different composition. As the composition was varied so that Mn: Zn 
atomic ratios became progressively greater than 3 the decrease at low 
temperatures became less marked until for alloys of about 15 and 10 
atomic % Zn the maximum disappeared and the magnetization became 
constant at low temperatures. 
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Usually the low temperature measurements were made as the specimen 
warmed after being immersed in coolant, the temperature being measured 
by a copper constantan thermocouple ; for reasons expressed in a later 
section measurements were also made as the specimen cooled to low 
temperatures and we were able to show that the same variation was 
obtained. It is important to stress the fact that the variation of the 
magnetization in the region of —45°c was completely reversible with 
temperature ; no thermal hysteresis was observed. 

It is interesting to note that, if we extrapolate the low temperature 
magnetization values and attribute the magnetic moment to the manganese 
atoms alone, in an alloy of composition Mn,ZnC the manganese would 
have an atomic moment of approximately one Bohr magneton per atom 
at 0°K. The calculated saturation intensity per unit mass of Mn,ZnC 
with atomic moment lps for the manganese atoms would be 
69 erg/g/oersted. 

Paramagnetic measurements above the Curie temperature were made 
on the alloy having composition nearest to the optimum composition 
Mn,ZnC. Unfortunately only a narrow temperature range was available 
for accurate measurement, for although measurements were made in a 
reduced atmosphere of argon it was impossible to prevent the loss of 
zine by vaporization and an extended range of temperatures could not be 
used without causing serious error ; as it was measurements were made up 
to 650°c. Invariably measurements taken as the specimen cooled gave 
smaller susceptibility values than those obtained as the temperature was 
increased. By comparing values obtained during the heating cycle with 
those obtained when the specimen was cooled quickly we have attempted 
to correct for this zinc loss, but it should be emphasized that the values 
obtained are nevertheless uncertain and can only show the trend of the 
behaviour at high temperatures. Subsequently resistance measurements 
made under a full atmosphere of argon showed that zine distills from the 
high zine content alloys at about 400°c and on this account it is thought 
that accurate susceptibility values will only be obtained by the use of 
positive pressure apparatus. 

In spite of these difficulties the trend of the variation of 1 /x with 7 
could be determined, and all our measurements showed the l/y—T graph 
to have slight but distinct curvature concave to the temperature axis. 
This variation is shown in fig. 4. 

Klectrical Conductivity 

Usually the sintered specimens were quite small (about 2 g) but in 
order to make electrical resistance measurements a sintered rod was 
prepared. A potentiometric method was used and measurements under 
one atmosphere of argon were made up to 550°c, The apparent resistivity 
of the sintered sample was 770 microhm cm at 20°c. The resistance 
varied linearly with temperature, the substance having an apparent 
negative temperature coefficient of 1-1 10-4 per °c. Above 400°c the 
resistance of the specimen began to increase and on cooling this increase 
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in resistance was maintained, although the temperature coefficient of 
resistance below 400°C was almost unchanged. 

On completion of these measurements the copper contacts and the leads 
not in contact with the specimen were found to be covered with a layer 
of bright yellow brass formed through distillation of the zinc. The 
initial composition of the sample for resistance measurements was 
Mng3.gZ0 46.4029, the apparent density and the x-ray density 4:2 and 
6-63 g/cm? respectively. 


Fig. 4 
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§ 5. Discussion 


Although an alloy of exact composition Mn,ZnC was not studied it is 
convenient to discuss certain aspects of the observed magnetic behaviour 
in terms of the structure of the optimum composition ; there are then 
three manganese atoms per unit cell occupying face centred positions. 
Extrapolation of the measurements at low temperatures shows ie 
Mn,ZnC has a magnetic moment of approximately three Bohr magnetons 
per unit cell at 0°K. We 

The interesting feature of the magnetic properties is the aan of 
the saturation intensity which occurs in the region of —45°c. The 
observed variation has two possible explanations, arising from structural 
changes or having a true magnetic origin in ferrimagnetism. The 
possibility of a structural change is considered first. The 
would decrease with temperature below —45°c if a new phase either we 
magnetic or less magnetic than that stable at higher temperatures was 
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formed below —45°c. A phase change at such a low temperature would 
most probably be a diffusionless transformation and be manifested by = 
abrupt rather than a continuous change in the magnetization. a 
type of phase change would also be expected to show pronounced ames 
hysteresis on reversing the temperature change and such hysteresis was 
not detected in the magnetic behaviour. Ps, 

Although x-ray measurements made at liquid nitrogen temperature 
showed the structure to be tetragonal we do not consider that the magnetic 
behaviour is to be attributed to a structural change. The distortion of 
the face centred cubic cell is very slight (c/a=0-9947) and we believe that 
this change of crystal structure may be explained in terms of a ferri- 
magnetic description of the alloy ; that is the structural change arises from 
the underlying magnetic properties of the manganese atoms rather than 
the magnetic properties changing on account of an alteration in structure. 


Fig. 5 


rf a T 


Three types of ferrimagnetic o—7 characteristics postulated by Néel. 


Néel (1948) has shown that the magnetization of a ferrimagnetic 
substance may vary with temperature in a number of ways dependent upon 
the distribution and interaction of the magnetic carriers. In particular 
fig. 5 shows those types of ferrimagnetic behaviour of interest in connection 
with the present work ; parts (a) and (0) of fig. 5 both show a maximum in 
the magnetization and differ only in the manner the magnetization varies 
as the absolute zero of temperature is approached. 

The observed behaviour of Mn,ZnC is in keeping with that of a ferri- 
magnetic substance having a magnetization temperature characteristic 
similar to that of fig. 5(b). For the alloys with Mn: Zn atom ratios of 
+-94 and 7-32 the magnetization no longer shows a maximum but becomes 
constant in the low temperature range studied. This variation is similar 
to the ferrimagnetic characteristic shown in fig. 5 (c). 

That the magnetization temperature characteristics of the alloys 
change from type (b) to type (c) of fig. 5 as the composition changes is in 
agreement with the general results of Néel (loc. cit.) who has shown that 
this is one of the trends to be expected when the composition or the inter- 
molecular field constants are changed. It is possible however that an 
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alloy of exact composition Mn,ZnC may vary in the manner of fig. 5 (a) 
and the type of magnetic characteristic change through type (b) to type (c) 
of fig. 5 because the curvature at low temperatures observed in the 
specimen of composition Mn,,Zn,gC.4) (Mn : Zn ratio 3-2) is slight and in 
this respect measurements at lower temperatures than —196°c are needed 
to establish accurately the low temperature behaviour. 

Ferrimagnetic substances, as shown by Néel and exemplified by the 
studies on ferrites, have a distinctive paramagnetic behaviour above the 
Curie temperature, the susceptibility y is given by 


Dee Ce Wyo 6 4p (1) 


C being the Curie constant, and 1/y 9, o and 6 parameters for a given 
substance. The eqn. (1) represents a hyperbola concave to the temperature 
axis. Although the measured susceptibility values of the alloy 
Mngo.9ZN49.4Co9 are uncertain the trend of the variation is in keeping with 
an equation of type (1), the observed curvature (except in the neighbour- 
hood of the Curie temperature) being concave to the temperature axis and 
not convex as is usual for a ferromagnetic substance. 

If the alloy Mn,ZnC is ferrimagnetic then the two sublattice model 
rigorously investigated by Néel should be applicable as it is for the 
ferrites. For these latter substances, however, their structure allows a 
non arbitrary division of the magnetic ions to form the sublattices and the 
state of the ions is known. At room temperature the manganese atoms in 
Mn,ZnC occupy equivalent lattice sites and their state of ionization is not 
known; the structure at low temperature, however, destroys the equiva- 
lence of the manganese atom sites and below a simple model is described 
which accounts qualitatively for the observed structure and magnetic 
properties. 

On account of the observed magnetic moment per unit cell of Mn,ZnC 
it seems justifiable to consider simple ratios. The simplest division is 
to assume that the manganese atoms on the side centres form one sub- 
_ lattice (A) containing two atoms per unit cell and the end centres form the 
second sublattice (B) possessing one atom per unit cell of Mn,ZnC. 

If, furthermore, integral atomic moments are assumed for the manganese 
atoms (as for instance in the Heusler alloys) the observed magnetic 
moment of Mn,ZnC can be obtained by the arrangements described in 
table 3. For ferrimagnetism the magnetization of the sub-lattices A and 
B are directed oppositely and table 3 gives the states of the ions A and B 
that would account for the observed moment. 

It is unlikely that manganese ions in widely different states would 
coexist in Mn,ZnC and choice d is not considered on this account ; in 
addition the number of free electrons available makes the choice c 
unsuitable unless the carbon atoms contribute electrons to the d shell of 
the manganese. Furthermore the presence of manganese ions in different 
states would distort the structural cell from that of the face centred cubic 
which was accurately observed at room temperature. The most probable 
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state for the manganese ions is therefore that with seven 3d electrons 
providing each ion with an atomic moment of 3 uB. 

It is unfortunate that the volatility of the zinc content in these alloys 
prevents the accurate measurement of the paramagnetic susceptibility 
over a wide temperature range and thereby precludes the determination 
of the intermolecular field coefficient n, x and B (see Néel, loc. cit.). It is 
not possible therefore to calculate the variation of the magnetization with 
temperature as described by Néel. 


Table 3 
State of 3d shell State of 3d shell Resultant moment 
ion A ion B per unit cell 

in Bohr magnetons 
a 3d6 3d° 3 
b 3d7 3d? 3} 
c 3d8 3d° 3 
d 3d° 3d° 3 


Using the ferrimagnetic description outlined above the magnetic 
properties may be described qualitatively as being similar to those postu- 
lated by Néel. The low temperature structure may be explained in terms 
of the attractive and repulsive forces exerted between ions with oppositely 
and similarly directed resultant atomic spins (see Zener 1951). These 
forces will be stronger the lower the temperature, the individual sub- 
lattices possessing greater magnetization the lower the temperature. If 
this is the case the result of these forces will be so as to expand the 
structural cell in the basal plane (taken as the 001 plane in fig. 1) and to 
contract the cell in the direction normal to this plane; if these forces are 
strong enough the original face centred cubic cell will be distorted to 
become tetragonal with axial ratio slightly less than unity. Using this 
description it is to be expected that the axial ratio would change pro- 
gressively as the temperature is lowered ; x-ray measurements to deter- 
mine whether this is in fact the case are projected. 

With regard to the role of the carbon atom occupying the body centre 
position, the variation of lattice parameter with zine content shows the 
zinc atom to have slightly greater atomic diameter than the manganese. 
These diameters are estimated to be 2-75, A and 2-79, A for the manganese 
and zinc atoms respectively. Applying Higg’s empirical rules (see 
Hume-Rothery and Raynor 1954) for the formation of interstitial phases 
in alloys of face centred cubic structure it is to be expected that the 
diameter of the carbon atom should lie between the limits 1:13 and 
1-634. The carbon atom must therefore be presumed bonded to its six 
equidistant manganese neighbours. 

The manganese zine carbon alloys all showed electrical conductivity, 
a particular alloy of composition Mng3.gZn4¢.4Co9 in sintered form 
possessing an apparent resistivity of 770 microhm cm. This value is 
probably much greater than the true resistivity of the alloy by a factor of 
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the order x 10 on account of the porosity of the sintered mass, it may there- 
fore be said to have electrical resistivity typical of a metallic alloy. The 
_ observed negative temperature coefficient of resistance is exceptional but 
this may be due to the porosity of the specimen causing resistance changes 
which mask a true metallic behaviour. 

If we accept that the manganese ions possess seven 3d electrons, the 
arrangement may be obtained by assuming that the two 4s electrons of 
each atom enter the 3d shell when manganese combines to form Mn,ZnC, 
requiring the valency electrons of the zinc atoms to be shared by the 
manganese and zinc atoms providing 0-5s electrons per atom ; this would 
account for the observed conducting properties. 
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Note added in proof.—Recent low temperature x-ray studies made by the 
authors have shown that the structure of the alloy Mn,ZnC changes in the 
temperature region where the magnetization is a maximum. The lattice 
parameters change in a continuous fashion indicative of a second order transition 
occurring in this temperature region. A preliminary neutron diffraction 
investigation kindly made by Dr. B. N. Brockhouse of Atomic Energy of 
Canada, Ltd. provides further evidence to support this point of view. In 
addition these neutron diffraction results indicate that at high temperatures 
the alloy behaves as a ferromagnetic substance but at low temperatures there 
occurs a gradual ordering of the atomic magnetic moments. It would appear 
that the second order transition previously mentioned is to be associated with 
a change from ferromagnetic to ferrimagnetic properties but in order to under- 
stand the detailed behaviour of the alloys further x-ray and neutron diffraction 
studies are needed. This work is in progress and the results will be published 
at a later date. 
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SUMMARY 


Three examples of negative heavy mesons which produce nuclear 
disintegrations at the end of their ranges in emulsion, are described. 
Certain characteristics of the parent and capture stars of these mesons 
are discussed. 


$1. INTRODUCTION 


THREE examples of nuclear disintegrations produced by slow negative 
heavy mesons in photographic emulsions have been found in the course 
of general scanning at Bristol and Dublin. Details of these are given in 
$§2, 3 and 4 respectively. Certain features of interest connected with 
these disintegrations and the stars from which the K -mesons were 
ejected are treated in § 5 in relation to some general characteristics common 
to events of a similar nature previously reported. 

The three events were observed in stacks of Ilford G5 stripped emulsions 
(flights $30 and $35) flown during the International Sardinian Expedition 
in Summer 1953. 


§2. Event K” (Br 1) 


(a) The Primary Particle 

The K-meson has a range of 12-23 mm in five emulsion strips. Mass 

measurements using standard procedures yield the following values : 

(i) By the «R method : My-=750+120 m, ; 
(ii) by the constant sagitta method, using the 
cell-size scheme of Fay, Gottstein and 

Hain (1954) : Mx-=590+130 m, ; 

(iii) by grain density v. range : My- = 840+100 m,. 


(b) The Secondary Disintegration 
On coming to rest in the emulsion the K-meson produced a character- 
isti¢ o-star with three branches, all of which ended in the stack, with 


* H. H. Wills Physical Laboratory, University of Bristol. 
+ Institute for Advanced Studies, Dublin. 
{ Physics Department, University College, Dublin (now at University of 
Rochester, N.Y.). : 
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respective ranges 5, 5920u and 17-27 mm. The first of these is too short 
to be identified ; the second is steep and from its appearance seems 
most likely to be a proton ; the third was produced by a proton with an 
energy at emission of 71-5 Mey. Thus the visible energy in this disinte- 
gration is about 110 Mev. 


(c) The Parent Star of the K -meson 

The parent star is of type 3+-0n, and of the branches, only that due to 
the K-meson can be identified unambiguously. One of the other tracks 
ends in the stack with a range of 3-72 mm and is very steep. The 
remaining particle had a projected track length of about 0-6 mm per plate 
and left the stack after traversing 25-7 mm of emulsion. From the rate 
of change of grain density along the track it appears that the mass of 
this particle is consistent with that of a proton or deuteron ; if a proton, 
its energy at emission was about 128 Mev. On the other hand, if one 
assumes that the track was caused by a charged hyperon, the time of 
flight before the particle left the stack was 1:8 10~1° sec, during which 
no decay occurred. 


$3. Event K” (Br 2) 


(a) The Primary Particle 

The track of the K-meson in this event was of rather poor geometry ; 
it had a range of 563, in two emulsion strips and an average dip of 
about 45°. In spite of this, however, it has been possible by careful 
gap-length measurements to show that the particle has a mass inter- 
mediate between those of the 7-meson and proton. 7 

Figure 1 shows a plot of mean gap-length versus range for the K’ -meson 
and for 8 protons and 1 deuteron which have roughly the same dip, lie 
in the same region of emulsion and were measured in the same layer as 
was the K-meson. The lengths of all resolved gaps were measured on the 
appropriate sections of the tracks and the mean gap-length thereby 
deduced. The measured mean gap-length on the K-meson at a residual 
range of 449 was 0-455-- 0-046 w whereas that expected at this range in 
these emulsions for a 7-meson is 0:655 4 and for a r-meson, 0:395y. It is 
clear then that the particle involved is almost certainly a K-meson, since 
the mean gap-width lies more than 4 standard deviations from that for a 


S210 a5, 
the corresponding 7-meson. 


The mass value deduced was 630775) m-- 
This rather low value could perhaps be due to an inelastic collision along 
its path; in fact there is a sharp deviation of 28° at a distance of 30, 
from its end. 
(b) The Secondary Disintegration 

The capture star produced by the K -meson has four branches, all of 
which end in the emulsion with respective ranges of 2576p, as [hs a6 m 
and 5:8. The first has been identified as a proton, the second is probably 
a deuteron or triton, the third probably a proton or deuteron while the 
nature of the particle producing the short prong is impossible to establish. 
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The kinetic energy of the charged fragments lies between 53 Mev and 
§1 Mey and is probably about 59 Mey. 


(c) The Parent Star of the K -meson 
The parent star is of type 12+0p. In addition to the K-meson, nine 
of the other eleven particles come to the end of their range in the emulsion 
and were identified as protons, deuterons or «-particles. Two particles 
escaped through the edge of the emulsion and were unidentified; the 
tracks in these cases were steep and grey. No evidence was found for 
the emission of another unstable particle from the star. 


Fig. 1 
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Plot of the mean gap-length versus range for the tracks of the K -meson, 
8 protons and | deuteron. 


§4. Event K (Db 1) 
(a) The Primary Particle 
The K-meson had a range of 19-60 mm in ten plates and the ionization 
in the neighbourhood of the interaction was consistent with that of a 
particle of mass ~ 1000 m, which was at rest or moving with negligible 
velocity. ' 


Estimates obtained for the mass of the primary particle by different 
methods were as follows : 


(i) Scattering v. range (constant sagitta) : My-= 780+100 m, ; 
(ii) mean gap length v. range 
from a proton ¢alibration : My-= 950+ 60m, ; 


from a 7-meson calibration : MR-= 1100+ 100 m,. 
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(6) The Secondary Disintegration 
On being brought to rest in the emulsion the K-meson produced an 
insignificant star of three branches, each 2-3 in length. Thus, no 
identification of these, the emitted particles, is possible. The visible 
energy of the star is of the order of 1 Mey. 


(c) The Parent Star of the K -meson 


The star from which the K -meson emerges is probably of type 4++1p 
or possible 4+2n. Three black tracks of ranges 40, 200u, and 12p, 
and one grey track, that of the K-particle, are arrested in the emulsion ; 
the two which are at plateau ionization are steep and may be followed 
- out of the emulsion stack. 


$5. Discussion 


The total visible kinetic energy in the parent star of event K (Br 1) 
is surprisingly small, being ~ 200 Mev. It appeared worthwhile therefore 
to examine whether this low energy in the primary star was peculiar 
to this event or was a general feature characteristic of the parent 
stars of other K -mesons. An analysis of similar events hitherto 
reported shows what may be an interesting, if not easily explicable, 
property which can be semmarized as follows : fig. 2 shows the distribution 
in N, and in n, for the parent stars of the various categories of 
heavy unstable particle events and from these histograms it would appear 
that the K -mesons are, in general, ejected from stars which have small 
values of NV, in contrast to those of the other heavy unstable particles. 
In fact, no well-established K -meson has yet been observed to originate 
in a star with N ,>12—a figure less than the median value of N, for the 
parent stars of the decaying heavy mesons. - 

The n, distributions, on the other hand, show no very significant 
differences for the various types of particles, though even in this case the 
median value of n, for stars from which K -mesons emerge, is the lowest 
of all the groups. From a comparison of the histograms of fig. 2 it would 
appear that the NV, and n, distributions for the 7-meson and K™ -meson 
origins are not significantly different. 

In event K (Br 1) the parent star is produced by a neutral particle. 
Tf one demands that the K -meson is produced in a reaction of the type 
a+N-->Y-+K, then from charge conservation and the known short 
lifetime of the neutral 7-meson, it follows that here one is dealing with a 
second generation process. In this connection it is of interest to consider 
the data summarized in the table, from which it appears that the ratio of 
charged primaries to neutral primaries for the parent stars of K -mesons 
differs from that for the parent stars of the decaying K-mesons. This 
may perhaps be explicable in terms of the preferential production of 
negative K-mesons by the charge exchange scattering of neutral heavy 
mesons. In view of the small statistics involved, however, it must be 
borne in mind that these apparent systematic differences may be present 
merely as a result of statistical fluctuations. 
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Figures 2 (a) and (b) show the distributions in N,, and n, of the parent stars for 
various types of heavy unstable particles. 

The cross hatched squares refer to events observed by Dahanayake 
et al. (1954) in which the particles were found by following out grey 
tracks from selected stars with n, >2, and not by scanning for the meson 
decays. Besides the bias against n,—0 or 1, there may also be a possible 
bias against stars with small NV, owing to scanning loss. The median 
N,, and n, values are indicated by vertical arrows. 

The data for hyperons and heavy unstable fragments are included 
here in view of their known association in the production of at least 
some types of K-mesons. 
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Perarie Kz K — single Unstable 

i charged secondary i fragments 
neutral 6 22 10 Lie 3 
charged 5 52 11 10 10 


In event K (Db 1) the small amount of visible energy associated with 
the absorption of the meson is striking. From this observation it would 
appear fairly probable, that in some cases, the capture of the K -mesons 
in emulsions will result in nuclear disintegrations with no observable 
secondary products. Examples of this kind have been reported previously 
by von Friesen (1954). It is possible to interpret such events, in which 
almost the entire rest mass energy of the K-meson is carried away by 
neutral particles, in terms of either of the following processes : 

(1) K +p— Y°+79+Q 

Q) K +p+n+B+@ 
where (=rest mass energy of the K-meson, B°=a neutral particle and 
the other symbols have their usual significance. Process (1) is the 
inverse of the ‘associated production mechanism’ established by the 
Brookhaven group. The production of neutral hyperons as a result of 
the nuclear capture of negative K-mesons has been well illustrated by the 
observations of Barker (1954) and De Staebler (1954). The kinematics of 
process (2) have already been considered by Friedlander e¢ al. (1954). 
If B°=7°, the energy and momentum considerations for processes (1) 
and (2) are very similar. If B°=v then the process (2) is analogous to the 
capture of «.~-mesons by all nuclei. There is as yet no direct evidence 
in favour of process (2). 
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ABSTRACT 


The @°-meson may be a scalar or a vector particle. The two possible 
interactions which account for its decay into two 7-mesons are considered. 
The cross sections (a) for the anomalous decay of #° into two 7-mesons 
and a y-ray, and (b) for the production reaction 7~+p — A°®+6°, are 
calculated, in the hope of finding a criterion to distinguish between the 
two spin values. 


§ 1. PossrpLe INTERACTIONS 


THE existence of a neutral particle 6° which decays into a charged 7-meson 
and one other light particle is now firmly established (Butler 1954). 
On the assumption that the second decay product is also a z-meson, the 
@ value is 214 mev. 

If the spin of the 6°-meson is assumed to be less than or equal to one, 
considerations of angular momentum and parity show that it must be 
a scalar or a vector particle. The interaction Lagrangian density, 
giving rise to just the observed decay is thus either 


La) f0( x)" (2) Grae ee ee oe, CT a) 
et LX) =ig0,2)| 8a x dle) — ED giey |, (1.2) 
yah [9,(a), ,(a”)] = (5, oh on x i=) Ategy Re ue ee 


Here @ and ¢ are 6° and 7 meson fields and M is the 6° mass. 

The vector Lagrangian (1.2) is, however ,incomplete. In fact if the 
electromagnetic interaction of charged mesons is also considered, gauge 
invariance of the theory requires an additional term involving the 6, field, 

L (x)= — 2iegd,(x) A ,(x)h* (x) (2). eer 3 (ees) 

We consider now the scalar and the vector interactions for renormali- 
zability. That (1.1) is renormalizable is well-known (Ward 1950). 
It is remarkable, however, that though the interaction of a neutral vector 
meson with a Fermi field is renormalizable (Matthews 1949), its direct 
interaction with a charged boson field is not renormalizable. In this sense 


interactions (1.2) and (1.4) may be considered as local approximations to 
the correct interaction. 
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§ 2. CRITERIA FOR DISTINGUISHING SPIN 
(a) Anomalous Decays 


At present a large number (about 20°) of anomalous decays are 
reported. Just as in the case of 7° decay (Dalitz 1951), one would expect 
that a small number (about 1°) of apparently anomalous cases are due 
to the alternative decay mode 


0° > at+a-+y. 

On account of the presence of L,? term in L,, in which this process takes 
place directly, it was thought that this effect might be more pronounced 
for the vector than for the scalar theory, and might provide a possible 
test for the nature of the 0°. An outline of the calculation is presented 
below. 

The first order Feynman diagram for 6°—»a+-+-7~ decay is shown in 
fig... 

Fig. | 


[=] 


The transition probability is w, where 


Here a= M?2—4y?, b=2M, and p is 7-meson mass. 
Figures 2(a) and (b) give the lowest order diagrams for scalar 27, y 
decay, while fig. 2 (c) gives the additional contribution from L,’. 


Fig. 2 


Bet ok 
(a) () (c) | 
The calculation of transition probabilities involves ae se over final 
states of three particles. If k,, ky, k are the momenta of the two ESE 
and the photon, respectively, the final integrations Pe peenon SORTEIOS 
w,andw,. Setting k,-+k,=2k’, w, for example equals {(k’+-k/2)?+-p7]*". 
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In the subsequent calculation this has been approximated to by 


[h’2 + h2/4tp2]V2p1 +k. k’/(e?+h/4+y?)]. . . . (2.1) 
Here k, k’ are magnitudes of k and k’ respectively. 
The results are f Anas 
UB y= (Z) (= -) Soh ai) f hy dls eos tes seat 


ele bs ae 

i) = FM eb) max” f° 
It may be noted that (2.2) has an infra-red divergence, so that a lower 
cut-off for observed phofon energies has to be introduced. Taking | 


kmin== 1 Mev. 2 2 2 1/2 
zs Why = UY kde. Yaad bab! te ORE 
Ht 4a} \4a] \3Mar 15 
and Wage ees bier : 
Tt (1.58) (=) Ph” oes 


For the vector case 
max 


oe 3 Z 1 /¢ "e2 Alyse ; ; 
hy lug) a ne 


w,” is the contribution from figs. 2 (a) and (6), w,” from fig. 2 (c) and w.,” 
comes from their interference. Here 
, lL fk(a—bk)3/*(M—k) M?(a—bk)>” 
A=5 ( (b—k)? ) k(M—k)(b—k)?’ 
2 (a—bk)3/? 
hh=9 (6—k)(M—k)’ 
1 k(a—bk)}? 
Is= § (ky 


w,” is infra-red divergent, and in fact gives the main contribution for 


ho even. > 
’ g° 6.8 
Wer y = (f) ele *) a. Xcel 2) gg een COREE 


Wen y m0 08 e a 


It thus appears that the ratio of 27, y decays (with y energies greater 
than 1 Mey) as compared with the total number of 27 decays is slightly 
greater for the scalar than for the vector case. The difference, however, 
is not significant. Also, at best, this ratio is of the order of 1%, so that 
it is not possible to explain most anomalous cases as y decays. 

We obtain a different picture, however, if the emission of only high 
frequency photons (say in the range 150—165 Mev) is considered. In 


this case 
dl hex OWES 2 M 
~\ ar) ar] \8Ma) \Sa50)> 28) 


min 


and 
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This means that events of this type are 10 times as frequent for vector 
as for scalar interaction, though the probability of 27, y decay occurring 
at all in this frequency range is 10-* times the probability of 27 decay. 


§ 3. THE PrRopuction REACTION 
A more promising criterion may be the study of the angular distribution 


of the reaction n+p > A960, 


This reaction is now considered to provide the production mechanism 
for 0° particles. Theoretically, however, one could not set up even a direct 
scalar renormalizable interaction for this reaction and it must be assumed 
to proceed via a (virtual) intermediate step +p —> n—A+6. The inter- 
action Hamiltonian for A> 6-+n would be Hi, vs, ,0°+c.c. for scalar 
6° and x, y,i%b,8,+¢.c. for vector 0° case. It may be noted that this 
three (rather than two) field vector interaction is itself unrenormalizable, 
and is postulated only by analogy. Let +k, —k be the momenta of the 
proton and z~ in the centre of mass system, k’, —k’, those of A® and 6°, 
and let angle ¢ be defined by k . k’=kk’ cos ¢. 

Then o(¢) = A+B cos ¢, for both scalar and vector 6°. 


Here A,=2 €ye3(k?+ €1€.)—m,?(€,€3+m Ms), 
B,=kk' (2k? +m 2+2€,€5), 
A ,5=3MyMo2Mg— € 1M 92€ 3+ 2(K?+ €1€y)€g€g— 2(k’?+ Egeq) eye 4M g?/M4" 
+4(k?+ €,€5)(k’?+ eseg)eney/m4”, 
By=—kk' [—2(h? + €1€9)—1m 9? + 4(k2 + €y€9)(k? + ge 4)/m4" 
2m.” 


ate ae (2+ ege,)]. 


4 


Numbers 1, 2, 3, 4 refer to p, 7, 4 and @ particles and ¢ stands for the energy. 
For Brookhaven 7 beam (beam energy 1 500 Mev) 


a,(p) cc 2-5+ cos 4, 
o,(¢) oc 15—9 cos ¢. 


It must be pointed out that the above cross sections have been calculated 
using first order perturbation theory. A more accurate calculation would 
replace y; in pion—nucleon interaction by the vertex operator A; (Dyson 
1949). In view of the unrenormalizability of the direct coupling for vector 
mesons, one would get into theoretical complications with indirect 
couplings, if the mesons did turn out to have spin one. 
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ABSTRACT 


In the course of structural analyses of aluminium-transition metal 
compounds, attempts have been made to find evidence for electron 
transfer effects proposed by Raynor. The results obtained are ambiguous 
if the intensities have not been measured on an absolute scale. Previous 
interpretations based on (f,—/F.) syntheses must be revised because 
deficiencies in the f-curves used in the calculations have been revealed. 
It is also shown that the effects of overlap, series-termination and limited 
counting area can alter appreciably the electron counts on F, syntheses. 
The ultimate possibility of detecting transfer is considered. 


§ 1. INTRODUCTION 


THE electron density distribution in a crystal can be derived by Fourier 
synthesis of x-ray diffraction data. This distribution can be regarded as 
the summation of a set of individual atoms, and analysis of data in terms of 
atoms involves the use of f-curves, which are theoretical calculations of the 
scattering due to isolated atoms. As the techniques of measurement and 
interpretation of diffraction data are improved, details of the continuous 
distribution are receiving more attention. This paper attempts to examine 
a particular set of results on electron distribution concerned with a pro- 
posed electron transfer in aluminium-rich intermetallic compounds, but 
the principles involved in the examination are of more general interest. 
The background to this work has been reviewed elsewhere (see Raynor 
1949, Taylor 1954). For the present purpose, we need only state briefly 
the problem to which the x-ray techniques have been applied. It has 
been proposed that in aluminium-rich compounds of the transition 
metals from chromium to nickel, electrons are ‘absorbed’ into the 3d 
shells of transition atoms to an extent which varies from 0-61 electrons for 
nickel to 4:66 for chromium. Attempts have been made to test this 
hypothesis by the use of diffraction data. Different and distinct methods 
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of treating the data have been applied by several workers (Douglas 1950, 
Robinson 1952, 1953, Nicol 1953) and these methods will be classified and 
distinguished in a review of the results in § 2. Specific details of the inter- 
pretation will be examined in subsequent sections. 


§ 2. Review or THE Metuops UsEep 


(a) The first technique uses only observed intensities in the preparation 
of electron density syntheses (#', syntheses). In the maps of electron 
density, the number of electrons associated with each of the atoms is 
obtained from the summation of the electron density within an area 
around the atom centre. Experimental measurements of the intensities 
used have not been made on an absolute scale, and this means that nothing 
is known about the existence or magnitude of any continuous distribution 
of electron density. All that can be measured from the synthesis is the 
relative weights of atomic peaks above the background. If it is assumed 
that this background is zero, then the total count over the peaks in a 
projection can be equated to the total number of electrons given by 
chemical analysis. Such a method applied to the results for Co,Al, 
(Douglas 1950) showed that the cobalt atom has 29 electrons instead of its 
normal complement of 27. 

(b) In the work of Robinson on f (AIMnSi) (1952) and Nicol on MnAl, 
(1953), a different method of scaling Fy counts was applied. The F, 
values were scaled in each of a series of regions of sin 6/A to have the same 
mean value as a set of structure factors (/’.) calculated with a Viervoll and 
Ogrim (1949) scattering factor for manganese and a James and Brindley 
scattering factor for aluminium. When the Fy values had in this way been 
placed on an ‘ absolute ’ scale, the F (000) term could be inserted from the 
total number of electrons per cell and electron counts were read directly 
from the synthesis. The results indicated an excess of electrons on the 
manganese atoms and also gave an appreciable (about 30 per cell out of a 
total of 400) count of electrons in the background between atoms. The 
interpretation for results on the 8 (AIMnSi) compound was uncertain 
because the measured chemical and (ideal) structural compositions differ 
appreciably, and there is in general no guarantee for these structures that 
all of the atom sites are always fully occupied. 

(c) Overlap and series termination may influence counts on F, syntheses 
and (f,—F.) syntheses have been used to eliminate their effects. These 
syntheses show how the electron distribution in the structure deviates from 
the ideal assumed in the atomic f-curves used to obtain the F.’s. The 
F, values have been scaled to F. values as in (b), so that the total count 
over the difference synthesis must be zero. Because the absorbed elec- 
trons are in 3d shells, their scattering is only appreciable at low angles, and 
an (f)—F,) synthesis using only low-angle data should reveal their 
presence. Syntheses of this type have shown for manganese in B (AIMnSi) 
an excess of one electron (Robinson 1952) and for manganese in MnA\l,, 
an excess of 1-5-2-0 electrons (Nicol 1953). In another method used by 
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Robinson (1952), certain low-angle reflections to which the transition 
atoms make a critical contribution were selected and their observed 
intensities compared with calculations for different assumptions about 

electron transfer. This is a variant of the (f)—f,) method which employs 
only particular reflections instead of the synthesis of all of them. 


§ 3. Tae Atomic Scarrerina Factors AND (f,—F,.) SYNTHESES 


Methods of type 2 (b) and 2 (c) depend for their reliability on the f-curves 
employed. Cochran (1953) has suggested a method of checking these 
using artificial syntheses. This check determines the number of electrons 
contained within a certain radius in the atom to which the f-curve corres- 
ponds. This has been applied to the scattering factors for aluminium, 
silicon and manganese: the procedure of Cochran was followed in detail 
and the. results shown in table 1 correspond to an atom radius of about 3A. 
They indicate that these f-curves correspond to atoms deficient in electrons, 
and in particular, that the Viervoll and Ogrim (interpolated) factors for 
manganese correspond to a relatively greater deficiency than the factors 
for aluminium. 


Table 1. Check by Artificial Syntheses on Numbers of Electrons 
Corresponding to Different Atomic f-Curves 


No. of electrons 
Atom zZ within 3:04 of 
atom centre 


Aluminium 13 James and Brindley 12-2 
Aluminium 3+ 10 5 Be 9-6 
Silicon 14 56 m 13-3 
Manganese 25 Viervoll and Ogrim 22-6 
3 (3d) electrons 3 ° Fs 3-0 
Mn + 3(3d) 28 . ie 25-6 


In the work of Nicol and of Robinson, the manganese Viervoll and Ogrim 
(1949) and the aluminium James and Brindley f’s were used. Viervoll and 
Ogrim’s results allow an estimate of the effect of extra 3d electrons to be made, 
and this was used by Robinson. 


To show how these results may affect previous interpretations, the 
results for one synthesis for MnAl, are discussed. Table 2 shows the peak 
weights in the asymmetric unit of the (hk0) projection, for an Fy) synthesis 
on an absolute scale for normal atoms (column (a)) and for an #, synthesis 
using the f-curves as in § 2 (b) (column (})). The difference is shown in 
(column (c)) but this would not be observed on an (/,—F,) synthesis 
because the coefficients used in (a) and (5) are first scaled te each other. 
This can be approximately adjusted by multiplying all figures in (b) by a 
constant factor to give the same total weight as for (a); the results are 
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given in column (d). The (#')—/’,) results of column (é) are then obtained. 
Comparison with Nicol’s values shows that his (/y—F) effects are larger 
than those given by the f-curve anomalies, but this calculation leads to the 
conclusion that the excess of electrons per manganese atom for this syn- 
thesis is now only 0:5. It will be shown later, however, that the real 
excess may be greater than is indicated by the above argument, for it may 
have been reduced in the operation of scaling Fy to F’, (cf. §§ 4 (c) and (d)). 


Table 2. Weights of Atomic Peaks in (hk0) Projection Syntheses of 


MnA\l, 
(a) () (c) (d) (e) (f) 
F, ie (F,—F.)| F. after | scaled Nicol’s 
synthesis | synthesis sealing | (Fy—F.) | (Fo—F¢)' 
to Fy result 
Mn 
(2 atoms) 50 45-2 4:8 48-6 1-4 2-4 
Al, 0 e 
(2 atoms) 26 24-4 1-6 26-2 0-2 0-7 
Al, 
(2 atoms) 26 24-4 1-6 26-2 0-2 13 
Al; 
(1 atom) 13 12-2 0-8 13-1 0-1 0-2 


The numbers of the atoms correspond to those used in fig. 3 (a) of Nicol (1953). - 


A similar correction would reduce the effects found by Robinson for 
B (AIMnSi), both on (#)—/.) syntheses and by comparison of low-angle 
reflections. It must be emphasized that this discussion has proceeded 
on the assumption that the background of electron density is zero. 
Because Fy is scaled to Ff, the difference synthesis can only show abnor- 
mality of peak-height ratios, and the interpretation is still ambiguous 
in the absence of absolute scaling. 


§ 4. THe F, SYNTHESIS 

(2) The ambiguity involved in Fy) counts assuming zero background 
(§ 2 (a)) is easily understood. This assumption for Co,Al, led to the result 
that cobalt has 29 electrons. _ It could alternatively be assumed that this 
atom has its normal complement of 27 electrons, that the number’ of 
electrons in aluminium is reduced to give the observed ratio of peak 
weights for cobalt to aluminium, and that the excess of electrons is distri- 
buted in a background density. For Co,Al,, this latter interpretation 
would indicate a background of 0-2 electrons per A%, or, in other terms 
would show that in a volume occupied by nine aluminium and two cobalt 
atoms there are thirteen electrons which are more than 0-74 from the 
atom centres, an alternative explanation which is physically quite 
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reasonable. The distance of 0-7 is the radius of the circle over which the 
count for an atom is made on the projection, and this radius is limited by 
_ the overlap of the atomic peaks in the projection. 

- (6) All counts made in this work have been restricted in the above way 
and even a radius of 0-74 usually means that circles from adjacent atoms 
are overlapping. Consideration of radial charge distribution for Hartree 
atoms (Hartree 1933) will show that in Cut, for example, there are about 
three electrons lying more than 0-74 from the atom centre. About half of 
these would project outside the counting circle of that radius. A similar 
estimate for aluminium cannot be made directly : the difference between 
U(r) curves for Al?+ and Al?+ (James, Brindley and Wood 1929) shows that 
the extra electron is almost entirely at more than 0-74 from the centre. 
Al?* has altogether 1-5 electrons outside this radius and it seems reasonable 
to assume that normal Al would have a greater proportional ‘ spread ’ 
than the transition atom. In that case electron counts on a restricted 
area would exaggerate the transition metal to aluminium ratio. Thermal 
vibration or the use of an artificial temperature factor would increase this 
effect. It is not possible to apply any quantitative estimate because the 
atoms in the structure are not isolated Hartree atoms. 

It is important to interpret correctly the appearance of an Ff’) synthesis ; 
this is a map of electron density (p) and a very small change in p at large 
distances from an atom centre may correspond to a total change of several 
electrons in a shell of large radius. The effect will, however, be swamped 
by the diffraction effects due to the large central peak of p ; because this. 
peak falls to a first diffraction zero at about 0-34 from the atom centre (a 
typical value for the work considered here), this does not mean that a 
counting area with even twice this radius will include all significant. 
contributions to p. 

(c) To illustrate the effect of these and other factors, some results 
obtained in the course of calculations on a projection of the Co,Al, 
’ strueture will be considered. This structure was solved by Bradley and 
Cheng (1938) using powder methods. Single-crystal (eye-estimated) 
intensity data collected by Dr. Newkirk in this laboratory are being used 
to determine more accurate parameters for the structure. At present, the 
R-factor for the (hki0) zone is 0:095. Counts have been made on Fy and 
several F, syntheses for a projection down the c-axis of the hexagonal 
unit cell: these calculations are made for a cell based on the orthohexagonal 
axes. The results are reproduced in table 3, and the positions of the peaks 
and counting areas are shown in the figure. All of the counts were made 
on syntheses calculated at 6° intervals ; a comparison for one synthesis 
calculated at 3° intervals showed that peak weight ratios could change by 
as much as 5°% on using the finer grid. The relationship of ait atom 
centre to the grid points of calculated p may affect the count, especially if 
the grid is too coarse. This factor may explain the difference between 
the cobalt peaks 5 and 6 which are related by hexad symmetry, but appear 
in different positions relative to the ortho-hexagonal grid. 

M2 
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The F, synthesis, scaled by assuming a zero background (method of 
§ 2(a)) would seem to show an excess of electrons on cobalt. An F, 
synthesis (F,,) drawn for Fermi-Thomas cobalt and James and Brindley 


Table 3 
Zero Zero Scaled to * Absolute 
background background F, zero scale ’ 
background 
Fy Ps di bigs 
Al 1 9:7 10-2 9-3 12:5 
2 12°3 13-3 12-0 12-8 
3 12-6 13-0 12-8 11-8 
Co 4 24-2 27:2 23-7 25-9F 
5) 29-4 28-2 29-2 25:3 
6 29-3 27-7 29-6 25-1 
*Co/Al mean ratio 2-36 2-20 2-37 2-09 


Normal Co — 27 
Ral Gy See ee 
(1) Normal Al 13 v3 
( Fermi—Thomas Co to 20cL ee 
James and Brindley Al 12-2 
+ Since the counting area for 4 is effectively a complete circle whilst the circies 
of 5 and 6 are cut off by chords, a change in the background adds more to peak 
4 than to the others. 


t osAU 


c-axis projection for CozAl;, showing positions of counting areas for peaks. 
Numbers correspond to table 3. OAB is the asymmetric unit. 
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aluminium atoms, without a temperature factor, shows individual 
variations between peaks of the same kind of atom comparable to those of 

the Fy synthesis and gives a mean Co: Al ratio a little higher than the 
f-carves would predict. If the coefficients of this synthesis are scaled to 
those of the Ff) synthesis, then the synthesis F,,is obtained. This scaling 
will include the effect of temperature factor, and low angle variations due 
to x;—%, resolution and extinction (neither of the latter can be effectively 
corrected). No F' (000) term is inserted and the peak weights are obtained 
by assuming a zero background as in F,. The results reproduce the 
abnormal Co : Al ratio, and the application of such scaling factors can thus 
make a critical difference to this type of result. The individual variations 
between peaks however, do show that the combined effects of limited 
counting area, of overlap and of series termination cannot be generally 
predicted, and that each synthesis will have its own features, so that 
results on J’) syntheses cannot be assessed without careful and individual 
investigation. 

(d) F., is the same synthesis as /., with an appropriate / (000) term 
(i.e. corresponding to the number of electrons per cell) inserted. Electron 
density values are then read directly off the synthesis and there is no 
assumption about zero background. This is the method of § 2 (b). There 
is in fact an appreciable background of electron density between the atoms 
in this synthesis and the ratio of peak weights is also changed as they are 
reckoned from a different background. The existence of this background 
on an F’, synthesis is explained because the / (000) term corresponds to 
normal atoms whilst the /-curves correspond to deficient atoms, or alter- 
natively, because the f-curves correspond to atoms whose electron density 
spreads appreciably beyond the counting circle. 

This method of scaling electron counts was used by Robinson (1952) 
and Nicol (1953) for /) syntheses, after the F’) coefficients had first been 
scaled to F, coefficients, so imposing on the fF’) values a mean angular 
variation corresponding to that of the #, values. In both cases an appreci- 
able number of electrons was found in the background, but it can be seen 
by comparison with the previous paragraph that this background is an 
artefact arising from the method of scaling. The introduction of this 
background would distort the ratio of peak heights. The scaling to PF, 
coefficients may give an approximately correct angular variation, but it 
cannot indicate the correct F (000) term—i.e. it can never show up a 
background density on an F')— F, synthesis. 


§ 5. Tae Derection oF ELECTRON TRANSFER 


Bijvoet and Lonsdale (1953) have questioned the possibility of finding 
electrons in outer atomic shells. They suggest that shells of large radius 
only affect spectra at low angles where no diffraction spectra exist. Ina 
discussion of the results of previous workers for LiH they criticized the 
extrapolation of measured /-curves to zero angle, and showed that the 


162 P. J. Black on the Measurement of 


difference between the scattering powers of the normal and ionized atoms 
at the lowest observable angle was negligible. Their conclusion was that 
no evidence for or against ionization in LiH could be obtained by this 
means. 

The extrapolation of f-curves to zero angle is not justified and the only 
possibility of detecting any ionic effect lies in the direct measurement of 
low angle spectra. For the LiH case the ratio Fy,,/Foo (where 
Fiy=fu—fa and Fopo=futsx) is 0-65 for normal Li and H, and has 
been measured independently by three workers as 0-50, 0-50 and 0:55. 
There is a significant difference, which is in the right direction for ionization, 
but it is not possible to say more unless one has correct f-curves for atoms 
inthis structure. In considering the difference in scattering effects between 
H and H-, Lonsdale and Bijvoet use an f-curve for H~ which is valid for a 
free ion but a calculation of its radial electron distribution shows that the 
second electron is distributed at 0-8 to 2-04 from the atom, whilst the 
interatomic Li-H distance is only 2-14. Waller and Lundquist (1953) 
have calculated the scattering in this crystal by correcting for the effect of 
overlapping between neighbouring atoms, and assuming complete 
ionization, obtained good agreement with observation for their calculation 
of the F1,,/Fy99 ratio. They do not indicate the extent to which this 
agreement depends on the degree of ionization assumed. 

One can assert that a complete set of Fourier terms will give a complete 
representation of the electron density in a structure—there is no question 
of low-angle * cut-off’ in the data—but it may not be possible to assign 
the (small) electron density between atoms to particular shells or particular 
atoms. If the Fy) term is not measured, the critical factor is the order of 
magnitude of the effect of any electron-transfer on the lowest angle spectra. 
Consider a structure with two different types of atom, separated by a 
distance d’. If an electron is transferred it should lie in or within a shell 
of radius d’/2, and the f-curve of this shell is given by 

f= | Uy)? a 
J 0 


ar 


where p.=47 sin 6/A. For a thin shell, this becomes 


fe sin (27 sin @/d) d’ 
~ (27 sin 6/A) a’ 
If d is the longest possible inter-planar spacing then the lowest angle 
spectrum occurs at sin @=A/2d, and the contribution of this electron is 
ie sin (7d'/d) 
va ah Nee 
For an NaCl type structure, f,=0-15 (the least favourable case). For 


the complex aluminium-transition metal structures, there are spectra 
trom planes of spacing 64 to 124, whilst d’=2-5A. Then f,=0-76 to 0-93, 
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and the effects should be detectable especially if the transfer involves two 
or three electrons. The complexity of the structure is, in this respect, 
an advantage. 


§ 6. CONCLUSION 


It is concluded that the existence of electron transfer effects in inter- 
metallic compounds has not been reliably established. It appears (see 
e.g. § 3) that the effect must be smaller than the transfer required by 
Raynor’s theory. If one can obtain accurate atomic positions for a 
fairly complex structure and place the intensities by experiment on an 
absolute scale, then phenomena involving two or three electrons should be 
detected. The best procedure would be to derive empirical f-curves for the 
two types of atom from the data, and then to calculate the numbers of 
electrons to which these curves correspond. There is no advantage in 
using calculated f-curves unless these correspond to some significant 
theoretical model. It may be useful for example to use an f-curve for 
the core (up to and including the 3p shell) of a transition atom which can 
be calculated theoretically and which one would not expect to be influenced 
by the structural environment, and use this in (/)>—F’,) syntheses to show 
only the 3d electrons. 


Lam grateful to Professor N. F. Mott and Dr. W. H. Taylor for provision 
of facilities and for their interest in the progress of this work. I wish to 
thank the Council of the Royal Society for financial support. I am 
indebted to Dr. J. B. Newkirk for permission to use his measurements on 
Co,Al,, and to Dr. Cochran and other colleagues for helpful discussion and 
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ABSTRACT 


We have attempted to fit experimental results for the scattering of 
electrons by heavy nuclei at 126, 154 and 183 Mev with curves obtained 
from various assumed nuclear models by way of a phase shift analysis 
carried out for the case of mercury. At these energies only nuclear 
scattering is important and the nuclear model completely determines 
the scattering. Charge distributions which reproduce the experimental 
curves either are homogeneous or have a slight dip near the centre of the 
nucleus and have their main decrease on the surface of the nucleus over a 
transition range of between 1:7 and 3:2x10-% cm. The root mean 
square radius of the distribution is determined to be approximately 
R=1-2x 10-8 x A’ cm, to a relative accuracy of about 5°%%. 


§ 1. Lyrropuction 


Ws shall report on phase-shift calculations of coherent Coulomb scattering 
of electrons by mercury. We consider only scattering from the nuclear 
charge distribution. Others of the effects usually considered have been 
discussed previously for this energy range (Schiff 1953) and seem to give 
only small corrections to the Coulomb scattering. At the energies 
considered, 125 and 180 Mev, the charge distribution of a mercury nucleus 
has a radius larger than the electron wavelength, and the scattering can be 
used to determine the size and shape of the charge distribution. 

We shall apply our analysis to the experimental data of Hofstadter, 
Hahn, Knudsen and McIntyre (1954), and shall show that the scattering 
is sensitive to the range of the charge distribution and to certain general 
features, especially the length of transition region in which the charge 
goes from its average value inside the nucleus to zero. 

A survey of nuclear models has been carried out by Yennie, Ravenhall 
and Wilson (1954), Ravenhall and Yennie (1954) and points to the general 
conclusions we draw here. We have extended the accuracy of previous 
calculations by Brenner, Brown and Elton (1954) which were inaccurate 


* Communicated by Professor R. E. Peierls, F.R.S. 
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at large angles as noted in that paper. The general agreement of present 
calculations is very satisfying considering the large number of phase 
shifts and high accuracy required. 


§ 2, PROCEDURE 
The procedure followed is similar to that outlined by Brenner, Brown 
and Elton (1954), which will be referred to as I, with some modifications 
which make the calculation more convenient. The notation of I will be 
employed. 
We are concerned with solving the pair of equations* 
ads k OM en yee 
ge tt R= 0, Seger m, SS = 0.0 omen (L.) 
where & is Dirac’s quantum number, i.e. k=j+-4 if j=l—4 and k=—j—} 
if j=/+ 5. Here j and / have their usual meanings. The solutions 
become asymptotically 


R~cos (r+Za log 2r—gkn-+-x,), ] (2) 
S~sin (r+ Za log 2r—fhm+y,). | 7 7 


The scattering cross section can be expressed in terms of y, (Mott and 
Massey 1949, p. 76). We first define 


ice) 


AC EO cep (ie) TP, (con 6) EPL (cos Oot el ns) 


vA PA) heal 
where we have reduced the expression of Mott and Massey to that for 
high energies ; i.e. we have set equal the two phase shifts corresponding to 
the same 7. Then the cross section is given by 


do mG 
ID lf Psec*s- Be Se th hy oc SCA) 


We split y, into 7,+6, where 7, is the Coulomb phase shift for a point 
charge and 5, is the ‘ anomalous phase shift’, which is different from 
zero for only the first few values of k. As explained in I, we express 
f(0) as fp(0)+[f(0)—fp(0)] where fp is the amplitude which would result 
from a point charge and f—fp involves only a few terms in a summation in 
k. We can easily determine 5, from a knowledge of R/S of (1) and of the 
regular and irregular point charge functions at some point beyond the 
nucleus. We have found a more convenient form than that given in I 
for the asymptotic series for the regular and irregular Coulomb functions 
which we denote by Rp, Sp and &;, S;. This is 
Sp=(X+U) sin o,+(¥+V) cos o;, | 
R,=(X—U) cos o,—(Y—V) sin o;, 
S;=(X+U) cos o,—(Y+V) sin o,, 
R,;=—(X—U) sin o,—(Y — JV) cos o,;, j 


(5) 


* We choose units such that /=c=—H=1, and neglect the ¢iectron mass in 
comparison with the electron energy. It has been shown in I that the errors 
made in this neglect are of order (m/E)?. 
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where o,=1+Za log 2r—hkr7+n, - - - . . ~ (5.1) 
and we can define X, Y, U, V by a power series expansion, 
Ae . 2s ot, ieee e 
X+i1Y=A,+71By+ oh (4 ittBy)+ WN Qkry (A,+27B,)+ sors 
i ; e e ; : 
k(U+1tV)= hp (41 +1B,)4 (2h) (A,+7Bs)4 Nk) (A,+71B3)+ ... 
(5.2) 
where 
A,=1,.5)=0; Aye By 0 ee 
Yap VS dal ee et ob, A, ,=(k?—s”) A,—2s Za B,, o> 5, eRe aOsaE 


Byi4=28 Za A+ (k?—s") B 
We have evaluated R,, Sp, R; and S, to five decimal places for k=1... 10 
and for r=12-5 and 14 for each k. We should be glad to furnish the 
numerical values of these functions upon request. To obtain five-place 
accuracy, we have had to evaluate only ~k terms in the series for the kth 
point charge functions, since the term (k?—s?) occurring in the recurrence 
formula becomes small for k~s. 

We have evaluated the point charge amplitudes to five decimal places 
by the method of Yennie et al. (1954) and consequently have improved 
greatly the accuracy of the calculations reported in I, where we had 
employed the amplitudes given by Feshbach (1952). We give a summary 
of our point charge amplitudes.* 


st 


Point Charge Amplitudes 


35° 40° 45° 
Real part of fp _ 3069 2-655 2-2675 
Imaginary part of fp —2-095 —]-2]1 —0-6488 

50° 55° 60° 65° 70° 13° 80” 
19261 1-6320 = 1-3811 1:1684 0-988] 0-8354 0-7058 
—0:2870 —0-0538  0-0951 0-1873 — 0-2413 0-2691 0-2789 
85° 90° 95° 100° 105° 110° 115° 


05956 0-5017 — 0-42155(5) 0-35302 0-29437(5) 0:24415 —«0-20114 
0-2766 0:2659 024981 023029 0-20882 018644. ~—s«0- 1 6394. 


120° 125° 130° 135° 140° 145° 150° 

0-16430 0-13278  0-10589 0-08303 0-06372(5) 0-04755 0-03417 

0:14189 0-12071 0-10075 0:08224 0-06537 0-05028(5) 0-03707 

These amplitudes agree with those of Feshbach to within one unit in 
the last place given by him at the angles at which they can be compared, 
except at the smallest angle, 45°. The method used by us does not give 
accurate results at small angles, but we expect the accuracy at large angles 
(90°) to be to the number of decimals quoted. 


* The G(@) tabulated by Feshbach (1953) is the same as our Sp(9). 
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For the calculations at 125 mev we have had to evaluate 7-8 anomalous 
phase shifts and for the calculations at 180 mev we have employed 10 
phase shifts. These phase shifts have been calculated to five decimal 
places.* Consequently, the theoretical curves are accurately defined as 
shown in figs. 1-7. High accuracy is important because 2 to 3 significant 
figures in the amplitudes are lost through cancellation at large angles. 


§ 3. RESULTS 


We show in figs. 1-7 results from nuclear models with more or less the 
correct range to reproduce the experimental results. The range of the 
distribution is defined in terms of the square’ (homogeneous) distribution 
which would reproduce the same low energy effects (Feshbach 1952) and 
this range, in turn, is used to define a nuclear radius parameter 7, in the 
following way : 

5. 47 (© 
red dR v4 yp 
Zl Orit: Or. 
Beary ee 10st ems 
(In this section c.g.s. units will be employed.) 
We have considered charge distributions of the following form : 
(2) Homogeneous with Gaussian type cut off 
1 


R2= 


(7) 


(b) Wine bottle 
Lr? /a? 
teat 1+exp [(?—22)/B?] 

The quantity 6? is related to the length of transition region over which the 
charge decreases from an average value to zero. In (a), for example, the 
charge goes from 0-88K to 0-12K as r? goes from «*—26? to 02 2B ; 
from this we can definea transition region of length («7+ 26?)"?— («?— 26?)¥? 
which is approximately equal to 26?/« for B<a«. We can relate a to Rk, by 
(6) and can relate K to the total charge Ze and « by requiring the volume 
integral of the charge density to be equal to Ze. We have chosen 
R.=1-2« Ax 10-13 em for all distributions calculated here. Thus, all 
distributions have the same second and fourth moments of the charge 
density. . 

In fig. 1 the results of various curves of type (@) are shown, for scattering 
at 125 Mey. The experimental points are for scattering at 126 Mev from 
Hofstadter et al. (1954). The experiments were carried out for gold and 
for lead, and gave identical results, within experimental error, for these two 
elements. We shall therefore use them without correction for mercury. 


i ions integrated out from the origin to 

* As a check the fitting of the functions integra | the: n te 

the asymptotic Coulomb functions has been performed at two different places. 

The phase shifts thus obtained differ from each other by only a few places in 
the last figure and should therefore be correct to the same accuracy. 


168 G. E. Brown and L. R. B. Elton on 


Absolute measurements have not been made and consequently the 
experimental points may be moved up or down by a common distance on a 
logarithmic plot. It is seen, however, from the large angle scattering that 
by suitable normalization of the experimental points, curves from charge 
distributions with transition parameter B?3«?/16 to «7/4 fit the experi- 
ments reasonably at this energy. (We shall see that somewhat shorter 
transition regions fit the 154 and 183 Mev data better.) 

Figure 2 indicates the effect of a wine bottle distribution. The 
transition region in this case is similar (but not exactly equal, if we use the 
above definition) to that of the curve in fig. 1 with 6?==«?/4. The effect of 
scooping charge from the centre and putting it farther out is to raise the 
scattering at intermediate angles, 60° to 90°. Although curves of the type 
(a) fit reasonably the experimental points within the limits of error of these. 
the general trend is easier to reproduce with the type (>) curve. 

We can, in fact, see what the effects of different modifications in the 
charge distribution are. (In all such changes we keep the 2nd and 4th 
moments constant.) Heaping the charge up at the origin, as in an expon- 
ential or Gaussian distribution, raises the scattering at small angles in the 
region 30—70°, as portrayed in fig. 2. Raising the charge at a distance 
from the origin, as in the wine bottle, raises the scattering at intermediate 
angles, 60 to 90°. Increasing the transition range decreases the scattering 
slightly at intermediate angles and markedly at large angles, >100°. 
A comparison with the results of Yennie et al. (1954) shows that these 
remarks are borne out by their work as well. 

It is necessary to calculate cross sections only at 125 Mey, as long as the 
points are fitted to sufficiently wide angles because theoretical curves 
extrapolate with energy in a known way, namely, the ratios of the 
differential cross sections to the point scattering cross section is the same for 
different energies H and XH’ at corresponding angles, or 


do(H, 0) /do(E,0)p — do(H’, 6’) /do(E’; 0’) p’ 
OMe = ao] aa <= ee aa 
where FL, EH’, and @” are related by 
(HR) sin 46=(H’R)™ gin 46’. Sh et ees een) 


This semi-empirical extrapolation rule was suggested by Dr. McIntyre of 
Stanford and applied by Yennie et al. (1954). To illustrate the utility of 
this rule, we exhibit in fig. 3 theoretical curves at 180 Mev and curves for 
the same energy obtained by extrapolating the 125 Mev curve. In fig. 4 
the experimental points at 154 and 183 Mev are shown and the curve 
obtained by extrapolating the curve passing through the centres of the 
126 Mev points. (The differences between theoretical curves for 125 and 
126 Mey and those between curves for 180 and 183 Mev are quite negligible. 
In the following figures 5—7 we have, however, performed the small extra- 
polation from 125 to 126 Mey and from 180 to 183 Mev to make a direct 
comparison with experiment.) The extrapolated experimental curve 
does not pass through all of the 183 Mev points, but the experimental 
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Fig: 1 


Differential cross section in em?/steradian 


30 60 90 120 150 


Scattering angle in degrees 


Differential cross sections for scattering at 125 Mev from nuclear models of 


the type 
1 


p= ET rexpl—) Fl 
where the following values of K and « are associated with the given f*: 
(yest a=w 7) 16, a=4-40, K=0-001595 
(ii) B? 07/8, a=431, K—=0-00171,; 
(iii) P?=02/4, oO: K=0-00216 ; 
(iv) B?=0?/2; a=3-27, K =0-00300. 


(Lengths are measured in units of /ic/H, energies in units of #.) All 
distributions have nuclear parameter ry=l-2 (in ¢.g.s. units). (Inset : 
charge distributions.) 
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Fig. 2 


Differential cross section in cm?/steradian 


30 60 90 120 150 


Scattering angle in degrees 


Differential cross sections for scattering at 125 Mev from: 


distribution with Gaussian cut off 
K : 
pe T exp [(?°—0*)/B 
with  p*=a?/4, a=3-91, K=0-00216, 
(ii) Wine bottle distribution 
K 1+-72/a? 


with fB?=«?/4, a=3'68, kK =0-00146(5). 


(iii) Exponential distribution 
p=K exp (-r/«) 
with a==0-99, K=0-0236. 


(i) Homogeneous 


(Lengths are measured in units of ic/H, energies in units of HZ.) All dis- 


tributions have nuclear parameter 71-2 (in c.g.s. units) 


distributions. ) 


. (Inset: charge 


Elastic Scattering of High Energy Electrons by M. ercury 171 


errors are large in the region of disagreement and by drawing the 126 Mev 
curve in a different way, but still within the range of experimental error, 
the discrepancies can be decreased. 

The fact that experimental curves do not extrapolate very well with the 
above rule may indicate the presence of small effects not considered here. 
Order of magnitude estimates of the effects usually considered (Schiff 
1953) seem to indicate that they are much smaller than the differences 


Fig. 3 


—Exact | 
Extrapolated 


Differential cross section in em?/steradian 


Scattering angle in degrees 


i i i ti ated at 180 mev for the models 
Theoretical differential cross sections calculated a : 
ee “sted as (i) and (ii) under fig. 2 and curves extrapolated from those. ot 


fig. 2 by the rule, eqns. (8) and (9). 


between actual and extrapolated experimental Evel If we assume oon 
our three specific calculations that all reasonable theoretioa! te 
extrapolate with the above rule, and that there are no evn ef pe 
other than those considered here, we must ascribe ane aries ae. 
above to experimental errors, and we would then oD ae aM ae 
by trying to fit all energies simultaneously . Even if iat eed ih a : 
smcorrect and the discrepancies arise from theoretica! inadequacies, we 
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have no reason to weight results at a particular energy more heavily than 
those at another. We illustrate, therefore, in figs. 5-7 comparison of 
theoretical and experimental results at all energies. 


Fig. 4 
! 


d 


Ou 


Differential cross section in arbitrary units 


30 60 90 120 


Scattering angle in degrees 


Experimental points at 84, 126, 154 and 183 mev and the curves obtained by 
extrapolating to 84, 154 and 183 mev the line drawn through the 126 
MeV points. 


We show in figs. 5 and 6 curves of type (a) with B?=02/8 and p2=a2/4 
respectively (transition regions of 1-72<10-em and 3-20 10-33 em 
respectively) compared with experimental points for all experimental 
energies, and in fig. 7 curves for the wine bottle distribution. The 
theoretical curves at 84 and 154 Mev in all figures and at 183 Mey in fig. 5 
have been obtained by the extrapolation rule. 
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We see that the general features of the nuclear charge distribution can be 
ascertained reasonably. The nuclear radius parameter is in the region of 
fo=1-2. For curves of type (a), 79 can be increased by about 0-03, whereas 
for type (b), 7) can be increased by about 0-06 or decreased by about 0-04. 


Fig. 5 


Differential cross section in cm?/steradian 


183 6 9 154 


30 60 90 120 150 
Scattering angle in degrees 


. < > 3 =A ieee Z > A 
Theoretical cross sections for the model given below fig. 1 with B?=«7/8, com 
pared with experiment at all energies. 


without destroying the fit with experiment. We would say oe bani 
results, assuming that the correct model is not radically piereny sags 
those we have investigated, that ry lies in the region 1-16 to 1:26. We sith 
further, that the actual transition region would seem to lie between the 
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transition regions of the curves, figs. 5 and 6, or between 1-72 x 10~ cm 
and 3-20 10-13 cm. Not much can be said about the charge distribution 
near the origin, since the scattering is relatively insensitive to this region, 
but a slight wine bottle shape fits the trend of experimental points more 
easily than a homogeneous distribution in this region. 


Fig. 6 


Differential cross section in cm2/steradian 


30 60 90 120 150 


Scattering angle in degrees 


mL. 5 . a Set . w é ; 
Theoretical cross sections for the model given below fig. 1 with 2 =a?/4, com- 
pared with experiment at all energies. 


= ; 

We have received results from Dr. Ravenhall and Dr. Yennie which 
point to similar conclusions in the case of type (a) curves, although the 
details of the charge distributions employed are different ; in particular, 
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an exponential cut off is used instead of a Gaussian cut off, and conse- 
quently, numerical values are somewhat different. 

We should like to thank the staff of the Cambridge Mathematical 
Laboratory and in particular the Director Professor M. V. Wilkes for the 


Fig. 7 


Differential cross section in cm?/steradian 


30 60 90 120 150 


Scattering angle in degrees 


: e eee ae 
Theoretical cross sections for the model listed as (ii) below fig. 2, compared 
with experiment at all energies. 


facilities with which this work was carried out. We should also like to 
thank Drs. Ravenhall and Yennie for communication of thew results 
plicati an resting correspondence, as well as Professor 
before publication and for interesting correspoi me ae 
Hofstadter, who sent us the experimental results before publication. 
er, sent 
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ABSTRACT 


The large-angle scattering of .-mesons in iron has been investigated by 
a method similar to that reported previously for lead. The measured 
large-angle distributions are compared with the theories of Moliére, 
Olbert and a modification of Olbert’s theory which is based on a single- 
scattering cross section inferred from recent electron-scattering measure- 
ments. At energies from several hundred mev to about 4 Bev the 
measured distribution, at large angles, is not consistent with present 
Coulomb scattering theories in which the effects of the finite nuclear 
dimensions are included. As in the lead experiment, the number of large 
angles observed is in rough agreement with Moliére’s theory for a ‘ point ’ 
nucleus. Some properties of the so-called anomalous scattering, which 
follow from the lead and iron experiments, are discussed. 


§ 1. INTRODUCTION 


Ix a recent paper (McDiarmid 1954, referred to as J), measurements 
were reported of the j.-meson scattering distribution in a 2 cm lead plate. 
It was found, in agreement with other workers, that the observed number of 
large angles for particles with energy in the region 0-6 to 4 Bev was an 
order of magnitude greater than expected from Olbert’s theory of Coulomb 
scattering for a ‘solid’ nucleus. In fact, the large-angle distributions 
fell very near to those expected from Moliére’s theory for a ‘ point’ 
nucleus. 

The purpose of the experiment reported here is to examine (by a 
method similar to that used in I) the large-angle ‘ tail’ of the scattering 
distribution in iron, and from this examination to determine, if possible, 
the dependence of the large-angle cross section on the scattering nucleus. 
Again the measurements are compared with the theories of Moliere 
(1948) and Olbert (1952), as well as a modification of Olbert’s theory 
which is based on a more realistic single-scattering cross section for a 
‘solid’ nucleus. The single scattering law assumed by Olbert is the 
‘point ’ nucleus cross section for angles less than A/F (A is the reduced 
wavelength and # the nuclear radius) and zero for angles greater than A/R. 
Although the multiple scattering distribution derived from this is reasonably 
rc i eee 

* Communicated by H. J. J. Braddick. 
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accurate for large thicknesses and heavy elements, it considerably 
underestimates the number of large angles in the distribution for small 
thicknesses and light elements. Hence a simple modification of this 
theory is given in the Appendix which can be applied to a more or less 
arbitrary single scattering law. 


§2, EXPERIMENTAL ARRANGEMENT 


The cloud chamber containing six 2-54 cm iron plates was operated 
underground at a depth of 26 m water equivalent. The expansions were 
triggered by the simple counter telescope shown in fig. 1. This counter 
arrangement, as compared to the one used in the lead experiment, consider- 
ably reduced the correction for bias against observing large scattering 
angles (see I). It also decreased, however, the number of useful pictures 
taken per hour by slightly more than a factor 3. 


Fig. 1 
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Experimental arrangement. 


The rest of the apparatus and the method of measuring the angles 
were the same as described in I. Again, a check on the accuracy of 
measurement showed that systematic errors in the angle measurements 


were less than -++0-2° and that the r.m.s. error in a single measurement 
of angle was 0:5°. 


a 
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§ 3. RESULTS 


_ Using the large-angle cut-off mentioned in I the r.m.s. angles of multiple 
scattering <@), were obtained for about 5000 particles which passed 
through the chamber and traversed the six iron plates. The frequency 
distribution of the r.m.s. angles was divided into the following three 
groups: group 1, 2-6° to 6°; group 2, 1-1° to 2-5°; group 3, <1°. For 
groups | and 2 the pf spectra were determined as described in I; these 
are shown in fig. 2. The large-angle scattering analysis then consists 
of comparing the number of large angles observed in each group with 
the number calculated by averaging the different scattering theories 
given in the Appendix over the pf spectra of the corresponding group. 


Fig. -2 
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pB spectra for groups | and 2. 


Some of the particles in group | lose an appreciable fraction of their 
energy while traversing the plates and therefore only angles occurring 
in plates 3 and 4 have been used in the large-angle scattering analysis. 
The fact that the particles have to traverse the remaining two plates 
imposes the lower energy limit shown by the dashed line in fig. 2. The 
particles in group 2 are of sufficiently high energy that ionization loss 
is not important, hence angles occurring in plates 2 to 5 have been used 
in the large-angle analysis. In figs. 3 and 4 the measured integral 
distributions are given for groups I and 2 respectively, along with the 
distributions expected from the theories of Moliere and Olbert ne 
ro=1-0, writing R=r)xX 10-1341/3 em). The present modification O 
Olbert’s theory lies between these two theoretical curves : For a 
angles and energies considered in both groups the difference bet ge 1e 
two theories is small and it is evident that the measurements agree 
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about equally well with either theory. This can be regarded as giving a 
rough check on the method of analysis used. 

The pf spectrum (fig. 5) of the particles in group 3 was obtained by 
subtracting the resultant spectrum of the particles in groups 1 and 2 
from the underground spectrum corrected for scattering in the apparatus 
(see I). The underground spectrum used was derived from the measured 
sea-level spectrum of Holmes, Owen, Rodgers and Wilson (private 
communication) ; although this spectrum is not very different from that 
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Measured and calculated integral scattering distributions for group 1. The 
data represent 226 traversals of a 2-54 cm iron plate. 


assumed in I it is presumably more accurate. The spectrum shown in 
fig. 5 should be reliable at values of pB above about 0-9 Bev/c, but below 
this there is some uncertainty due to possible errors in the correction 
applied for scattering in the apparatus and the fact that here the spectrum 
is derived from the difference of two nearly equal quantities. 

Again, for particles in this group angles occurring in plates 2 to 5 have 
been used in the large-angle analysis. In fig. 6 the measured and 
calculated distributions are given and it can be seen that the measurements 
are quite inconsistent with either Olbert’s theory (with rg9=1-0) or the 
present modification of this theory (with r>=1-1). As in the lead experi- 


ment the measurements are in reasonable agreement with Moliére’s 
theory for a ‘ point * nucleus. 
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Measured and calculated integral scattering distributions for group 2. The 
data represent 1568 traversals of a 2-54 cm iron plate. 


An uncertainty exists in the theoretical distributions shown in fig. 6 
due to the uncertainty, mentioned above, in the low energy end of the 
pB spectrum of fig. 5 and the fact that low energies contribute appreciably 
to these distributions. However, this uncertainty is not great enough to 
affect appreciably the results, e.g. if the number of particles below 
0:9 Bev/c is uncertain by a factor 2 (which is certainly an upper limit) 
then Moliére’s distribution is uncertain by less than + 20°, at all angles 
and the other distributions by less than + 80°, at all angles. 


Fig. 5 
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The contribution to the expected distributions due to incoherent (i.e. 
inelastic) scattering can be estimated from the theory of Gatto (1953). 
Although this is relatively more important in iron than in lead it is still 
too small by a factor of nearly 10 to account for the number of large 
angles observed. In fact, adding this contribution to the * tail ’ of the 
theoretical distribution (‘ present theory ’) in fig. 6 increases it by about 
30% in the region of 6=7°. In connection with incoherent scattering 


it can be mentioned that no secondary particles were observed to be. 


associated with any of the large-angle events in either lead or iron. If 
these events are due to incoherent scattering of any kind (Coulomb, 
short-range forces, etc.) then nuclear disintegration products might be 
expected in some cases. Whether or not these low energy secondary 


Fig. 6 


ifexe) 


fe) 


ialaaiil 


Moliére 
‘4, - 


ya 


| 
| 
| 
/Present 


7 theory 
R= 


Olbert 
BLO 


2 eT PL eae | 


T 


Number greater than @—~- 


fe) 


4 5 6 7 8 9 
8 (degrees) —* 


Measured and calculated integral scattering distributions for group 3. The 
data represent 18 420 traversals of a 2-54 cm iron plate. 

particles are observed depends on the experimental conditions and it is 
estimated from the number of large angles observed that, if disintegrations 
occur, the number observable in both the lead and iron experiments would 
be considerably less than one. Therefore the fact that none was observed 
is not evidence against any type of incoherent scattering. 

The maximum energy to which the present measurements refer is 


probably considerably less than might be supposed by looking at the pp 


spectra for groups 3 of this and the lead experiment. Although the mean 
energy of the particles in both groups is about 8 Bev, as much as half the 


—————— 
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contribution to the theoretical scattering distributions (Moliére’s) in 
figs. 6 and 7 of I, comes from particles with pp less than about 2 Bev/c 
and only 5 to 10% from particles with p8 greater than about 4 Bev /c. 
_ If we assume that the measured scattering has the same energy depen- 
dence as Moliére’s theory (the fact that agreement is obtained with this 
theory in all groups of both experiments is evidence that it does) then 
follows that the measurements refer mainly to particles with pp less 
than 4 Bev/c. At energies above this the only information given is an 
upper limit to the large-angle scattering cross section. 


§ 4. CONCLUSIONS AND Discussion 


The conclusions drawn from the lead and iron experiments are the 
following :— 

(1) At energies from several hundred Mev to 3 or 4 Bev the jy-meson 
scattering distribution at large angles in lead and iron is not consistent 
with present theories of Coulomb scattering (coherent or incoherent) by 
a ‘ solid ’ nucleus. 

(2) At these energies, in both lead and iron, there is an ‘ anomalous’ 
large-angle ‘ tail ’ to the distribution, the magnitude of which is in rough 
agreement with that predicted by Moliére’s theory for a ‘ point ’ nucleus. 
In order to account for this anomalous component by the present 
modification of Olbert’s theory the value of the nuclear radius would 
have to be decreased by a factor 3 (from the value corresponding to 
aeeick ). 

(3) At energies of a few hundred mev the present experinients are not 
capable of distinguishing between the ‘point’ and ‘solid’ nucleus 
theories ; at most they set an upper limit to any anomalous scattering 
component present at these energies. 

Although more measurements are required, particularly at low energies, 
it is still possible from the experiments performed to date to indicate 
roughly some of the properties of the anomalous scattering. Within the 
statistical limits of the experiments (in particular the present ones and 
that of Whittemore and Shutt 1952), all the data are in agreement with 
the theory of Coulomb scattering for a * point ’ nucleus. This would seem 
to imply the following: (1) the scattering cross section varies approxi- 
mately as 1/(p8)? or 1/B? (# energy) ; (2) the differential distribution varies 
as 1/6” with n=3 to 5; (3) the cross section per nucleus varies as Z? or 
A2 but very probably not as Z or A. . 

Ag mentioned in I it does not seem possible at the moment to interpret 
the measured distributions at large angles. Several authors have suggest- 
ed that the anomalous scattering indicates a non-Coulomb interaction 
between fast y-mesons and nucleons. The difficulty with this suggestion 
is that any new scattering potential postulated would have to be consistent 
with the experiments of Conversi ét al. (1947) and Fitch and Rainwater 
(1953) on the interaction of stopped p-mesons with nuclei. Thus, such a 
potential, as well as predicting the properties of the anoinalous scattering 
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given above, would presumably have to approach zero at low energies or 
else give rise to a very short-range force concentrated at the centre of the 
nucleus. Neither seems very likely. 

On the other hand an explanation in term of Coulomb forces alone 
must be consistent with the electron scattering data of Hofstadter et al. 
(1953, 1954). Although the momenta, and hence the reduced wavelengths, 
of the electrons and y-mesons differ greatly, the maximum momentum 
transfer observed in both experiments is of the same order (~ 250 Mev/c). 
In the Born approximation the scattering cross section depends only on 
the momentum transfer (the maximum value of this is of course limited 
by the momentum of the particle), and hence in this approximation no 
difference would be expected in the two experiments (at least in the 
coherent scattering). The results of a phase-shift analysis might be 
expected to depend also on the momentum of the particle (the higher the 
momentum the more phase shifts have to be included in the analysis). 
However, such a dependence would not be expected to be great enough 
to account for the difference between the electron and u-meson experiments. 

Possible the most significant difference between the electron and y-meson 
experiments is the fact that only coherent scattering is measured in the 
former, whereas both coherent and incoherent are measured in the latter. 
As mentioned above, it is not possible to distinguish between these two 
types of scattering in the present experiment and this appears to be the 
case in all ».-meson experiments to date. Thus there is no experimental 
evidence against an explanation of the anomalous scattering in terms of 
incoherent Coulomb effects. The existing theory of incoherent Coulomb 
scattering (Gatto 1953) predicts the observed energy and angle depen- 
dence but is inconsistent with the magnitude and possibly the Z depen- 
dence of the measured scattering. In this connection a more exact 
theoretical evaluation of the interaction used by George and Evans 
(1950), in their analysis of stars produced by p-mesons, might be useful. 
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AP PEN DIX 


The following is a simple modification of Olbert’s theory (1952) which 
gives the multiple scattering distribution corresponding to a more or less 


ee oe LU 


Large-Angle Scattering of -Mesons in Tron 185 


arbitrary single scattering cross section. For reasonable assumptions 
about the scattering cross section the method gives a distribution which is 
accurate to about 20%, at large angles. The particular single-scattering 
law used in the present application is one inferred from a Born approxi- 
mation calculation by Rose (1948), for a uniform nuclear model (U.N.M.). 
The Born calculation for a U.N.M. gives pronounced maxima and minima 
in the scattering cross section ; these are not found in either the exact 
phase shift calculations for 125 Mev electrons (Yennie et al. 1953) or in 
the measured scattering cross section of these particles (Hofstadter et al. 
1953, 1954). However, a smooth curve lying along the maxima of ‘the 
cross section given by Rose (for r7=1-1 in R=r, x 10-1848) does fit the 
électron data reasonably well and it is such a single scattering law that 
is considered here. 

In terms of the scattering variable x, defined by Olbert and which 
refers to projected angles, the probability of a deflection between « and 
x+dz in consequence of a single collision in a thickness ¢ of scattering 
material can be written approximately 

6-5 1° [l—exp (—4C0G?x*)} 

P(e, Rh) da @ Te0@ zi dat, ee ANZ L) 
where the effects of screening, which are not important for the present 
purposes, have been neglected and where @ is defined by Olbert and 
O=0:510-7(Nt/A) Z4R* for B~1 and R in cm, N=Avogadro’s 


Fig. 7 


number, t=thickness in g/cm?, A=atomic weight, and Z=atomic 
number. For 2cm lead C=1:75X 10-3 (rg= 1-1) and G=7-15 for 8. ~.1. 
For 2:54 em iron C=3-02 x 10~° (79=1-1) and Gt 3 for B ~ I. 
Apart from inherent uncertainties in assuming a Cross section given 
by a smoothed-out Born calculation, eqn. (A.1) is approximate only in 
that the factor 6-5/8 replaces a complicated function af & obtained when 
the conversion from total to projected angle is made. This function 
varies from slightly less than 1 at x=0 to § at large x 80 ones replacing 
it with an approximate average value introduces only a slight waueie 
Following Butler (1950) we now consider a single-scattering law we bus 
at ay into two parts A and B as shown in fig. 7, where v, is such that the 
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probability of a single collision giving a deflection «>, in one traversal 
of the scattering plate is small. Then to a good approximation the 
multiple distribution corresponding to this single scattering law can be 
written H(x)=H,(x)+H,(x), where H,(x) is the multiple distribution 
for A (i.e. expresses the probability of a resultant deflection x due to all 
single collisions (a ) and H,(#) expresses the probability of a deflection x 
as the result of one single collision greater than a, and all single collisions 
less than a). The distribution function H(x) is approximate in that it 
does not include the effect of more than one single collision greater 
than 2. Thus the smaller the probability of such a collision occurring 
the better will this approximation be. 

The multiple distribution due to A, for the case of the single scattering 
probability for a ‘ point ’ nucleus, is simply the distribution f(a, x9) given 
by Olbert. To a good approximation this is also the multiple distribution 
due to A for the single scattering probability for a U.N.M. given by 
eqn. (A.1) provided x, is chosen small enough so that these two single 
scattering probabilities do not differ by more than 20 or 30% at v=ap. 
It is easy to show that such a value of x, is given by 

1:3 

(36CG?)1/4 ° 
Therefore, if we choose a, as given by eqn. (A.2), then Olbert’s distribution 
will be a good approximation to H,(x), but at the same time a, must be 
so large that the probability of a deflection x>2, due to a single collision 
is small. In general both these conditions are satisfied simultaneously, 
e.g. for 2cm lead and 2-54¢m iron (ry9=1-1) eqn. (A.2) gives a~ 1 
and ~ 3 respectively and the probability of «>a, from eqn. (A.1) is 
less than 0-01 and 0-001 respectively. 

The function H,(x) gives the probability of a resultant deflection 
e=2'+x" due to one single collision x’>x, the probability of which is 
given by equ. (A.1), and a deflection x” due to all single collisions less 
than x», the probability of which is given by Olbert’s distribution. Thus 
HM ,(x) is simply the joint distribution of these two probability functions, 
namely, 


Cp 


(A.2) 


A (x)= 


f(a" S=2—2"', &)Po(x’, R) dz’, 
/ —@ 
where f(v"=x—vw’', a) is Olbert’s function, for the appropriate value of 
18. ye the variable replaced by x—a’, and Py(x’, R) is given by eqn. (A.1) 
for 2’ > | v)| and is zero for x’ < | x9 |. 
Finally we have 


A(x) du=f(x, x) dx+da | f(x" =2—z', x5) P(x’, R) dx’. (A.3) 


No attempt has been made to evaluate the integral in eqn. (A.3) 
analytically. Given Olbert’s distribution for the appropriate value of 


vy, H(x) may then be obtained for a particular case by a numerical 
computation. 
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In fig. 8 the distribution function H(x), for 2m lead and Pol ies 
compared with Moliére’s distribution for a ‘point’ nucleus and Olbert’s 
- distribution for r7=1-4 and 1:0. These distributions as shown have been 
multiplied by a factor 2 and thus refer to absolute deflections, i.e. positive 
and negative. It can be seen that Olbert’s distribution, even for the 
smaller value of 7, is only a fair approximation to H(x) at large x, and 
that the distribution obtained simply by adding the single-scattering 


Fig. 8 
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probability eqn. (A.1) ( é 
small-angle Gaussian approximation 
factor of 5 in the region x=3 to 4. vee —— 

In fig. 9 the distribution function H(«) is given for 2-54 Pa : sai 
yo—=1-1. Comparing this with the other distributions shown it 1s ee : 
that Olbert’s greatly underestimates the number of large deflections 
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(i.e. by a factor of 10 at ~=5) but, as in the case of the © point ’ nucleus 
distributions, a fair approximation to H(x) would be simply the Gaussian 
approximation at small a plus the single-scattering probability at large x. 
In general, the smaller the slope of the single-scattering curve the better 
is this approximate distribution given by the sum of the Gaussian and 
sin gle-scattering tail. 

Also shown in figs. 8 and 9 are the relative contributions due to inco- 
herent scattering for large energy transfers estimated for the theory of 
Gatto (1953). 
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Coulomb scattering distributions for 2-54 em of iron. 


Recently Cooper and Rainwater (1954 and private communication) 
using an extension of Moliére’s theory, have obtained a multiple scattering 
distribution which can also be applied to an arbitrary single scattering 
cross section. These authors give a distribution for 2.cm of lead for a 
single scattering law very similar to the one used here (for r7=1-1) and 
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hence a comparison is possible. The scattering variable used in bho 
which, for the case of 2 cm lead and 7,=1-1, is related to Olbert’s variable 
_£ by x=0-754/ho. To convert from the Cooper-Rainwater distribution 
to the one given in fig. 8 it is necessary to multiply their distribution by 
a factor 1/0-75 in order to renormalize and also by a factor 2 because the 
distributions in fig. 8 are for absolute angles. Because of slight differences 
between the single scattering law used by Cooper and Rainwater and 
eqn. (A.1), their distributions (single and multiple) at large angles have 
to be multiplied by a further factor 9/12 x 6-5/8 before being compared 
with those of fig. 8. The factor 9/12 arises from a difference in the 
assumed form factors (see Schiff 1953) and the factor 6-5/8 from the fact 
that Cooper and Rainwater tke the same single scattering law for projected 
as for total angle. Therefore, converting from ¢/¢, to x and multiplying 
their distribution (at large ¢/¢)) by the above factors gives the points 
shown in fig. 8. The agreement is seen to be excellent. 


REFERENCES 


Butter, 8. T., 1950, Proc. Phys. Soc. A, 63, 599. 

Converst, M. E., Pancornti, E., and Picctont, O., 1947, Phys. Rev., 71,209. 

Cooper, L. N., and Rainwater, J., 1954, Phys. Rev., 95, 1107. 

Firou, V. L., and Ratnwater, J., 1953, Phys. Rev., 92, 789. 

Gatto, R., 1953, Nuovo Cimento, 10, 1559. 

Grorae, E. P., and Evans, J., 1950, Proc. Phys. Soc. A, 63, 1248. 

HorstapTer, R., Fecuter, H. R., and McIntyre, J. A., 1951, Phys. Rev., 92, 
978; 1954, Ibid., 95, 512. 

McDrarmoip, [. B., 1954, Phil. Mag., 45, 933. 

Mo rere, G., 1948, Zeits. f. Naturforschung, 3a, 78. 

OLBERT, S., 1952, Phys. Rev., 87, 319. 

Rose, M. E.,' 1948, Phys. Rev.; 78, 279. 

Scuirr, L. I., 1943, Phys. Rev., 92, 988. 

Wuurremore, W. L., and Suutt, R. P., 1952, Phys. Rev., 88, 1312. 

YEnNnNIE, D. Re Witson, R. N., and RAvENHALL, D. G., Phys. Rev., 92, 1325. 


SER. 7, VOL. 46, NO. 373.—FEB. 1955 ) 


[ 190 ] 


XXIII. Collision Processes in Meteor Trails 


By H. 8S. W. Massey and D. W. Sipa 
University College, London* 


[Received November 22, 1954] 


SUMMARY 


Extensive data on the probability that an evaporated meteor atom 
will produce an electron in an ionizing collision before being brought to 
rest are now available. To obtain from these data values of the ionization 
cross section for impacts between meteor atoms and atmospheric molecules 
the momentum loss cross section for elastic impacts of this kind must be 
known. These latter cross sections are calculated for roughly represen- 
tative collisions in which the interaction between a meteor atom and a 
gas molecule is taken to be that between calcium and neon atoms calculated 
statistically. Values for the ionization cross sections are then derived 
and found to be consistent with the recent theoretical considerations of 
Bates and Massey (1954). 

The diffusion of ionization in a meteor trail is also discussed in relation 
to known observed data on the mobility of positive ions in gases. 


$1. INTRODUCTION 


THE study of the optical and ionization phenomena associated with the 
development of a meteor trail has progressed so far that it is possible to 
provide a quite detailed interpretation of the events involved in terms of 
certain atomic collision coefficients. The values of these coefficients 
averaged over the various meteor constituents may be obtained from an 
analysis of the meteor data and the question arises as to whether they are 
consistent with atomic collision theory and with direct measurement. 
The ionization produced by a meteor is due to the impact of atoms 
evaporated from the meteor with atmospheric atoms and molecules. 
From the meteor data it is possible to derive the chance that a meteor atom 
will produce an electron before it is brought to thermal equilibrium with 
the surrounding atmosphere. For comparison with collision theory and 
with direct experimental data it is necessary to relate this chance to the 
appropriate collision cross sections. The chance can be expressed roughly 
as @,/Qq where @; is the mean cross section for ionization and Q, for 
momentum loss (or diffusion) in a collision between a meteor atom and an 
atinospheric molecule. If Q, can be calculated Q; may be obtained and 
compared with expectation. Similar considerations apply to the emission 
of light by a meteor, the chance that a meteor atom will emit a quantum 
with frequency in the observed range being given by Q,,,/Q, where Q,., is 
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the mean cross section for an inelastic collision which leads ultimately to 
emission of an appropriate quantum. 

In this note we report the results of a calculation of Q, for a roughly 
representative collision, the impact of calcium atoms with neon atoms and 
discuss briefly the derived values of Q,. 

For the above application the cross section Q, is required for atoms 
with energies of 100 ev or more but the rate of decrease of ionization 
density in a meteor trail is determined largely by the rate of diffusion for 
positive ions with thermal velocities. These ions are mainly ionized 
meteor atoms and their diffusion coefficients may be estimated quite 
closely from observed data on the mobility of alkali metal ions in nitrogen. 
The comparison of these estimated values with those derived from meteor 
observations is discussed. 


§ 2. PROBABILITY OF IONIZATION oR Excrration By 4 Merror ATOM 
IN TERMS OF COLLISION CROSS SECTIONS 


Let J(H,6)dw be the differential cross section for an elastic collision 
between a meteor atom of kinetic energy # and a gas molecule, in which 
the direction of relative motion is turned through an angle 6. If the 
masses of the colliding systems are nearly equal the kinetic energy lost 
by the meteor atom as a result of the collision is given closely by 


4H(1—cos@). 
Hence in proceeding a distance 6s in an atmosphere containing n gas 
molecules per c.c. the amount of energy lost in elastic collisions is given by 


InESs i (1—cosO)I(E,0)dw=3nH@,(B)8s . . . . (A) 

where Qa(#)= | (1—cos 6)I(E,.0) dw, re a we (2) 
is usually known as the momentum loss cross section. We now have 

GF =— QB) —nBQ.—nB Qn - - - + (8) 


where Q,, Q,, are cross sections for ionization and excitation collisions 


, zs Th = 
between meteor atoms of energy H and gas molecules and #;, #,,, the 
mean energies lost in these respective collisions. If 54, is the number of 


electrons produced in ionizing collisions in passing the distance 6s we have 


aN, 
ae ms 
Using (3) we have : 
ae =—Q,/(4HQg+H Qi: tb caVex)- = Ms i ‘ Fs (4) 


This gives for 8, the chance that the atom will produce an electron before 
te 3 5 5 . : § ont 
its energy is reduced below the ionization threshold (2#;,), 


oh = ee) Kas 
B=| _ Q(B)dH/AEQa+ EQ: tbe) » - + - (5) 
J 28; 


OZ 


192 H. S. W. Massey and D. W. Sida on 


In practice the denominator can be put equal to }#Q, without serious 
error. This is because Q, is at least comparable with Q; and @,,, when 
ESE, or E,, and is much greater when H=E,. Actually @, will in 
general be small for values of # considerably greater than 2H;. If £, is 
the value of # below which Q, can be neglected we have approximately 


E 
pao | te ee 
|», QaE)E 

A similar expression holds for the chance of excitation. 

In practice, of course, there are many kinds of meteor atoms involved 
and at least two kinds of atmospheric molecules. We must then regard 
the cross sections Q;, Q, which appear in (6) as suitably weighted means 
over all the collisions of a particular kind which can occur. It is also 
assumed that f is less than unity so that consideration does not have 
to be given to the changes in the cross sections which result after ioniz- 
ation of the meteor atom occurs. The interpretation of ionization effects 
is essentially simpler than that of excitation for the ionization due to 
a meteor is directly observed whereas excitation can be only studied 
indirectly in terms of light emission. Many complicating factors are 
thereby introduced. 

To carry out a definite calculation we take calcium as a representative 
meteor atom. It is not possible to evaluate @, for collisions of calcium 
atoms with either N, or O, molecules. The detailed calculations were 
therefore carried out for collisions with an atom; that of neon. This 
has the advantage of possessing a closed shell of electrons for which 
the method of calculation is most suitable and an atomic number close 
to that of nitrogen and oxygen. It also has the advantage that experi- 
mental verification of the calculated cross sections is at least a possibility. 
In the formula (1) it was assumed that the colliding systems had approxi- 
mately the same mass. This feature will still be retained by taking as a 
representative collision one in which the interaction is that between a 
calcium and a neon atom but with the masses equal to that of calcium. 


§ 3. CaLcuLaTION OF Momentum Loss Cross SECTION 
It is well known that for encounters between systems of atomic mass 
classical mechanics is adequate exeept for very small angles of scattering. 
Owing to the weighting factor (l1—cos@) which appears in the formula (2) 
for @, these angles make a negligible contribution. We therefore have 
(Massey and Burhop 1952), when the masses of the two systems are equal, 


Pale Daa Dp 

ary: 
where 0= —2| —., 
ee ic 


2 22 / 
BIG 1 Se vy 


V(r) is the interaction energy and 7, is the largest positive zero of (7). 


Collision Processes in Meteor Trails 193 


The interaction energy between two atoms may be calculated to a good 
approximation by a statistical method due to Lenz (1932), Jensen (1932) 
and Gombas (1935). This treatment assumes that the electrostatic 
potential and electron charge density of the two systems are very little 
distorted by the interaction. The total potential and charge density are 
then obtained by superposition. The total potential energy is then given 
by the usual electrostatic formula and the total kinetic energy, which 
must include allowance for the increase due to the Pauli Principle, by the 
statistical formula for a Fermi gas 


4 / 3 \5/3 h2 
i= 5 (=) on Pi 


and a term correcting for exchange relations between electrons with 
parallel spins 
SS N\ue 
Lae ed awe 24/3 
ase 9 (=) pies 


m is the electron mass and p the electron charge density. 
In this way the interaction energy is obtained in the form 


SEE eff aee eet a — 5 | (Varvolra)+ Velrado(rade 


T, To 


Vir)= 


(477) 1/8 
5e5/8 


(Se) me | Elo ptr) (otra)? Entra} Ide 


% e ie ee | Cetra) + plrs)}*—{olr)}#9—{plra)} lar. 


2\87/ (47)18 
(8) 


4,2 are measured respectively from the nuclei of the as atoms, of 
respective charges Z,e, Z,e so that r;—r,=r. V,(r;), V (ra) are the 
respective electrostatic potentials due to the electron clouds in the two 
atoms, p(7',), p(7') the respective electron charge densities. ' 
This statistical formula is probably most satisfactory for use in 
discussing collisions between systems with closed shells. This is not quite 
- true in our case but the part played by the two outer electrons in the 
calcium atom will be very small when the kinetic energy of relative motion 
is as high as 100 ev as in the meteor encounters, and the effective inter- 
action is that between Ca++ and Ne for which (8) should be a good 
"OX ion. 
ean fields for Ca (Hartree and Hartree 1938) and Ne "hg 
1933) are available from which p and V, in (8) may be apm’, F peas 
these data V(r) was calculated from (8) using numerical SION, n 
is illustrated in fig. 1. Further numerical integration conld nee BS 
carried out to obtain I(H, 6) from (7). Calculations were carried out for 
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E=50, 100, 500, 1000 ev and the results are illustrated in fig. 2. Finally 
(, was obtained for these three energies, the values being given in table 1. 


Table 1. Momentum Loss Cross Section for Collisions between Ca and 


Ne Atoms 
E (actual) # (assuming masses equal) 
r Qa 
oe y (in units 79") 
50 75 4:39 
100 150 3°20 
500 750 1-34 
1000 1500 0-77 
Fig. 1 
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Internuclear distance in units a, 
Interaction energy between calcium and neon atoms, calculated by statistical 
methods. 
§4. Ionization PRopasiniry ror Mrerror AToMs 
If we assume initially that the ionization cross section for impacts 
between a meteor atom of kinetic energy H anda gas molecule is effectively 
independent of H for FE > E,, the threshold energy in (6), we have that 
8=20 E: dH 
“Je EQA(#)” 
Over the energy range of interest Q@, can be represented with sufficient 
accuracy by the form 
Qa=Vo(log B/C), 


where Q)=0-165 wa)? and C=12800ev. We find the following values 
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of Q,/8 for different initial energies and tw 
Wie aiccaay 2. 2 o reasonable values for the 


Table 2 
Initial energy of meteor atom Q;/B in units 7a)? 
ev H,=20 ev H,=10 ev 
1000 0-35 0-33 
500 0-54 0-49 
100 1-61 1-25 


It will be seen that, to be consistent with our assumption that @; 1s 
effectively constant over the energy range from E, to 1000 ev the 
probability of ionization 6 should increase with the velocity v of the 
meteor atoms as v" where n=—1-3 or 1-2 according as #, is taken as 20 or 
10 ev. The most recent observed data (Kaiser 1954) give n=0-5--0:5. 
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Angle of scattering (C.M. System) in radians 


Differential cross sections for elastic collisions between calcium and neon atoms. 
Curves I, II, III, IV refer to collisions in which the neon atoms are at 
rest and the calcium atoms have kinetic energies of 50, 100, 500 and 1000 


ev respectively. 


This suggests that it is likely that Q; increases somewhat as v decreases. 
According to the theoretical considerations of Bates and Massey (1954) Q, 
varies approximately as v1 over this energy range so that the calculation 
of 8 was repeated with this assumption. It was then found that »—0-87 
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(E ,=20 ev) or 0-57 (H,=10 ev), either of which is certainly not in conflict. 
with the present observed data. 

Assuming that the theoretical variation of Q@,; with v is correct ie 
value of Q, comes out to be Q,°(1000/H)"? Byzay? where @;°=0-35 or 0-33 
according as H,=20 or 10 ev respectively. is the kinetic energy of the 
meteor atom in ev and f, is the value of the ionization probability for 
meteor atoms of 1000 ev energy. £, is not known very accurately but a 
value of 0-2 seems to be a reasonable estimate. This gives absolute values 
of Q,° which fall within the range of values expected theoretically. Thus, 
according to Bates and Massey (1954), the ionization cross section due 
to impact of atoms A and B is given approximately by 


Q,=(3-8 X 10-5 g MYR 3/BYt) may? 


where g is a statistical weight factor of order unity, / is the reduced mass 
of the colliding systems (with the proton mass as unit), R, is the 
separation at which the potential energy curve along which A and B 
initially approach intersects a suitable one for interaction of A and B+, 
FE is the energy of relative motion in ev and t is the lifetime in seconds for 
auto-ionization of the complex AB when the separation is less than R,. 
Values of t between 10-18 and 10-1 sec are reasonable. Assuming these 
and taking g as unity we find that the meteor data require that R,, should 
lie in the range 0-3a) to 2-0a), which is not unreasonable. 

We may therefore summarize the situation as one in which the theory 
is consistent with the observations but there is a great deal of flexibility 
in the comparison which can only be removed by laboratory experiments 
on the ionization of atoms by atom impact. They are very difficult to 
carry out but would seem to be worthwhile. 


$5. AMBreOLAR DIFFUSION COEFFICIENT IN A Merror TRAIL 


The rate of decay of the radio echo from a meteor trial gives a measure 
of the rate of decrease of electron concentration in the trail. This decay 
is due to ambipolar diffusion of ions and electrons, to recombination and 
to electron attachment. Analysis of the evidence suggests that the latter 
modes of decay are usually unimportant (Kaiser 1953), the main loss 
being by ambipolar diffusion. From a study of the decay rates of echo 
from the Arietid meteors Kaiser finds that the most probable value of 
the ambipolar coefficient D at the height concerned (92-5 km), where the 
particle density is 7 x 101%/c.c¢., is 2:3 x 104 ¢.g.s. units. 

An ionized meteor atom of mass M/ will come into thermal equilibrium 
with atmospheric molecules in a time ¢ of order 


(2M)V2 -& dB 
nm Jy, BPQA{(E) 


where #;, is the initial kinetic energy of the ion and £, the mean thermal 
energy of a gas molecule. Substitution of calculated wale of Y, shows 
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that this time is between 10-4 and 10-3 sec at 92-5 km altitude. For 
electrons the corresponding time is of order 


M (1 \'2 Ft dH 
8m z, E820) (HB) ” 


n 
where Q,° is the momentum loss cross section for encounters of electrons 
with gas molecules. This may be taken as about 10-15 em? and is very 
roughly independent of energy. If H, SE, we find that a time of about 
0-05 sec is necessary for the electrons to come to thermal equilibrium at 
92-5 km altitude. This is not far below the echo decay times observed 
for the Arietid meteors but it is probably quite a good approximation to 
assume that both the electrons and ions have come to thermal equilibrium 
with the surrounding atmosphere. Under these conditions and, provided 
there is no magnetic field present, the ambipolar diffusion coefficient 
D=2D,, where D,is the diffusion coefficient of the positive ions. When 
such a field is present the situation is more complicated. If fluctuations 
in ion and electron concentration are ignored the diffusion of electrons 
transverse to the magnetic field of strength H is reduced by the factor 
v."/(v,.?-+-w?) where v, is the electron collision frequency and w=eH/mc. It 
has been pointed out by Kaiser (1953) that the effect of this factor would 
be to reduce the electron diffusion coefficient D, transverse to the field to 
a value comparable with D, at a height of about 100 km. When D,=D, 
the ambipolar diffusion coefficient D becomes equal to D, instead of 2D,. 
We therefore expect that for the Arietid meteors the mean ambipolar 
coefficient D should be between D,; and 2D. 

D, can be written in terms of the ionic mobility /,; by the relation 


page! 
Dat. ([), 


where x is Boltzmann’s constant. The mobilities of alkali ions in nitrogen 
of such purity that no clustering takes place have been measured by 
Powell and Brata (1932). It is unlikely that the mobility in oxygen 
will be very different. Using the observed values we may determine D; 
for the conditions prevailing for the Arietid echoes. Allowance may be 
made for the low temperature (205°K) at 92-5 km as the ion diffusion 
coefficient may be assumed to be proportional to the absolute tempera- 
ture for constant gas density. We find for Nat, Kt, Rb* and Cs* 
ions respectively that D,=2-37, 2-10, 1-86 and 1:383><10% C.g.8. units 
respectively. In view of the fact that for singly charged ions in the 
same gas the mobility varies approximately as (1+ M,/M,)"?, where M, is 
the mass of a gas molecule and I, of an ion, it is to be expected that: this 
range of values of D; includes all those likely to arise in a meteor trail. 
The agreement with the mean value derived from the Arietid meteor 
echo decay, which is between 1:15 and 2:3x 104 ¢.g.s., can therefore be 
regarded as reasonably satisfactory. 


ys 
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XXIV. Oscillatory Motion of a Viscous Liquid in a Thin- Walled 
Elastic Tube—I : The Linear Approximation for Long Waves 


By J. R. Womersitey 
Physiology Dept., St. Bartholomew’s Medical College* 


[Revised MS. received October 19, 1954] 


§ 1. INTRODUCTION AND SUMMARY 


THe problem of determining the motion of a liquid in an elastic tube 
when subjected to a pressure-gradient which is a periodic function of the 
time arises in connection with the flow of blood in the larger arteries 
(Helps and McDonald 1954, Womersley 1954). Attempts have been made 
in the past to measure the rate of flow in the aorta and femoral artery 
of the dog and rabbit (Shipley, Gregg and Schroeder 1943) and to relate 
these observations to the varying pressure. In the absence of an adequate 
mathematical theory, these were not very successful. More recent 
determinations by direct observation of the motion through the trans- 
lucent arterial wall, using high-speed cinematography (Helps and 
McDonald 1954, Womersley 1954) have been accompanied by measure- 
ments, not only of pulse-pressure, but of pressure-gradient. Fair 
agreement has been shown between the observed rates of flow and a simple 
solution for oscillatory motion of a viscous liquid in a tube with rigid 
walls (Womersley, in press). 

In this paper the corresponding solution for a thin-walled elastic tube 
is given, it being assumed that the effect of the inertia terms in the 
equations of viscous fluid motion can be neglected. An approximate 
correction for the effect of the inertia terms will be presented in Part IT 
of this communication. 

It is shown that, when the liquid contained in the tube is viscous, 
the pressure-wave cannot be propagated without distortion. Not only 
is the motion damped, but the wave-velocity rises as the frequency 
increases. For constant frequency, the wave-velocity rises as the 
viscosity of the liquid decreases, tending to an asymptotic value equal 
to that for a perfect fluid given by Lamb (1898). It is also shown that 
the longitudinal oscillation of the walls of the tube (caused by the viscous 
drag on its inner surface) is important in determining the rate of flow, 
which may be 10°% greater than that in a rigid tube under the saine 
pressure-gradient. 

0 es ee 
* Communicated by Sir Geoffrey Taylor 
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§ 2. THE SOLUTION FOR THE RiciIp TUBE 


The simple solution for the oscillatory motion of a viscous liquid in a 
rigid tube, under a simple-harmonic pressure gradient, was given by 
Lambossy (1952) who gave the formulae for velocity and viscous drag. 
He did not apply his result to the motion of the blood in arteries, being 
solely concerned with the effect of the viscous drag on the frequency 
response of pressure recording instruments. The author obtained the 
same result independently, in a different form, and derived the expression 
for the rate of flow. This result has been used to predict rates of flow 
from observed pressure-gradient (Womersley, in press). This solution is 
repeated here for completeness. Let R be the radius of the tube, w the 
velocity along the tube, Ae‘”’ the pressure-gradient, 4. the viscosity of 
the liquid, p the density of the liquid, v=/p the kinematic viscosity. 

The equation of motion is 

OF a Owe a Lice ae 
PTR OE gamcnintere tat t he Ce a we (1) 

Let y=7/R and w=ue'"! and let the non-dimensional quantity R/(n/v) 

be denoted by «. The equation for w is 


Ae LGUs wo) e cuca Ie 3 
dye ydy (2) 
Ake! J 9(at?/2y) 

and therefore w= Set mat a et eint (3) 


In Lambossy’s paper he separates this into real and imaginary parts 
in terms of ber and bet functions. He also replaces R?/x2 by its equivalent 
vin. The multiplying constant in (3) then has the simpler form A/inp, 
which is convenient for calculation, but does not show the relationship 
of the result to that for steady flow. It is more instructive to express 
the results in terms of modulus and phase, rather than in real and 
imaginary parts. Following the notation of McLachlan (1934) we write 


J o(xi?/?)—= M(x) exp 76, (a) 
and therefore (3) above may be written in the following form : 
If the real part of Ae'"’ is M cos (nt—d) then 
ARM. 


pe re 


w= — 


sin (nt—d+€,) er ice pg! 


where J/,’ and e¢, are defined as follows : 


We have es = M, (ay) 
Toloal*y)/ J (ni8?) = PAR exp i {0g(cy) —09(2)} 
Let M,(ay)/M (x) be written ho, and write 89 =F (%)—Oo(ay). Then 
My'=/(1-tho§—2hy 6089), , - ss ee 


tan ep=h sin 85/(1—hy cosdg). . . . . . . (6) 
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The rate of flow, Q, is given by 
“1 
Q=27R? | wy dy gr ogee wee mm Od) 
Jo 


and integration of (3) gives at once 


TEL eA 2S ,(at3/2) 

Qa sal! game oO (8) 

so that for a pressure-gradient M cos (nft—¢) the rate of flow is 

7M Rt My S. 
Q= —— .—} sin (nt—d+e,)) N aiee: a6 G) 
fe va 
where, if 2M j/aMy=hyo, 5y9= (37/4) —0, (0) +0,(c), 

M yy = V (lL +h y9?—2hyy 008 845), | 
C10) 


tan €49 =Ayo SiN 849/(1—hyo Cos 849). j 


The quantities M,)'/«? and €,) have been tabulated by the Mathematics 
Division, National Physical Laboratory, for the range 0<«<10, and are 
published (Womersley, in press). 

We make the following observations on this formula. First, the 
importance of the non-dimensional parameter «. This determines the 
character of the motion, and seems to vary very little (at corresponding 
points) in different animals. For example, in the femoral artery of the 
dog, which has a diameter of 3mm, the pulse frequency being 3 per 
second, the value of « is 3-34. In the corresponding artery in the human 
thigh, the larger diameter is offset by a lower pulse-rate, and the value 
of « is about the same. The variation in phase-lag across the tube 
(somewhat analogous with the ‘ skin-effect ’ in electrical conductors) is 
considerable even at normal pulse-frequencies. Velocity-profiles for a 
range of values of « have been given by Lambossy (1952). The effect 
of « on the rate of flow is shown in fig. 1. The graph of phase-lag against 
x is self-explanatory. The graph of Qynax/Qsteaay gives the ratio of the 
modulus of the flow (i.e. its maximum value in either direction) to that 
of the Poiseuille flow for the same pressure-gradient. The values of « 
which are of interest in the study of the motion of the blood in the large 
arteries come on the steeply-falling parts of both curves. It is of interest 
to note also that both flow and phase-lag approach their asymptotic 
values quite slowly. For values of «>10 there are simple asymptotic 
expressions for both My,)'/x? and €y which can be derived from the 
asymptotic expansions for the Bessel F unctions. They are 


M9 a2=1fo?—2/a3+/2+ 1/a*, 
Ga = 2/0 4/ 2+ Liga 19/244/2a?. 
If the observed pressure-gradient is expressed as a Fourier Series, 


: RAE annroc 
then the rate of flow can be calculated term by term, using the appro 
priate value of « for each harmonic, and the result compared with the 
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observed rate of flow. This has been done, with fair agreement with 
observation (Womersley, in press) when one recalls the drastic nature of the 
assumption that the artery is a rigid tube, and that the formula contains 
no disposable constants. 


Fig. | 
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§ 3. THe Equations or Morron or THE Evastic TUBE AND THE 
DERIVATION OF THE PuLSE VELOCITY 


In Lamb’s paper (1898) “On the Velocity of Sound in a Tube as 
Influenced by the Elasticity of the Walls ”’, he considered a compressible 
inviscid fluid in a thin-walled elastic tube, took a wave-equation for the 
pressure inside and obtained a ‘ frequency equation’ for the velocity in 
terms of the dimensions of the tube and the elastic constants of the 
material. For a gas in a metal tube the correction is negligible, but he 
remarks “ for the other extreme, e.g. water in an india-rubber tube, we 
have two values of c, (hH/2ap,)"? and (B/p)¥2”. In these expressions h 
is the thickness of the tube, a its radius, and p its density. H is Young’s 
Modulus, and B is H/(l—o?) where o is Poisson’s Ratio. Po is the 
density of the liquid: 
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Lamb adds that the former value of c is the same as the velocity of the 
pulse-wave as deduced theoretically by Résal (quoted by Korteweg in 
_ Liouville’s Journal in 1876). It is apparent from Lamb’s results that if 
the wavelength is long compared with the circumference of the tube, the 
velocity is independent of frequency. We show below that this does not 
hold if the velocity of the liquid is taken into account. Let R be the radius 
of the tube, h its thickness, p the density of the material, p, the density of 
the liquid. Take the axis of z along the tube, and let ¢ be the displace- 
ment along the axis, € be the radial displacement. For the elastic constants 
take B to be H/(1l—o?) as above, where # is Young’s modulus and c 
Poisson’s Ratio. Then the extensions are 


a6 along z, and & radial. 


z R 
The tensions are P along z, and Q radial, so that 
PaBh{S +oqh. en een Lee 
Q=Bh {5 +055 ees kee ee 2) 
Lamb gives the equations of motion 
ph = a Zz i Coe wee eee eee Ls) 
phot =p—% Sie aa SEE) 


where p is the excess pressure on the inner surface, but if the liquid is 
viscous there will be a surface traction on the inner surface of the tube due 
to viscous drag, and eqn. (13) must be modified. If w is the longitudinal 
component of the velocity of the liquid and w its radial component, the 
surface traction per unit area will be 


ow dou * 
— + — BAG e). Soe 15 
bb ( a: =), (15) 


where jx is the viscosity of the liquid. Making this modification to eqn. 
(13) and inserting the values of P and Q from (11) and (12) the equations of 
motion of the tube become 


as y 2 0& 
ee (RS) +elgetRe) «09 
PrimmER oy 0) 02). pee Bee 
gees t “(az z) int ca Sele ots (17) 
pe poe oR Oe RP a 
together with the boundary conditions for the motion of the liquid 
U= a, AeA yc ek (138) 
w= gs Atel ee mes aes = (C18) 


ot 
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In bis ees notation the equations of motion of the liquid are 


Ou 1 Op ou. lou Cu “wu 
rites te a= Soa. OF Eee aa] (20) 
ss ss CU al cn C7 OU eee © 
toe TY eG, be el ra oe 
together with the equation of continuity 
Ou uu. dw 
ne Pe pe 2 
a ts bap 0 iy ee ee 


If a trial solution is taken in which all the variables vary as 
exp [in (t—z/c)], we shall be concerned with the type of motion in which 
u/c, w/e and nR/c are all small. , It is clear from the equation of continuity 
that w, the radial component of velocity, will be of order nR/c compared 
with w, the velocity along the tube. It follows that the non-linear terms 
are in general of order 1/c compared with the main linear terms, and we 
therefore neglect the inertia terms in the first approximation. Let 


p=P, exp [in (t—z/c)], 

Uu=uU, exp [in (t—z2/c)], 

wW=W, exp [in (t—z/c)]. 
Equations (20) and (21) reduce to 


Chase a ena u, 1 FR dp, 

a Sched CALS AA 3 * 
igi? ra PP Re 2 ove weed eras 
dw, Idw, 2 ey 
dy y dy ae ow, Pe asia i eee 


where the terms in 0?w/dz?, @u/dz* have been omitted, since n?R2/c? is 
small (see below). 
The equation of continuity is 
poe a _ ink D5 
int jini of pe ROD a oan eee (25) 
If now it is assumed that p,=A/J (ky), where & is to be determined, the 
equations of motion can be integrated to give 


J y(a03/2y) ink® A, 


w=Cy Tae?) cu Boe e_fpolky), AA) 
, Ji(ar?y) RB A, 
U1 =C a> aay a?) a” ae ol (ky); sae ee a 


where C, and C, are arbitrary constants, and the constant J 9(«23/2) has 
been introduced (by analogy with the simple theory for the rigid tube) for 
convenience. If these values of u, and w, are inserted in the equation of 
continuity it should reduce to an ern, From (26) and (27 ) 

inR mR  So(at?/2y) 72n2R3 A 

pag his ae C Oe ape J (og3/2) (a 13/2) i Cun ae) —}2 JT (ky), . (28) 

ld J (at?/2y) Rk? A 

"3/2 (1 peed s 
wage yar? Os RACH lie ge ts yz pp olky). -  sarteey 
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For the equation of continuity to be satisfied these must be identical, and 
therefore k=inR/c and C,/C,=inR/cui?/2, so that J,(ky) becomes I o(rR/c) 
_and J,(ky) becomes iJ ,(nR/c). 

At the highest frequencies likely to be of interest in the pulse-wave, 
n* R?/c? is small compared with «2. As an example, for the sixth harmonic of 
the pulse-frequency of the dog, n?R2/c2~6 x 10-4, and the ratio of this to the 
corresponding value of «? is about 9x 10-°. It is therefore justifiable to 
replace i?«?—k? in (26) and (27) by 73a2, and in view of the smallness of 
n/c to use the approximations J,(nR/c)=1 and I,(nRy/c)=nRy/2c. This 
is the same degree of approximation that is implicit in omitting the terms 
Pw/dz and d?u/dx? from the equations of motion. Inserting these 
approximations in the expressions for the velocity-components, _ 


ad g(ae?/*y) AY 


=f | 30 
Wy C7 J 9(a03/2) I PoC ? ( ) 
ink 2S (a23/2y) As 
= — = a Sn |e fe (31 
“I= Fe [ 1 BPI ,(ai82) 4" oe aD 
At the inner surface of the tube, i.e. when y=1, 
Ay 
AG ——— . . . . : . . . . . . 32 
wy=Cy+ ee (32) 
dD Dib ise aye. links : 
ea aa pn Ve 5 . . 3 (33) 


where F',,.(x) has been written for 2J ,(ai?/*)/a1?/2.J9(a1?), the function of 
«% which appears in the expression for the rate of flow in the simple 
theory. The value of dw,/dy at y=1 is also needed, to substitute in the 
equations of motion of the tube. Making the same approximations, this is 


DaaA n?R? A, RA 
Se oO) ns ae ae (34) 


If it is now assumed that é=D, exp [in (t—2/c)], 


C=—E, exp [in (t—2/c)], 
where D, and Z, are arbitrary constants, the boundary conditions for uw, 
and w, become 


ae wee (35 
un Bia Oq Pol 3 . : . Q A o ( ) 
; link A, 


and the equations of motion of the tube become 


Ai Bio*( ine al $ BF) 
See) 
Coe. L Aiea 
Athy = — 28) — SP. iat Flt) + 2 pee C2 
Bie nt Jo [gee e138) 
ee +l s)Py 
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In (38) the term in du/dz in (16) has been omitted, since it is of order 
nh? /c?, | 

Equations (35)-(38) are four homogeneous equations in the arbitrary 
constants A,, C,, D,, H,. Eliminating them will give a ‘ frequency 
equation ’ which will determine the wave-velocity c in terms of the elastic 
properties of the tube and the non-dimensional parameter «. The result 
of the elimination is 


0= 1/po¢ 1 0 —in 
inR/2poc? inR/ 2c —in 0 

Ihp 0 n2—B/pR? ino B/cRp 

ee ee ee 

2° p hR’ pec’ Ty CK p pc 


The first step in the evaluation of this determinant is to make all the 
elements non-dimensional. The term —ipovx? F'yo/2phR is replaced by its 
equivalent —inp,RF 1 9/2ph, since va?=uR?. If now the first column is 
divided by poc, the second column by inR/c, fourth column by —in, and 
the fourth row by in, the determinant becomes 


i 1 0 1 

\ Fy zi 0 
Rpolhp 0 — B/pc® —oB/pe* 

i) —F ,,Rpo/2hp —oB/ pc? 1— B/ pc? 


At this stage the following further approximations have been made. In 
the third row, third term, ?R?/c? has been neglected in comparison with 
B/pc?. This is justified, since (B/p)'/? is the velocity of shear-waves in the 
material of the tube, and this will be greater than the pulse-velocity, so 
that B/pc?>1, and n?R?/c? can be neglected in comparison with it. The 
first term of the fourth row is of order n2R2/c2x2 and can therefore be 
neglected. If now we write k=hp/Rpo, x=kB/p,?, and reduce the order of 
the determinant by replacing the second row by (Row 1—2 Row 2) and 
the third row by (Row 1— Row 3) the equation becomes 


0=) 1—F 5, 2 1 
i a l+ox 
—if —or k—a 
which reduces to eae ( 
In this equation /<1 since the tube is assumed to be thin-walled, and in 
the application to the arteries p= py. 
When « — oo, i.e. as the viscosity of the liquid +0, F 19 > O and the roots 


of this equation tend to equality with those of the corresponding equation 
which can be deduced from Lamb’s equation for the compressible fluid. 
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An equation which is essentially the same as (39) has been given by 
Morgan and Kiely (1954) in a somewhat different notation. Thev give 
. two approximate solutions, one for * small viscosity ’, by using the asymp- 
totic expansions of the Bessel Functions, and the other for ‘large 
viscosity ’, by using the power series, but do not seem to have appreciated 
the dependence of the wave-velocity on the single non-dimensional para- 
meter, «. 
The roots of (39) are given by 


(l—o*)a=G4/{[G—(l—o?)H} .. . .°. . . (40) 
where 
1+1—o 
G= Tes +(k/2--o—4h), 2. <2 2 (41) 
ei TOP he Panam tay ae oe vee 8 (42) 


In terms of the notation used in (9) above, 1/(1—F'4))=exp (—te)/M 49’ so 
that the quantities required to compute the roots are already available. 
Since F(x) is complex, x is always complex, and the motion is either 
damped or unstable. The sign of arg w is determined by that of arg G, and 
since arg G is always negative, the motion is damped. If we write 


(l—o)e/2e=aX—iY . . J... . (48) 


and denote by c, the velocity for the perfect fluid, ie., write 
Co=(hE/2Rp,)2, then if c, is the wave-velocity, Co/¢y=X and over a 
distance of one wavelength the amplitude will be reduced in the ratio 
exp (—27Y/X). The quantities of practical interest are, therefore, Ise. 
which gives the ratio of the velocity to that for the perfect fluid, and 
exp (—27Y/X). It is hoped to publish full tables of these, together with 
other quantities required for calculating the longitudinal velocity and rate 
of flow of the liquid, in a later communication. Table 1 gives their values, 
together with the reduction in amplitude over a distance of 10 cm, for the 
first four harmonics of the pulse of the dog, in the femoral artery. They 
have been calculated for k=0:1, i.e., for a wall-thickness of one-tenth of the 
radius, for o=} and o=0. In this table, f;=exp (—27Y/X) and 
fo=exp (—27Y2/XA), with z=10 cm. 


Table 1 
o=0 
c/Cg | hi fo | 
0-842 0-132 | 0-917 | 
0:876 0-294. 0-900 | 
0-894. | 0-381 | 0-883 | 
0-906 | 0-442 | 0870 
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The amount of damping seems rather large. A diminution in amplitude 
of 56°, over so short a length of 10 cm one would expect to have been 
observed and remarked upon, but until more accurate observations are 
available it is not possible to say with certainty that this degree of damping 
is greater than that which exists in the arterial system. In practical 
observation it might well be masked by the change in shape of the pulse as 
it travels. Figure 3, which shows the relationship between « and 
exp (—27Y/X), illustrates how much greater the damping is for the smaller 
values of co. Figure 2, which shows 1/X =c/cy against «, demonstrates that 


1-O7 


ee ee 
X= ORT 


arma 


2 i) 4 5 6 fi 8 9 10 


a 
for values of « greater than 3, which are those of greatest practical interest, 
the change in velocity with « is quite small. It is somewhat larger the 
smaller is the value of o. Since the wave-velocity varies with frequency, 
energy will be transmitted with the group-velocity, and it is the group- 
velocity that will be observed when the pulse-velocity is measured. If a 
wave-motion is described by the form exp [i(nt—mz)], the group velocity 
c, 18 given by c,=dn/dm. In terms of « and X this may be written 

Co + € a dX 
= 21-3 F) pt ee ae 
‘g 1 Aes LU 
The analytical form for dX/d« is unsuitable for computation. Figure 4 
shows a graph of log X against log «, and since 
adX  d(log X) 


X dx d(log a) Ge. 


2ry 
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it is possible from this to estimate its magnitude. It’will be seen from the 
graph that for those values of « which apply to the femoral artery 
(3<«<7), d(log X)/d(log ~) is about —0-045, so that for these conditions 
c,/¢, is about 0-98, so that the difference is only 2%, and over the range 
x=3 to «=4 is certainly never more than 2}°%. Until accurate measure- 
ments of pulse-velocity are made over short lengths of artery, this effect 
is not likely to be worth taking into account. 


§ 4. THe Motion OF THE LIQUID 
In this section the details of the motion of the liquid at a particular 
value of z, i.e., over a short length of artery, will be studied. If the origin 
of z is taken at the mid-point of this short length, the longitudinal velocity 


A 1 J 9(a13/2y) int ; 
w= {= +0, BAC e : ; 5 3 - (46) 


where the value of C, is to.be determined from the boundary-conditions. 

From (35), (36) and (37), D, and #, can be eliminated. If, in the resulting 

equation, 7 is written for the ratio of C, to A,/poc, the value of 7 is given by 
2 (1— 20) 


a Pag=2a\or haba (40 
where z is the root of eqn. (39) and then 
Arf, Folat**y)) sn 


The first point to be noted is that for a given pressure-function, the 
pressure-gradient is inversely proportional to the wave-velocity, so that 
the velocity of the liquid, and therefore the amount of flow, will be 
proportionately greater in that which is more elastic. The constant, 


A,/po¢ may be written ) 1 ( 1 ) 
Fy aoe ce Te ee ae A ii 
Pg A el in 


and since «?=nR?/v this may be written 

A] ppt Ay Rate ey ey eee oe 
where A,’ is the quantity that will be given by Fourier Analysis of an 
observed pressure-gradient. This establishes the identity of the constant 


in (49) with that in the simple theory for the rigid tube. The average 
velocity across the tube is 


w= 


A 
={linF pet. 5 1.4 2 wee fen 
Pole 


This can be put in the form 
mee Age a aaa 
pe ee 


eo 


sin (né—d-+-ey)') ss. 


for a pressure-gradient M cos (nt—¢), and the values of M 10 and e459. 
compared with those of M,9’ and ¢4) for the rigid tube. For c=} and 
o=0 the expressions for 7 and 7 Fy are particularly simple. 
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When o=3, n= —2/a(1—F,), 


Spite — Wet e(l= Fis). a's. 5! (82) 
When o=0, n=(1—2/x)/F 4, 

eae oe are te eS (53) 
Using these formulae, and the corresponding one for c=}, the values of 
M9" and €9” have been calculated for c=}, t, 0, for the first four harmonies 
in the femoral artery. The results are given in table 2. It will be seen 
that if the wall is ‘ perfectly elastic ’, i.e., suffers no volume change in 
stretching, the average velocity is some 10° more than that for a rigid 
tube under the same pressure-gradient. At the other extreme, if the wall 
is composed of a material which does not form a ‘ waist ’ on stretching, the 
velocity is slightly less than that for arigid tube. The phase-lag for o=tis 
greater than that for the rigid tube, but for the smaller values of a is less, 

and this difference in phase-lag becomes smaller as « increases. 


Table 2 


32-11° | 0-6379 
4-72 0-7945 | 17-86° | 0-7767 | 20-29° | 0-7378 | 22-38° | 0-7436 | 19-95° 
5:78 0-8275 | 13-71° | 0-8117 | 15-77° | 0-7807 | 17-48° | 0-7839 | 15-82° 
6-67 0-8471 | 11-59° | 0-8279 | 13-58° | 0-8076 | 14-80° | 0-8096 | 13-49° 
| 


Table 3 
o=3 o—0 
Amplitude Phase- Amplitude Phase- 

os ratio difference ratio difference 

: i 0-122 74:0° 0-527 26:2° 
Ne 0-166 52-9° 0-226 16-2" 
5:78 0-190 38-0° 0-145 12-6" 
6:67 0-257 28-6° 0-104 10-5 


It is also of interest to calculate the longitudinal velocity at the wall. 

is J A | 
This is Aa an) oft 
Pof 
so that its ratio to the average velocity is (1-+)/(1-|-4# 10): i 
and phase of this are calculated, they will demonstrate the ratio of oi 
magnitude of the velocity at the wall to that of the average velocity, and 
the phase-difference between the two. Table 3 shows these quantities, 
for c=}, and c=0, for the same values of « as in tables | and 2. 


2 


Tf the modulus 
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These figures emphasize the critical role of Poisson’s Ratio in deter- 
mining the details of the motion. The phase-differences in the above 
table are shown with their correct sign in relation to the amplitude ratio, 
and therefore show that the velocity at the wall ‘leads’ the average 
velocity. The most striking point about this effect is its magnitude, 
which is greater than might be expected. It is not possible to estimate 
it from the experimental results available, but it would seem that the 
attempt to find experimental means of measuring the effect would be 
worth while, since it would form a critical test of the theory, and would 
throw light on the elasie properties of the arterial wall under dynamic 
conditions. 

The formula for the rate of flow might at first sight be written 


A 
Q= — 7 R(1+9F yo)e™ Thoth ee eee 


Po 
but this is approximate only, since at any time the value of the radius 
is not R, but R+€, and € varies with the time. A better approximation 
should be given, therefore, by 
2 ; 
= oR (1+ =) (Layers al oe 
Now 2é/R=w/c, a useful formula for rough use in estimating whether 
the arterial expansion observed in an experiment is consistent with the 
observations of pressure-gradient and average velocity. Even (55) is 
not fully corrected for variations in radius. In the original equations, 
& has been regarded as constant from the point of view of motion of the 
liquid. In fact the boundary of the liquid is at R+€. If we replace 
Rk by R-+-€ in (20) and assume that w is still essentially a function of y 
(i.e. that the lines of laminar flow expand and contract with the artery) 
(20) may be written 
Ow ldw R? se 2E\ dw ; 2€\ R? dp 
oy? Ty ey iv ( ) B) a = ( =m ae 
the inertia terms and the term in 0?w/dz? being omitted. 

If we seek a solution which is of the same form as that for a constant 
radius, we can imagine 0p/dz represented by a Fourier Series in n(¢—z/c), and 
corresponding Fourier Series for w and 2¢/R (which can be represented 
in terms of w) substituted in the equation. The products of the Fourier 
Series can be multiplied out, and a set of equations for the velocity 
components obtained by collecting up corresponding terms. The wave- 
velocity is not the same at all frequencies, however, and on multiplying 
two periodic terms together there will be exponential factors in z ‘ left 
over’ as it were, which would disappear (being equal to unity) in a 
system with constant wave-velocity. A typical term is 


(56) 


exp {in(t—z/c,)} x exp {in(t—z/c,)} 
which has to be compared with exp {2in(t—z/c,)} where c, is the wave- 
velocity for a frequency n/27, and c, the wave-velocity for double this 
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frequency. The question is, therefore, how far does exp{2in(1/c,—1/c,)z} 
differ from unity ? A glance at the graph of c/cy against « shows that 

_for x>3 the variation in c is small, and in fact for the first four harmonics 
of the pulse of the dog, all these ‘ stray ’ exponential terms have wave- 
lengths of ten metres or more. 

In the next section, therefore, a simple perturbation method, along 
the lines described, is used to solve (56). It is admitted that this is not 
necessarily the most important correction to be made to the linear solution 
given above. But it is the simplest to do, and provides a ‘lead-in’ to 
the more important correction for the inertia terms. And though at 
first sight less important in principle than the inertia-term correction, it 
may well be equally important in magnitude, for it may cope with fairly 
large arterial distentions, such as occur near the heart, without the 
mathematical difficulties that arise in the analysis of finite strain. 


§ 5. Frrst-ORDER CORRECTION FOR VARIATION IN DIAMETER 

A solution is sought of eqn. (56) when the pressure-gradient is a Fourier 
Series of four harmonics, together with a constant term which will be 
assumed to give a Poiseuille flow, the static expansion of the tube (which 
would give a tapering effect) being neglected. This is justified, since 
this constant term is small, being less than one-eighth of the largest 
oscillatory terms. The detailed solution will be developed for two 
harmonics only, since this illustrates the method adequately without 
wearisome repetition, but the numerical example quoted has been worked 
out for a series with four harmonics. Since (56) is non-linear, we may 
no longer write the pressure-gradient as 


AA, OFA 5 Ota. oe 
and take the real part, or half the interaction terms will be lost. It is 


necessary to start from the pressure-gradient in real form, and write 
down its exact complex equivalent. Assume that 


Op/dz—=M,+M, cos (nt-+¢1)+M, cos (2nt+o). - ~ (57) 
If we define Ay, A,, Az by 
A,=M,, Ay=3M, exp (*4:); A,.=3M, exp (tp) 
apfdz—=A, +A, e™*+A,*e "+A, ernie, (88) 


using a standard notation for conjugate complex quantities. Taking the 
value of w already obtained as an adequate approximation to use in 
substituting for 2é/R in eqn. (56) we have 


then 


142¢/R=1+ = (Gy eft C,* ei, e%!1O,* emt) (59) 
i) 


: Av? a 
in which O,=(X,—1¥) $5 (14+-1F »): peers... (60) 
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In (60) and (61) the suffixes indicate that the values of « and quantities 
depending upon it appropriate to the harmonic. If (58) and (59) are now 
inserted in the right-hand side of (56) we have, omitting the constant 
factor R?/p, 


(+2¢/R) 2 =A, b A, eittt A* ent A, ezintt 4 o* ea2int 


A E 
ae =a (G4 e? ntt Ci * e tnt, e22 nttO%* etn 


= 


1 : 
= (A, ent Ai e-? RCT ett O,* e-intt OC, e2int 
0 


+-C,* e7?nt) 
1 f Ff ; F : 
+ om (A, e2intt 4 * ee Cs e'tttO,* e-? ntt O, e2int 
0 
ees th Bane er ae aes te ee gt ee 
Now consider the term in 1+ €/2R on the left-hand side of (56). Assuming 
W=Up tw, e' "tt w,* ett we e2tnt_t w.* e-2int 


this becomes 


R DEN OU ee . : z 
ope 1+ R OF =734 (w, e'tt_w,* ent) + 2734?(we e2t ntt wo* e-2in ty 
1342 : ; 
a Pi (wy ef tt_y,* eC, ett, ett, e2int 
0 


+C,* e722 ft) 


203q? 


; (w,e?*"?'_w,* e-2t ay) CO e? ntt@* ef", e2int 
cn i oe 
+C 2 © ayy 
If now (60) and (62) are inserted in (56), and corresponding powers 
of erst are collected, a set of equations for wp, W,, Wy", Wy, W.* can be 
written down. The terms independent of e’"! give the equation 


d?w ldw A, hk? 1 R2 ‘ 
7 : wy oo we aa FACT 4 Car sts) 
0 


a 
ia {7 _* ; * 
y Y ayy . Y ; : « 
4 em (07, —Cyw,*+2(C,*w.—Cawe*)} .. . . . (63) 
and from the terms in e'”!, e2'”"¢ give 


dw, , ldw, R2 


2 
+1302, = za ni + 


(4(0,+4,0;*+4,*C,) 


dy | y dy Herat ies 
nae Rua. : 
—- rn (—C QW, +20 1°), 5 o z A . in (64) 
dw, . ldw, A,R? FR 


er ecrsmeer 2) | Ore 8h oe 
apy dy 1" 


; ; tan 
Tog 0 eater pela? (GOS 


An approximate solution, correct to order 1/c) can be obtained by 
inserting the known forms for w, and w,*, wy and w,* in the right-hand 
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sides of these equations. If this is done the integration can be carried 
through and the result expressed in terms of functions already known. 
The equation for w) becomes 


dwy  ldwy  AgR? 
dy?! y dy 


ih 
4p gee + 4G. Oat) 
Hea BS de { ae | 


Cy | para? 11J,(ai8?2) 
It will be seen that, with the exception of the main term, the constant 
terms cancel, so that this becomes 
dw  ldw, A R? R 
dy? © y dy Cope 


IAC eee es 
{0,*4 im PDP soya, ae To(ai-3P) f° 


(66) 


In this equation only the first pair of terms has been retained, the other 
pair being identical in form. 
Since for any Bessel Function of zero order 


d? iL a oe 5 io 
(i a -a) Tilly) = 1p) 


the solution of (66) which is zero at y=1 is 


Apher 1 Re OA a) 
DS tar rrr head ree gle a2 ; J o(ai3/) 
] R? OUAt. J (ae-8/2y) 
(oe ele ae > < |— —__-; >. 3 67 
Co fe m1" 5342 1 J (ai-3/2) (67) 


If in this we replace 


Jo(ai?!2y) 
he 1—F 
1 oy} (ous!) by 10 


7—3/2, 
and {1— cee by 1—Ff,,* 


) bod 7 2) 


we shall have the expression for the average velocity. 
Inserting the vaules of C, and C,*, and writing 


X,+1Y,=N, exp (204), —1 =k exp (71) 
the contribution to the average velocity from this pair of correction 
terms is a eR! 
4 j fee MeN ee (68 
My) - My - Copeat COS (€49’ —€19 +21 +71) (68) 


and there will be another contribution from the other pair of terms, 
identical in form, but with values appropriate to the second i 
This correction is in opposition to the main se When = uhis 
expression can be simplified, for then ky=3N 7M) and y,=204—€49 - 
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Substituting these values, and replacing R?/u0? by its equivalent 1/pn, 

(68) becomes, when o=3, 

—.— ——— COS (303—€49 Js eC CC (8D 

C4 pn? Ny COs (35;—€ 9 ) (69) 
In the equation for the first harmonic, the constant terms cancel each 

other in the same way, and the equation reduces to 


Pw, , 1d R 1 R? 1 R y Fol%st*7y) 
Pw, bec ity Bre wero A ean va re 2A y*ny * F053) 


dy? ' y dy 
1 R? *4 J (a gt*/*y) 


ae 70) 
Co a J o(%9t/?) ( 


The solution of this equation can be written down at once by repeated 
application of the well-known result : 


dw  ldw Vike Jolly) 
a aye ty dy TEMA) 


dy?” y dy 
then the solution of this which is zero at y=1 is 
— A fJolly) — Solky) 5 
SES? Saha Te) ae ae 


w (71) 

In imposing the condition that all the correction terms vanish at 
y=1, a further approximation is being made. Physically, it enforces 
the condition that the motion of the wall is due to the main terms only 
and the correction terms have no effect. Since the correction terms are 
small, this approximation may be adequate. For complete consistency, 
the arbitrary constants in the expressions for average velocity which 
are substituted in the equation should be left ‘ floating’ and the fully 
corrected solution substituted back in the equations of motion of the 
tube. The ‘ frequency-equation ’ would then be non-linear, and the 
pulse-velocity would depend on the particular form of the pressure 
function. The same situation will arise if a similar method is used to 
calculate the inertia-term correction. 

It is convenient to have a shorthand notation for 1—J(« ,yi®/?)/J9(«i3/2). 
This will be written Fy(«,y). Its conjugate Fy*(«,y) will also be needed. 
The full solution of (70) using this notation, is 


R?A, 1 RAC 
W4= (7 Co us Taig) Fol) 
Hy Gea ier ; 
dite (Poles) —Fo®(oy)} 
No ieee ; 
coaidge  Polery)—Fo(aay)}. - - . (72) 


The expression for the average velocity can be written down at once 
from this, by replacing Fo(x,y) by 1—F) wherever it appears, except 
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in the main term, where it is to be replaced by 1-++-7,F 49. If this is done. 
and the values of C, and C’,* inserted, we then form the expression | 
Wy eft w,* enint 
which will be the first harmonic of the average velocity corresponding to 
the pressure gradient in its real form (60). We have 
W eft w,* e-int 
M, 1 M,M 
ee ee Nee " 0 ib 
Mp” M9" sin (nt-+$1+€49")— Co np? 
1 MM, 
4c, np" 
Leen 
6, np? 


NM yo’ Myo’ cos (nt+4,) 
KN 1M oy’ M 9 {cos (nt-+,)—cos (nt+5)} 


KN M5 {Mi cos (ni+d4)—Mo) cos (nt-+y;5)} 
(73) 
where we have written 
$4=$1—91+€19 + €10 > 
te=$2—$1—92—Y1 +29 +10» 
Ps=po—2€ 49; 
ty=b2—$1 81+ Y2—€10 E10» 
$s=Py+€29 —€10 + 
The main term in this expression is shown negative. This is correct, 
since the pressure gradient has been taken positive. The correction for 
the second harmonic can be obtained in the same way. It is convenient, 
when computing the corrections for a practical example, to divide through- 
out by the modulus of the main term, and this has been done for the 


correction to the second harmonic, which is shown below, with the common 
factor 1/npcy omitted, and without the main term 


M,N 2M 29' cos (2nt+-#) 


1 Ma. 
+ gure, v7 {M y9' cos (2nt+-7)—M 9’ (cos 2nt+y,)} . (74) 
20 


where pe=b2—S2t+€ 20 +10 > 
bp=g= 261-811 F 10 +10 - 


When c=} these expressions can be simplified slightly by substituting 
for K, and K,, but this substitution is hardly worth making when 
r numerical calculation, each correction is expressed as a 
s of cos nt and sin nt, and these are summed and applied 
as proportional corrections to the components of the main term. 

As an example of the magnitude of a typical set of corrections, a 
complete calculation has been done for an experiment by McDonald 


computing. Fo 
set of coefficient 
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(Womersley 1954) on the femoral artery of the dog. Fourier Analysis 
of his pressure gradient record gave 


— a —0-159-++0-774 cos (nt-+0° 39’)+1-317 cos (2nt—82° 45’) 
—0-743 cos (3nt-+ 26° 30’)—0-414 cos (4nt—16° 39’). 


These coefficients are in mm of mercury per centimetre. The conversion 
constant to bring them to absolute units was included in the common 
factor 1/c np. 

It is not possible to make an accurate estimate of 1/cgnp until accurate 
measurements have been made of the pulse-velocity over short lengths of 
artery (see §6, below). The author was permitted to examine the film 
record of this experiment, and made a rough measurement of the maximum 
arterial expansion on the projected image. This gave a maximum 2€/R 
of about 6%. Since 
Sin |a—iY | 

0 


ae 


2€ 


R 


and the maximum average velocity was 88 cm/sec, this gives 


Cy=1310 cm/sec 


approximately. The pulse-velocity, estimated from records taken on 
other experiments, suggested a value of cy of about 850 cm/sec. 

In the calculation of the corrections, cg=1000 cm/sec was used, with 
o=}. Table 4 gives a list of the coefficients for each harmonic, corrected 
and uncorrected. The effect on the steady component, which is not 
shown in the table, was to reduce it by about 2 cm/sec, rather more than 
14%. This result is particularly important to the physiologist. It 


Table 4. Hxpansion-Correction 
Fourier Coefficients for Calculation of Average Velocity 


Uncorrected Corrected 


m | Coeff. of cos mnt | Coeff. of sin mnt | Coeff. of cos mnt | Coeff. of sin mnt 


1 +19-08 +33-14 +19-07 +28-74 
2 —31-78 +14-89 35:13 +16-05 
3 — 8-79 —10-58 —10-38 —10-99 
+ — 0-44 — 5:86 + 1-01 + 3-29 


means that in a system such as this, in which the oscillatory terms are 
large, and the steady component small, Poiseuille’s law cannot be used 
to calculate the steady flow from the pressure gradient—nor can the 
constant pressure gradient be inferred from the observed rate of flow 
or average velocity. Figure 5 shows the curves of average velocity 
over a complete pulse-cycle, plotted from the coefficients of table 4. 
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The observed average velocity is also plotted, for comparison. It will 
be seen that the general shape of the curve is brought nearer to that of 

_the observed curve by the correction. A rough calculation shows that 
to reduce the forward maximum on the calculated curve to the same as 
that for the observed curve it would be necessary to assume a wave- 
velocity of less than 500 cm/sec. This would not be consistent with the 
amount of arterial expansion observed. 


Fig. 5 
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§ 6. ConcLUSION—THE NATURE OF THE APPROXIMATIONS 


In the approximations so far made, the neglected terms are of order 
n*R?/c? or smaller, with the exception of the inertia terms in the equations 
of motion of the liquid. The largest of these might seem at first sight 
to be of order nR/c, but in fact they are larger than this, and their neglect 
is the major blemish on the theory outlined above. It is proposed Be 
remedy this in Part II of this communication, and to make the approxi- 
mation complete to order nR/c. Should it be necessary o ee ae 
approximation a stage further, and include terms of order Q ae u 
term in d?w/dz2 must be included, and the simple anions 
I,(nRy/c)=1, I,(nRy/c)=nRy/2c can no longer be used. The papers 
for the wave-velocity will then involve two non-dimensional pat amoters, 
x and nR/c. The larger values of nR/c which make this aee J a 
not likely to be important in the study of the circulation of the aes 
In the larger animals increased values of R and lower values otf ¢ are 
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offset by a lower pulse-rate. The range of values of « required to cover 
all conditions in the larger blood vessels in most animals is not very 
great, probably 1<«<12, though more accurate information is needed 
to check this general impression. 

This does not mean that the study of the motion at higher frequencies 
is without interest. On the contrary, the question of the limit of stability 
of the motion of some importance, for the light it may throw on the 
incidence of turbulence in a straight pipe. No pipe is completely rigid, 
and something may be learnt about the incidence of turbulence by 
considering small perturbations of a steady pressure gradient, allowing 
for the elasticity of the walls. 

To test the theory experimentally, the first requisite is an accurate 
measurement of pulse velocity in a short length of artery, say 3 or 4 cm, 
over which the diameter is reasonably constant. A method already 
exists with sufficient potential accuracy, but does not seem to have been 
used. The counter chronometer is now a standard piece of equipment, 
and such an instrument, working at one megacycle frequency, switched 
on and off by a pair of standard capacitance manometers, would be 
capable of measuring pulse-velocities to an accuracy better than one part 
in a thousand. Whether such accuracy could be achieved in practice 
is doubtful, owing to the rounded shape of the pulse-wave, making precise 
definition of its position difficult, but a manually adjustable d.c. gate, of 
the type used in some digital computers, could be used to define the crest 
of the pulse, and release a sharp-edged pulse from a trigger-circuit. Such 
measurements of pulse-velocity would be by far the most accurate ever 
made, and a detailed survey could be made of variations in pulse-velocity 
throughout the circulatory system. 

The degree of damping could be measured in two ways, either by 
accurate recording of the form of the pulse-wave at two points a known 
distance apart, or by simultaneous recording of the pulse-wave and its 
derivative. If A is a complex Fourier Coefficient of the pressure and A’ 
that of its space-derivative, then A’=—inA/c, and if c is complex, i.e. 
if co/c=X —7vY, then A’=A(X—iY) . in/cp, so that the phase-shift between 
the two will depend on the amount of damping. For measurements of 
this kind to be successful, the sensitivity of pressure-recording instruments 
must be raised. At present they are designed to measure pressures of the 
order of 100 mm of mercury, with a stable zero. The limits of the space- 
derivative measurement are of the order of +3 mm. If the amount of 
damping predicted by the theory proves to be greater than that observed, 
the simplest modification to take account of it would be to postulate 
different properties for the artery longitudinally and radially. To reduce 
damping, the artery would have to be freer longitudinally (so that its 
interaction with the moving mass of the fluid would be less) and stiffer 
laterally, a combination which is not compatible with homogeneous 
elastic properties. Such an elaboration is not justified until experimental 
proof is forthcoming that the degree of damping predicted is, in fact, too 
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large, and the viscosity of blood is known with greater accuracy than at 
present, for the accuracy of the estimate of « depends upon this. 

It will be shown in Part II that the corrections for the effect of the 
inertia terms, although appreciable, do not destroy the general picture 
given by the simple linear theory outlined above. 
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SUMMARY 


In measuring the optical constants of an absorbing medium two pieces 
of information are required and may be obtained by examining the re- 
flection of plane-polarized light at the principal angle of incidence. Since 
this angle is very close to 90° for metals in the infra-red, the value of 
measurements at smaller angles is examined. It is shown that there 
are serious disadvantages if reflecting coefficients only, such as R,?, 
R2 or the ratio p?, are measured. The ratio p of the reflecting coefficients 
of the p- and s-components and the relative phase difference 4 provide the 
most suitable parameters to measure at angles of incidence less than the 
principal angle. These are directly related to the modulus and the 
argument of the complex refractive index n—vk. 

Calculation of optical constants from the measurements is always 
tedious. A chart has been prepared whereby the real and imaginary parts 
of the dielectric constant may be obtained from p and 4 with an accuracy 
of 1%; the preparation and use of these graphs is discussed.  Illustra- 
tions are provided by quoting experimental results from 2 to 12 yu 
obtained with nickel, copper, silver, aluminium, mercury and antimony. 


§ 1. INTRODUCTION 


Wuen plane polarized radiation is reflected obliquely it becomes elliptic- 
ally polarized. The relation between ¢, the angle of incidence and y the 
angle of refraction is given by Snell’s law and the elementary conception of 
a real refractive index is extended for an absorbing medium so that an 
adequate description is given by a complex refractive index (Matossi 
1949). Thus we have 


sin d=n sin y=(n—7h) sin y=n(1—ix) sin x. 


The analytical justification of a complex refractive index is obvious yet 
its use requires care ; for example the surfaces of constant phase in the 
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absorbing medium do not coincide with the surfaces of constant amplitude 
(Schuster 1904, Ditchburn 1952). Two pieces of information are required 
to describe the behaviour of an absorbing medium and therefore two 
measurements are necessary. The real and imaginary parts of the 
dielectric constant, n2—k? and 2nk, have a more direct significance than 
n and k particularly in dispersion. 

In recent years some attention has been given to measurement of the 
optical constants of absorbing materials and various techniques have been 
developed (Bock 1945, Simon 1951, Avery 1951, 1952, Robinson 1952. 
Robinson and Price 1953, Bor 1952, Archard, Clegg and Taylor 1952, 
Conn and Eaton 1954 a, b). This revival of interest is in part a new 
application of improved experimental methods and in part stimulated by 
the successes of the band theory of solids which provides a well established 
framework for interpretation. No attention has been given to the 
peculiar difficulties presented by the study of metals in the infra-red since 
the work of Forsterling and Fréedericksz (1913). These difficulties have 
in the past proved such a stumbling block that, to the best of our know- 
ledge, no reliable measurements with any metals of high conductivity have 
been reported at wavelengths longer than 5. The purpose of the present 
paper is to indicate why experimental methods hitherto adopted fail when 
-applied to metals in the infra-red and how the difficulties may be overcome. 
It emerges that the best parameters to measure are pand 4. pis the ratio 
of the reflecting coefficients, k,, and &,, of the p- and s-components of the 
reflected radiation and J is the relative phase difference so that 


p exp (14)=R,/R,=cos 6+ x/cosd—y. . . . . (I) 


The conventions adopted in defining R,, R, and 4 are discussed briefly 
in an appendix since the accepted usage in optics is not self-consistent. 

A laborious calculation must be made before n?—k? and 2nk are obtained 
from measurements of p and 4 at each wavelength. A great deal of time 
is saved if these operations are carried out graphically. This also ensures 
that the significance of changes in the experimental parameters is immedi- 
ately grasped and the experiments can be pursued more intelligently. The 
most suitable graphical display is considered and extensively illustrated by 
practical examples. 


§ 2. MEASUREMENT OF THE OPTICAL CONSTANTS OF METALS IN THE 
INFRA-RED 


The optical constants of metals are best determined by studying 
reflected radiation. Examination of transmitted light is restricted by 
strong absorption of the conduction electrons to very thin specimens and 
the question arises to what extent the properties of such specimens accord 
with those of bulk material. In studies of reflection the relevant thickness 
is directly related to the ‘ penetration depth ’ and care must be given to 
the preparation of representative surfaces. 


Q2 
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The two pieces of information required must be obtained from measure- 
ments which depend on different functions of n and k. If both measure- 
ments are associated with essentially the same function of n and k the 
experiments cannot be successful. #,?, or R,R,*, and Rk, or RR,*, are 
plotted against ¢ in fig. 1, from which it appears that measurement in the 
neighbourhood of the principal angle ¢ offers great promise and much 
attention has been given to such measurements. If metals are to be 
examined in the infra-red however, ¢ usually exceeds 85° (Foérsterling and 
Fréedericksz 1913). For this spectral region the methods reported by 
Simon (1951), Avery (1951, 1952) and Conn and Eaton (1954 a) together 
with the method proposed by Abelés (1952) are of value in investigating 


Fig. 1 
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Intensities reflected by various media as a function of the angle of incidence. 


materials of lower absorption when ¢ is smaller but become very difficult 
when dealing with typical metals. The second method of Conn and 
Katon (1954 b) is relevant but laborious for routine measurements. A new 
technique must be sought and the first step is to assess the benefit 
of working near ¢ and the loss if ¢ is less than 4. 

Ratios of intensities such as R,?, R2 and p? are readily measured in the 
infra-red ; we consider R,? and p? since the former obviously yields 
essentially the data obtained at normal incidence and the latter provides — 
the additional information in a simple experimental form. Ry, the 
reflecting coefficient at normal incidence, is given by 
mo 12+ ke 1+R2 n?+k2+1 2 
TET EB a coho eae ° Kei ee aa 
Since both n? and k? are very much larger than unity it follows that 


R2= 
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Measurement of Ry yields (n?--k?/2n). Relevant values of the optical 
constants are listed in tables 1 and 2. 


___ Using the Fresnel expression discussed in the appendix it is easy to show 
that : 


— ’ —tan dsind 
Recap (3,\— oo" and Patna danas ; 
where s=/(n?—sin? d) or s?=n?—k?—sin? 6—2ink. 


If R,*, p* and s* are the complex conjugates of R,, p and s respectively, so 
that A,A,* or RZ is the coefficient of reflected intensity, we then have 
1+K,R,*  ss*-+cos*¢ 
1—R,R,* — (s+s*) cosh’ Cale oe De He 

1 * * 2 4 sin2 

ee + pp _ 88 -+tan me g (5) 

l—pp* (s+s*)tan¢dsing 
Equations (4) and (5) are exact and express measured parameters A and B 
in terms of the optical constants. If sin? ¢ is negligible compared with n?, 
we may write ss* as n?+-k? and s+-s* as 2n. Under these circumstances 
cos? ¢ is negligible compared with n?+-k2 and the above relations are 
equivalent to 


A= 


2A cos d==(n?-+k?)/n Ce te heer, Merten LO) 
and 
2B tan ¢ sin d= (n?+k?+tan? ¢ sin? ¢)/n 
=2A cos d+tan? ¢ sin? d/n. St See as) 


Equation (6), when compared with (2) shows that measurement of the 
reflected intensity of the s-component at any angle of incidence provides 
the same information as measurement at normal incidence, a conclusion 
which indeed is obvious from fig. 1. The second form of eqn. (7) reveals 
that unless the term tan? ¢ sin? ¢/n is comparable with 2A cos ¢ measure- 
‘ment of B gives no additional information which is significant, 
Conversely, provided ¢ is large enough to ensure that these terms are 
comparable, measurement may be made with profit. At the principal 
angle these two terms are equal since, with the approximations made 
above, 

tan? J sin? d6=n?-+-k? (Drude 1887). Plane een) 


These conclusions, based on eqn. (7), are illustrated in fig. 2. pp* or p? 


is plotted against ¢ for the various cases listed in tabie 1. The continuous 
curves are obtained by calculating p® using both terms on the right hand 
side of eqn. (7). If the second term is ignored so that the right hand side is 
2A cos ¢, the dotted curves are obtained ; these are therefore calculated 
from a knowledge of (n?--h?)/n only and it is necessary we work at such an 
angle of incidence that these curves differ significantly. If, as in the 
method of Avery (1951, 1952), the required information is to be obtained 
by measurement of p? at two angles of incidence, one of these must be near 


the principal angle. 
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Table 1 
Wavelength Material n k Reference 
‘03 u* — Cu (sputtered 0-43 5:6 Forsterling and 
: ipa a Fréedericksz (1913) 
2 5:10 pe Cu (sputtered) 2:98 28-4 3 
3 1-00 Ni (bulk) 2-62 5:26 — Ingersoll (1910) 


* These values of n and k are conservative since values much larger than 
j j 1 “ @ » ‘ 
these are commonly obtained in the infra-red ; see table 2. 
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p” or P,,?/R.? plotted against ¢ for the various cases listed in table 1. 


The 
numbers on the graphs refer to the order in table 1. 


The measurements which lead to these conclusions are of RZ and p? ; 
the difficulty arises because the moduli of both depend on the function 
(n?--k?)/n. To measure one of these is sufficient and p” is to be preferred 
for experimental reasons. The obvious associated parameter is the argu- 


ment of p, that is 4. It follows from eqn. (3), making the approximation 
that sin? ¢ is negligible compared with n2, that 


tan 4=— 2k tan ¢ sin d 


n*+-k?—tan? d sin? ¢ 
hence 


2 cot 4 tan ¢ sin 6=—(n2+-k? 
Equation (9) should be compared with eqn. (7). 


tan’ d sin? d)/k. . . (9) 

Since we are compelled 
to work at angles less than ¢ the essential quality of the new equation is 
that cot 4 is controlled by a new function (n?+-k?)/k. If we ignore the 
second terms in (7) and (9), that is assume that tan2 ¢ sin? d is much less 
than (n?-+-k*), an approximation which demonstrates the essential quality 
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of these relations but which is too inaccurate to use in computing, we have 
that 
n~(n?+-k?)/(2B tan ¢ sin ¢), Seep are we lL) 


aon pe) 2corAtandsind). . . . .. . (9a 
Hence 


: 1 
|x |~2 tan ¢ sin 6 / ( +tan? 4) and arg n~JB tan J. 


Measurement of p and 4 offers a means of measuring the real and 
imaginary parts of the refractive index even if ¢ is significantly less than 
o. Any two experimental parameters which give the required information 
must necessarily be related to the modulus and argument of the complex 
refractive index and therefore to p and 4. The conclusion is therefore 
that p and J, or their equivalents, must be determined if measurements in 
the infra-red are to be successful. The experimental methods developed 
as a consequence of this are discussed in the following paper. 


§ 3. COMPUTING THE OPTICAL CONSTANTS 


To deduce numerical values of n and k or of n?—k*® and 2nk from 
measurements of p and A requires laborious calculation. Two families of 
graphs have therefore been prepared whereby, knowing p and 4Z, values of 
n?—k? and 2nk may be obtained by inspection with an accuracy of the 
order of 1%. The use of these charts* is not of course limited to work in 
the infra-red. Moreover their value is not restricted to experiments in 
which p and 4 are measured: in many cases transformation of the 
experimental variables to p and A is simple (Conn and Eaton 1954 b). 

From eqn. (3) 


stan ¢ sin d6=(1+p exp ((4))/(lI—pexp(24)) . . (10) 
Hence ; 
2 ee cane! = Dp nth 
s _” k sin’* ¢ 2ink = € +p exp (0 7 SAWS OLON 
tan? ¢ sin? ¢ tan? ¢ sin? d l—p exp (t4) 
(11) 
By making use of complex conjugates it is readily shown that 
1+p?+2pcos4  1+y _ 2p cos 4 Pe 
ee ete AE SS herve y= : 12) 
eee ai 2, 008d I=y ee oi ( 
The phase angle is given by 
2p sin 4 ae 2x ee or, 
tan Q= es =z so that sin 2Q= Te and cos 2Q= eee (13) 
From eqn. (11) we then have 
n2—k?— sin? . : l+y 1—# 
 (n2— k?) = —tan? P . cos 2Q= ——.—.__ (14) 
X (n b)= tan? sin? d d l—y 1+ 
—2nk we I+y 2% (15 
and = Y (2nk)= eer =tan? P . sin 20= are D) 


*This is reminiscent of the well-known Smith charts ha in studies of 


electrical impedance. 
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where y and x are the simple functions of p and 4 given respectively in 
(12) and (13). These equations, which are exact and which can be derived 
_ from the relations of Drude (1887), have been used to prepare families of 
contours or graphs. 

The charts are shown on a greatly reduced scale in fig. 3. The abscissae 
are 4, the relative phase, over the range 0° to 120°, which may be read on 
the scale used in practice with an accuracy of 0:2°. With the convention 
adopted (Appendix A), 4 is negative. The ordinates are the amplitude 
ratio, p, extending from 0-100 to 1-000; these may be read to 0-002. 
Since the practical values of n?—k? and 2nk vary over a wide range, it is 
convenient to retain tan? ¢ sin? ¢ in eqns. (14) and (15) as a scaling factor. 
When ¢ is 51° 48’, 73° 10’ or 84° 19’, tan? ¢ sin? ¢ has the values 1, 10 or 100 
respectively ; it is therefore convenient to use one of these angles of 
incidence. If sin? ¢ is negligible, the presence of this scaling factor is the 
only influence of fd. One set of contours are the loci of constant values of 

Y=—2nk/tan? ¢ . sin? ¢ ; 
these are the closed loops which terminate at the top left hand corner of 
the diagram with the particular value of Y indicated on each. The second 
set of contours are very roughly orthogonal to these and are loci of constant 
values of 
X =(n?—k?—sin? ¢)/tan? ¢ . sin? ¢. 

Below the heavy zero line X is positive; above, it is negative. The 
measured values of p and 4 define a point; the values of X and Y are 
readily obtained by interpolation from the lattice of these curved 
contours. Over most of the chart, which measures 120 cm by 90 cm, X 
and Y may be obtained with an accuracy which varies between 1% and 2%. 

In preparing these curves subsidiary graphs were first drawn on a large 
scale ; these were of log,, X and log,) Y against p for various values of 4 
and against 4 for various values of p. In those regions where log) X and 
log,) Y become infinite, subsidiary graphs were drawn of X or Y against 
p and A. The accuracy of these must be somewhat greater than that 
available on the final chart and in all cases it was $%. In this way a 
large number of individual values of X and Y were plotted and the loci 
defined by particular values of X and Y were then drawn. / 

The practical accuracy for a range of materials and wavelengths is 
illustrated in table 2. The quoted values of the real and imaginary parts 
of the various dielectric constants have been measured in this laboratory 
and will receive fuller consideration in reports to be published shortly. 
In selecting samples for this table, an attempt has been made to illustrate 
the use of these charts over a wide range of p and 4. 


§ 4. GENERAL COMMENTS 


The merits of such a chart are not restricted to the speedy determmation 
of X(n?—k?) and Y(2nk) following measurement of p and 4. i he 
° ; ia ats 

accuracy of either n?—k? or 2nk may be assessed by inspection. ‘Thus 


230 ~=s J. R. Beattie and G. K. T. Conn on Optical Constants of 


when 4 is about 90° and p about 0-4, the precision with which p is measured 
fixes the accuracy of n2—k? while 4 controls the accuracy of 2nk. This 
may be contrasted with a case in which 4 is about 40° and p is about 
0-3 when n?— k? is primarily affected by 4 and p controls 2nk. If particular 
precision is sought in either n?—k? or 2nk, a rough estimate of X and Y 
ensures that attention is focused on the relevant measurement. From 
the wider spacing of the contours it is clear that the highest accuracy is to 
be expected if 4 lies between 60° and 120°, that is in the region of the 
principal angle. Given the experimental errors in p and 4, the errors in 
the values of n?—k? and 2nk can be readily estimated. 


Table 2 


yay 2nk n2— ke? b 
Specimen Ainw| ¢ i? a chart cale.| chart  cale. $ 


Ni (bulk) 3-0 _| 84°19’ | 0-699 | 80-2°| 89-5 | 88-5 |—103-5 |—101-8 85°6" 
9-5 | 84°19" | 0-870 | 36-8° | 380 | 384 |—780 |—775 | 88°3 
Cu (bulk) 7-0 | 77°O’ | 0-973 | 12-3°| 380 | 382 |—1460 |—1457 88°32 
Ag (bulk) 4-0 | 84°19’ | 0-948 | 49-5° | 66-0 | 65-3 |—463 |—463-5 87°21 
Al (evap.) 2-0 | 84°19’ | 0-911 | 69-3° | 41-5 | 41-5 |—201-5 | —203-6 86°3 
Hg (liquid) 6-0 | 83°30’ | 0-526 | 44-2°| 300 | 297 |—2 —3 86°41 
Sb (evap. 
metallic) 6-0 | 84°19’ | 0-219 |111-3° | 53-0 | 53:3 | 52-5 51-8 | 83°26’ 
73°10’ | 0-506 | 18-2° | 54:3 | 54-1 | 52-0 53-1 
51°48’ | 0-806 | 5-2° | 52-0 | 52-4 | 51-0 52:3 
12-0 | 84°19’ | 0-496 |101-2° | 71-0 | 70-9 |—17-6 |—17-3_ | 83°22’ 


The specimens referred to as ‘ evap.’ were prepared by evaporation iu vacuo; 
the figures for antimony refer to the metallic phase. The figures listed in the 
column labelled ‘ calc.’ were the result of exact calculation. 


For metals in the infra-red n?—k? is almost invariably negative and 2nk 
is large. It follows from the former conclusion that p is greater than 0-42, 
and from the latter that 4 is small unless 4 is large. Moreover if 4 is 
small, p is large ; the accuracy with which n2—k? is obtained rests largely 
on 4 and of 2nk on p. The measured value of J may be increased and p 
reduced by studying multiple reflections. 


§ 5. CONCLUSIONS 


The initial stimulus of this investigation was to seek a satisfactory basis 
for measurement of the optical constants of metals in the infra-red region. 
Since p and 4 are essentially the modulus and argument of the complex 
refractive index these provide two different and distinctive pieces of 
information. 'This conclusion is of importance in any case where measure- 
ment at the principal angle of incidence is difficult or impossible. 

The labour of preparing a computing chart such as that described is 
considerable ; a supply of working copies has been prepared and is avail- 
able. With such a chart it is not only possible to convert experimental 
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measurements rapidly to figures of the real and imaginary parts of the 
dielectric constant but appropriate steps can be taken to ensure that the 
measurements are made at suitable angles of incidence. 
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APPENDIX A 


The Phase Difference A and the Fresnel Coefficients 


The expressions of Fresnel for the reflected amplitudes in the p- and 
s-directions, parallel and perpendicular respectively to the plane of 
incidence are usually written in the form (Jenkins and White 1951) 

R,=tan ¢—x/tan d+y and R,=—sin d— x/sindty. . . (A.1) 
In this form there is a difference of sign interpreted as a change in phase 
of z for the s-component while there is no change in the phase of the 
p-component if + y is less than 90°. At normal incidence however, the 
components are indistinguishable. 

Derivations of R, and &,, based on the elastic solid theory (Preston 1912) 
or on the concept of electrical impedance (Schelkunoff 1943) yield the 
same sign, contrary to eqn. (A.1). The simple experiments of Drude and 
Nernst (1892) show that light incident normally on a silvered surface is 
reflected with a change in phase of approximately 7; the reflecting 
surface isa node. Both components must indeed undergo the same change 

in phase at normal incidence since Haidinger fringes are obtained as an 
unique set of bands. Moreover experiments at grazing incidence such as 
those of Wiener (1890) or those made with a Lloyd’s mirror (Titchmarsh 
1941) show that light of both components suffers a phase change of 7 
in this case also. 

The discrepancy between these conclusions and eqn. (A.1) arises because 
of an inconsistency in definition. The positive vectors E; and E, illustrated 
in fig. 4 lead to the expressions (A.1). For the s-component both vectors 
are perpendicular to the plane of incidence and drawn in the same sense, 
a convention which conflicts with that adopted for the p-component when 
the angle of incidence is zero. This inconsistency is removed by adopting 
E,. in fig. 4.as the positive reflected vector associated with the p-component. 
It follows that geviy: 

R,=—tan d—y/tan +x and k= —sin d— x/sin o+y . 
and that 


(A.2) 


p exp ((4)=R,/R,=cos ¢+x/cos —x . (A.3) 
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Fig. 4 


PLANE OF 
INCIDENCE 


[llustrating the conventions adopted in defining the positive electric vector 
of a reflected ray. 


Fig. 5 
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Phase changes at various angles of incidence, (a) for a dielectric, (6) for an 
absorbing material. 


Fig. 6 


@) Sot 


The disposition of the positive reference vector when the angle of incidence is 
(a) small, (6) near grazing incidence. 
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With this convention both the p- and s-components suffer a phase change 
of 7 on reflection by a dielectric at normal incidence : the phase difference 
Aiszero. For an absorbing material it is easy to show that the reflecting 
coefficient at normal incidence and the associated phase changes are 

Rw=— _—— a and tan Ops —2k/(n?+k?-+1). te A4) 
A is zero as before and in the infra-red where n2 and k? are both large 
compared with unity, 5, , is but slightly less than z. 

Figure 5 shows the phase changes at various angles of incidence for a 
dielectric and for absorbing materials. It appears that while the s- 
component suffers a phase change of 7 at grazing incidence, the p-compon- 
ent does not. It is obvious, however, from fig. 6 that the positive sense of 
the reflected p-vector leads to destructive interference at the surface of 
the mirror; figs. 6 (a) and 6 (b) show the positive reference vector at 
angles of incidence which are (a) small and (b) near 90°. The direction of 
reference moves through 90° as the angle of reflection changes from 0 to 
90° so that the positive sense of the reflected vector is inverted when 
compared with that of the incident vector. This inversion in sense of the 
reference vector must be borne in mind. 


ASP HeN. DEX as 
The Use of Data from Measurements of Reflectivity 
If measurements of reflected intensities are used to determine the real 


and imaginary parts of the dielectric constants, we have from the exact 
eqns. (4) and (5), 


| s?—2ACr(s)+C?=0 sh ye PRS es) 
and 
[s?}—2B7r(s)4+-7?=0. . 7. 2 3. 4 4 (B.2) 


In these equations r(s) is the real part of s, C is written for cos 4 and 7 for 
tan¢.sind. The solutions of these equations are 
Ink=2l (m—P?)1/? and n?—k?= 21m -+sin® ¢, 
where | 
m=|s |? =T0(AT—BC)|(BT—AC) and J=r(s)=4(T?—C?)/(BT —AC). 
These expressions although new in form are equivalent to those given by 
Abelés (1952). They imply measurement of R2 and p? at a particular 


value of ¢. f ee ae 

The convenient experimental method of Avery (1951, 1952) requires a 
solution of B2 for different values of f. This is possible only by graphical 
means or by assuming that sin? ¢ is the same for both values of 4 (Hodgson 


1954). 
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SUMMARY 


For metals the principal angle of incidence in the infra-red is near 90° 
and it has been shown that the real and imaginary parts of the dielectric 
constant are best obtained by measuring p and4; p is the ratio of the 
reflectivities of the p- and s-components and J the relative phase difference. 
pand 4 may be obtained without using birefringent materials by measur- 
ing two ratios of intensities, a technique which provides a cross-check. 
The experimental methods and arrangements developed to measure a 
number of metals in the spectral range 2 to 12, are described ; this 
range can be extended where necessary. Use of multiple reflections is 
discussed and errors are examined. Measurements are reported of 
evaporated aluminium mirrors since these are frequently used. 


§ 1. INTRODUCTION 


In the sixty years which have elapsed since the definitive work of Drude 
on the measurement of optical constants numerous authors have reported 
results of experiments on metals but very few of these refer to the infra- 
red region. The measurements of Forsterling and Fréedericksz (1913) 
extend to 5, those of Ingersoll (1910), Bock (1945) and Bueche (1948) 
to 2-5. Data on the reflectivity at longer wavelengths are avail- 
able, largely the work of Hagen and Rubens (1904), but such measure- 
ments of a single parameter are of limited interpretative value. The 
figures of Forsterling and Fréedericksz have necessarily received most 
general acceptance but as Wilson (1936 and 1953) has observed there 
are serious discrepancies with theory. 

The success of the modern theory of metals and developments in ex- 
perimental technique both emphasize the need for a study of the optical 
constants of metals in the infra-red. The problems peculiar to such 
measurements are examined in a paper (Beattie and Conn 1955) which 
shows why previous experiments have been confined to relatively short 
wavelengths and the manner in which these difficulties may be overcome. 
The measurements required of radiation reflected by a specimen are p, 
the ratio of the: eflectivities R, and FR, of the p- and s-components, and 
A the relative phase difference. The technique is thus an adaptation of 
one of Drude’s methods to suit present limitations m the infra-red. 
Since suitable birefringent materials are not available optical elements 
which give plane polarized light must be used. 


* Communicated by Professor W. Sucksmith, F.R.S, 
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In the present paper the experimental methods which have been 
developed are discussed and sources of error are examined. Practical 
experience has been obtained in studies of silver, copper, nickel, antimony 
and mercury over the range of 2 to 12 u and there is every prospect of 
extending such measurements, where necessary, to longer wavelengths. 
Detailed reports on these metals are now in preparation. Practical illus- 
trations are given in the light of this experience and by reporting measure- 
‘ments of the optical constants of aluminium mirrors prepared by 
evaporation in vacuo since these are widely used. 


§ 2. EXPERIMENTAL METHOD 


Plane polarized radiation is obliquely reflected by a specimen and 
passes through an analyser to the monochromator and detector. The 
intensity transmitted by the system ASP of fig. 1, [(p, 4,4) for azimuths 


wp» and #4 of the polarizer and analyser respectively is of the form (Conn 
and Katon 1954 a), 


I (bp, by) =1) (sin? Yp sin? ys, + p? cos? bp cos? J, +4p sin 2s» sin 2, cos J). 
Pp Ya 0 P ATP P Yat oP P ‘A 


(1) 
For metals of high conductivity R, differs little from unity and is con- 
veniently incorporated in J)=}H?R?,; EF is the amplitude of the 
incident light. 


Optical path from the source N to the entrace slit of the monochromator So. 


A, 8S and P represent respectively an analyser, a specimen and a polarizer 
(see text). 


_ Measurements of the transmitted signal are required at three different 
azimuths at least to obtain p and 4. At a fixed yp the intensities / , and 
/, are measured when the analyser is perpendicular and parallel to the 
plane of incidence so that 

1 =I (Pp, 7/2)=Iy sin? pp ; [,=I (hp, 0) =I yp? cos? isp. 
Hence p=tan pp 1/(I,/I,). P+. 64 Un 
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bp is chosen to ensure accurate measurement of p- To measure 4 the 
third term on the right hand side of (1) should clearly be as large as 
- possible ; the azimuths should therefore be +7/4 and thus 


1,=1(+7/4,+7/4)=1,(1+ p?-+ 2p cos AVA ep 52 (3) 
Then cos A= {2I,/T,—(1+p?)}/2p. 2. 2 2... (4) 
Also 1,21 (40/4, 7/4)=1,(1+ p?—2p cos 4)/4 (5) 
and cos A= {(1+p?)—21,/I,}/2p. (6) 


In principle either (4) or (6) may be used to measure J but neither is 
satisfactory. Since both p and J are required to provide two pieces of 
information such as the real and imaginary parts of the complex refractive 
index, it is desirable that the measurements are independent which (4) 
and (6) are not. If (4) and (6) are used together it follows that 


1 
00s A=3 (p+ =) {lp L)(lp+1a)} to b-Oae 2 shit 7) 


With the values of p usually encountered 4{p+(1/p)] does not depart 
seriously from unity; for example if p is 0-8 the value is 1-025. In 
general 


é il 
if 0-40<p<l1 then 1-45>4 (o++) elt 
p 


Equation (7) offers a means of measuring 4 largely independent of p, 
a feature which is particularly relevant to the small values of 4 found 
at long wavelengths when p tends to unity. 

Equation (1) is symmetrical in %, and %p so that the polarizer and 
analyser are interchangeable. This is a matter of some experimental 
importance since a prism monochromator shows such pronounced 
polarization that it has been used to estimate dichroism (Hyde 1948). 
With the simple prism instrument used in these experiments the ratio of 
the transmitted intensities at 2. for radiation polarized in the vertical 
and horizontal planes is about 0-70. The difficulty is eliminated if the 
azimuth of the polar between specimen and monochromator is fixed, 
that is, it fulfils the role of polarizer (%p). The first polar is then treated 


as the analyser (y,). 


§ 3. AZIMUTHS OF THE POLARIZER AND ANALYSER 


tp and ys, are measured with respect to the principal plane of the 
specimen and the origins of the two azimuth scales can be located 
geometrically with an accuracy of about 1” of arc. Greater precision is 
essential if systematic errors are to be avoided. The zero errors opp 
and 8, of these scales can be measured by examining any convenient 
specimen at such an angle of incidence that 4 is approximately 7/2 and 
at such a wavelength that the signal is large. Nickel has been used at 
an angle of incidence of about 84° and a wavelength of 2-23 these 
values are suitable for hand-polished, electro-polished or buffed nickel. 
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In the first experiment to measure 5), the polarizer is set at a nominal 
azimuth of 7/2 and intensities are compared for two settings +,’ of 
the analyser. Inequality of the two signals reveals a systematic error 
in the setting of the analyser. 

If pp=7/2—Spp and py=$,’— dp, 
eqn. (1) gives to a first order approximation, 

I(7/2—8 hp, fy' —8y4)=I, sin? by’ {1—2 cot fy’ (64,—p cos A . dyhp)} 
and — ) 


I (7/2—8p,—,’— 8b, )= Ty sin? 4’ {1+2 cot f,’ (844—p cos A . dpp)}. 

disp disappears if 4 is 7/2 so that 
dys, =(1—«)/4 cot by’ 

where « is the ratio of the measured intensities. If %,’ is between 30° 
and 60°, 84, may be measured to +0-1°. 

After applying this correction so that di, is zero, dp may be measured 
in a number of ways. The most convenient is probably to examine 
discrepancies in measurements of p?. Following eqn. (2) we define 


B=1(7/4—dphp, 0)/I(7/4— dybp, 7/2) 

y=1(—1/4—84fp, 0)/1(—n/4—8ypp, 7/2). 

Ignoring second and higher order infinitesimals it follows from eqn. (1) that 
B= p'(1+45yp) and y=p2(1—45pp) 

hence dsp=(B/y—1)/8. 


In practice 5, may also be measured to -+0-1°. It is a merit of this 
technique that errors in azimuth may be measured using specimens under 
examination. This is indicated above by quoting figures for nickel. 


and 


$4. Tok REFLECTANCE RATIO 
In examining the use of eqn. (2) to measure p, errors in azimuth will 
be neglected since the most important source of random error is fluctuation 
of the signal 6/. 
pe ahi | | / ees gae2 2 | 2) / 
dp/p= 301 (1/1, + 1/L.)= 361 (cosec? bp+-sec? yp/p*)/Ip. 
The smallest value of 5p/p occurs when \/p=tan yp; then 


. ive 
Sp/p=d1 ¢ ab -) [2-Te 


This optimum condition requires a change of yp with each specimen and 
it is experimentally more convenient to set %p=7/4 when 


Splp—51(1 + 5) /Ty WF SONI 


Taking as a practical figure 1/500 for d5//J), eqn. (8) indicates that 5p/p 
is about 0-4°%, when p is of the order of unity and 1-5°% when p is as low 
as 0-4; if the optimum condition is obeyed the corresponding figures 
are 0-4°%, and 1-2%. In the infra-red p is relatively large (Beattie and 
Conn 1955) and tends towards unity as the wavelength increases. 
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§ 5. Tue Revative Pass 


The influence of noise may b 38 iti Spee 
. y be assessed by writing unity in ec 7 
for 3|p+(1/p)]so that 2 ; eae) 


, cos A ==(1—J,/I5)/(14-1 4/13) and 84=8(L4/I5)+/(L4/1)/(L 414/15). 
Since errors in J, and J, must be added it follows that 


Saas) VG) iG) @ 


The second form of eqn. (9) is suitable for 4<7/2 when we write 

: d1/I,= 1/500 

as in §4; the third form is used if 4>7/2 when 8///, is set equal to 
1/500. In effect the assumption is made that either J, or ,, whichever 
is relevant, gives a full scale deflection. 


Fig. 2 
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Error in the relative phase 4 plotted against 4. 


The error in J is shown as a function of 4 in fig. 2 where 62 is given by 
§4=1/(5001/a) with a=(1—cos 4)/(1+-cos 4). 

The need to work at reasonably large angles of incidence is obvious since 
the accuracy deteriorates when J is less than about 15°. This is inherent 
since it is cos J which appears in eqn. (1), but multiple reflections may be 
used to overcome the limitation. This technique is discussed below. 
In practice 85° is the maximum angle of incidence which can be used 
and as the wavelength increases 4, of necessity therefore, tends to zero. 


§ 6. GENERAL COMMENTS 


Four measurements of the signal are thus made at each wavelength, 
yp being 7/4 in all cases. Two of these, 4,9 and 7/2, give p using 
eqn. (2); two, $,—-+7/4, give A using eqn. (7). pand 4 are parameters 


R.2 
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of an idealized ellipse, that is of the elliptically polarized radiation re- 
flected by the specimen. Measurement in practice gives values averaged 
over a range of wavelengths and over a range of angles of incidence. If 
the radiation is brought to a focus on the specimen the surface area may 
be reduced without loss of flux. Since only three of the four measured 
intensities are necessary it is possible to check that the range of wave- 
lengths and of incidence angles is not excessive. From eqns. (2), (3) and 
(5) it is easy to show that 
I,t1,=1,+1). 

It emerges from experience that this equation should be obeyed to 
within 1 or 2%. 

The symmetry of eqn. (1) is used from time to time to check the zeros 
of the azimuth scales. These may alter by accident or because of sagging 
or cracking of the selenium films which form the polars. The relation 


(bp, 4 )=L(—pp,—,) or I(+)=1(—) 
should be obeyed. If it is not obeyed, 
T(b’p, Bb )=L(bpt+ depp, by + dyq) 
=L (bp, fa)+(OL/dxp)dpp+ (SL/dy54 dy, (10) 


L(—P'p, — Wy )=L(— ppt dippp,—ha + dt) 
=L (bp, ba) —(OL/dybp)dpp— (51/dys4 )dy4. (11) 
By averaging the signals /(+) and J(—) the influence of these errors in 
azimuth is eliminated at the expense of doubling the number of readings. 
This means of dealing with such errors, not that described in § 3, was 
used in measurements of liquid mercury. 


and 


§ 7. EXPERIMENTAL DETAILS 


The source of radiation, a Nernst glower N in fig. 1, is inside a water 
jacket to reduce thermal fluctuations of intensity. Light emitted 
obliquely from the edges of a filament is partially polarized (Latie and 
Martens 1907, Czerny 1924); a measurement of the ratio of intensities 
of such light, polarized in vertical and horizontal planes, gave about 
0-80. The beam is restricted to radiation from the centre of the source 
by focusing a sharp image on the variable slit S, and to verify that the 
light is then substantially free from polarization, the width of 8, should 
not influence p?. Alternatively the specimen is removed and the apparent 
value of p for ‘ free space ’ is measured ; this should not differ significantly 
from unity for all wavelengths. The second alternative is more thorough 
but the polars and the mirrors M, and M, must be rearranged. 

The efficiency of the polarizing elements must be high. A percentage 
polarization of at least 99-5°%, is obtained with each polar using six films 
of selenium inclined to the beam at the Brewster angle (Elliot, Ambrose 
and Temple 1948). The behaviour of such polarizing devices has been 
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examined by Conn and Eaton (1954 b) who showed the importance of 
eliminating multiple reflections between films. In the present arrange- 
_ ment the effect of these reflections is removed by carefully centring the 
image on the entrance slit of the monochromator, S,. Since the azimuth 
of the analyser is changed from one reading to the next the beam must 
be clearly defined by a stop to prevent fouling by the metal frame of the 
polar. This is the purpose of a circular diaphragm on the mirror M,, 
imaged by M, on the aperture of the analyser. The convergence of 
light on the specimen, which is placed at the focus conjugate to S,, is 
controlled by a rectangular diaphragm at M,. 

A turntable, free to move about a vertical axis and graduated in 
degrees and minutes, supports a specimen ; the plane of the front surface 
is defined by three ball contacts. This simple arrangement is adequate 
to examine light reflected once. Without increasing the angle of incidence, 
the accuracy may be improved at longer wavelengths, when p tends to 
unity and 4 becomes small, by studying m reflections between a pair 
of specimens. p” and mJ are then measured. Figure 3 (a) illustrates an 


Fig. 3 


i (a (b) 


Arrangement for the study of multiple reflections. 


arrangement whereby the number of reflections may be changed without 
altering the optical path (see Jamin 1847), and since two SUS UIO Ne are 
excessive for accurate measurement of the relatively low values of p in 
the neighbourhood of 2 p, fig. 3 (4) shows the case which permits euudy ie 
single reflections. In both cases the mirrors are enor on ge e 
carriages both of which may move on a parallel slide ral to he 
turntable ; the mechanical quality of the slide and of these i 
clearly important. The radiation is brought So a focus poeicen the 
mirrors, the half-angle of the converging cone being, about 4° or ep 
: Serious errors arise if ¢ is greater than 85 
1 
since the angle of incidence affects the calculated meee of oe pa 
imaginary parts of the dielectric constant by Cage erica ie Reet 
tan? ¢ sin? d (Beattie and Conn 1955). At eroauel ppaier ae i ~ conte 
ation would be essential but the width of the source anc a a i ee a 
meter slit vitiate attempts to define ¢ more secutataly a 4 + oe 
angle of incidence is chosen and by altering the soba am : Pi ces 
the number of reflections is increased to compensate for the diminut 


aperture ratio of about 1/50. 
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of 4 and (1—p) at long wavelengths. Multiple reflections have been of 
value for measurements of silver, copper, nickel and aluminium beyond 
6; for wavelengths less than about 12, and all save metals of the 
highest conductivity single reflection and a single specimen one or two 
inches long are adequate. 

Since there is no cogent reason to anticipate that the optical constants 
of metals change rapidly with wavelength in the infra-red, a mono- 
chromator of modest resolution has been used hitherto. As in the 
arrangement of McAlister (1941) the radiation passes twice through a 
60° rock-salt prism of base 3 cm. Over the range 2 » to 14 1 the resolution 
is of the order of 0-1 1 but a slit width as large as 0-25 » or even 0-5 uw at 
long wavelengths has been used because of the intensity lost by defining 
apertures and polarizers before the light enters the monochromator. 
The instrument has been calibrated and checked from time to time by 
using the well-known atmospheric absorption bands and cells containing 
benzene and chloroform. The effects of radiation scattered within the 
monochromator can be serious particularly on measurements of J,. At 
long wavelengths J, is small and it follows that a large proportion of the 
signal at 10, say, may arise from scattered radiation of shorter wave- 
lengths. Although attention has been given to this in the design of the 
spectrometer (Williams 1948), precautions must be taken against the 
possibility that about 1°, of the signal at wavelengths greater than 
9 arises from radiation of short wavelength. For measurements at 
wavelengths longer than 6 1 a selenium powder filter is used (Pfund 1933) ; 
this filter, indicated by F in fig. 1, transmits 4°, at 2. and 90% at 6p. 
At wavelengths greater than 9 « a mica sheet is substituted for the opaque 
chopper blade commonly used. As the result of these precautions it is 
believed that the data obtained are substantially free from the effects 
of light scattered in the monochromator. 

The beam is interrupted at a frequency of 2-67 cycles per second by a 
chopper C in fig. 1. The e.m.f. of the evacuated Schwarz thermocouple 
is magnified by aresonant galvanometer amplifier tuned to this frequency 
by adjusting the feed-back ; this is followed by an electronic amplifier 
(see Beattie and Conn 1954). The output is presented, rectified and 
smoothed, on a linear scale of 500 divisions. 

The longest wavelength at which measurements have been made to 
date is 13. This limit is set by intensity of the source and diminishing 
values of 4. J, decreases and accuracy deteriorates beyond 10. At 
shorter wavelengths changes of 1°% in (n?—k®) and 2nk can usually be 
determined. Since ¢ is fixed, cross-checks ensure that only minor errors 
affect p and 4 but absolute accuracy is affected by measurement of the 
mean value of ¢ ; an error of 0-1° for instance at 85° causes a systematic 
error of 4%, in (m?—k?) and 2nk. Though there may in consequence be 
some uncertainty in the absolute values, the variation of the dielectric 
constant with wavelength is not affected. This source of error may be 
reduced by decreasing 4 which is possible using multiple reflections, 
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§ 8. OpricaL Constants oF ALUMINIUM MIRRORS 


Mirrors prepared by evaporation in vacuo of aluminium condensed on 
glass are very widely used, particularly in the infra-red. It is therefore 
important to record the optical constants of such mirrors which are 
frequently used as a reference in measurements of reflectivity. Moreover 
in experiments such as these under discussion appreciable polarization 
may be introduced by mirrors ; for instance M, and M, in fig. 1 are not 
used on their optic axes and may therefore sontmbutd significantly to 
apparent polarization from the source. Measurements of such specimens 
reveal the orders of magnitude encountered in examining metals of high 
conductivity and give practical point to many of the comments made 
above. 

Evaporated aluminium was condensed at 5 to 104.v. per second on 
two identical glass plates at a distance of 20 cm from the source. The 
residual gas pressure within the chamber during evaporation was of the 
order of 10°-°mm Hg; a liquid nitrogen trap was used between this 
chamber and the oil pumps. The thickness of the film, as measured 
subsequently by interferometry, was 1550-50 A.u. The optical constants 
of these surfaces were measured from 2, to 12 within twenty-four 
hours of their formation. The angle of incidence was 84-2° and up to 
five reflections were used in covering this range. The results are listed 
in the table. 


Optical Constants of Mirrors of Evaporated Aluminium 


Ain p p Ae aa 2nk n k lige 
2-0 0-885 60:8 —268 76 2°30 16-5 0-968 
2:5 0-885 50°7 —401 130 3:22 20:3 0-970 
3-0 0-883 43-0 —565 213 4-4] 24:2 0-971 
4-0 0-894 34-6 —881 362 5:97 30:3 0-975 
5:0 0-899 28:5 —1285 604 8-19 36:8 0-978 
6-0 0-898 24-6 —1670 930 11-00 42-4 On9 Ty 
7-0 0-900 20:8 —2185 | . 1440 14-63 49-0 0-978 
8-0 0-907 18-5 —2730 | 1870 17-00 : 55-0 | 0-980 

* 9-0 0-909 16-2 —3305 | 2590 21-1 61 3 0-980 
10-0 0-910 14:5 3890 | 3405 25-4 é 67:3 pes) 
11-0 0-910 13-6 —4270 | 4025 283 ; te ] td = Neate 
12-0 0-910 12-3 4700 f 5130 33-6 | 76-4 0-931 


244 J. R. Beattie on Optical Constants of 


The absolute accuracy to be attached to values of n2—k? and 2nk, of the 
order of 4°%, is determined by the error in ¢ and therefore in tan® ¢ sin? ¢. 
Since & is much larger than n it follows that the error in k& is of the 
order of 2°, and, from the figures for 2nk, the error in n is also of the 
order of 2°. This error increases to about 5% in the region of 11 ~ to 
12 since errors in measurements of intensity become significant. The 
associated error in R, is }°%, or less, depending on the values of n and k ; 
the figures for R,2 have therefore been given to three significant figures. 
The functions n2—k? and 2nk are listed in the fourth and fifth columns of 
the table for convenience of reference since these are of greater value in 
interpretation than n and k. It may be noted that n?—k? is large and 
negative ; this is contrary to the assumption made by Hagen and Rubens 
(1904) in deriving an explicit expression for Ry. 

Pressure, rate of evaporation and temperature at which the metal is 
condensed all influence the quality of mirrors of evaporated aluminium. 
On exposure to air an oxide film about 15 A.v. thick forms very rapidly ; 
subsequent growth of this layer lasts about a month when the oxide 
layer is about 404.vu. thick (Hass 1947, 1952). The optical constants 
therefore depend upon a number of factors; aluminium oxide (Al,0,) 
in the form of sapphire is not transparent at wavelengths greater than 
about 6 (Chasmer et al. 1951). For the present purposes of illustration 
a fuller examination of these effects is clearly out of place. 

The polarization which is introduced by using such mirrors may be 
calculated using eqns. (5), (7) and (9) of the paper quoted above (Beattie 
and Conn 1954 a). If ¢ is small it follows from these equations when 
terms of the order of 44 are neglected that 

l— p=2nd¢?/(n?-+k?) and 4 = 2kd?/(n?-+-k?). 
Hence when ¢ is 4°, it follows that 1—p is of the order of 10-4 at 2 
and 4 about 2 minutes of arc. These quantities decrease as the wave- 
length increases. It is clear therefore that aluminized mirrors may be 
employed with greater freedom than has perhaps been the custom when 
polarized infra-red radiation is employed. 
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XXVIII. CORRESPONDENCE 
The * Boundary Effect’? in the Thomas—Fermi Model for Atoms 


By R. A. BALLINGER and N. H. Marcu 
Department of Physics, The University, Sheffield 


[Received October, 1954] 


THE unmodified: Thomas—Fermi (TF) method leads to a formula for the 
total binding energy of an atom which is not in good agreement with 
experiment. However, in an important paper, Scott (1952) has recently 
examined the reasons for the discrepancy, and has concluded that 
corrections must be applied for the so-called ‘ boundary effect ’ and for 
exchange. One of us (March 1953) has examined the method employed 
by Scott for calculating the exchange energy, and has shown how a 
modified electron distribution, in which some account of exchange is 
introduced, can be set up. In this work, the ‘ boundary effect ’ was not 
considered, and no modification of Scott’s binding energy formula was 
called for, the only effect on the energy being the introduction of a term 
of order Z, whilst terms of order Z4/3 are neglected in Scott’s treatment. 
The purpose of the present note is to make some remarks concerning 
the ‘boundary effect’, and to suggest a modification. to be made in 
Scott’s formula. Scott has, effectively, estimated the correction to be 
applied to the TF formula by comparing the results of the TF approxi- 
mation with the exact wave mechanical results for the case of particles 
moving in a Coulomb field. In the following we use a somewhat different 
method from that employed by Scott. If one considers Z particles, 
moving in a Coulomb potential Z/r (we use atomic units throughout), 
then the total energy when the lowest }Z levels are doubly filled and we 
have n completed shells is 
Fo DP eng Fe Sh ge 
where v is related to Z by 
n(n-+-1)(2n-+1 Z 
oot rele es 
We are now interested in an asymptotic form for when Z is very large. 
It is easy to show from (2) that 


n= (3) 1/37; 1/3__4 as (3)? 7-13 ca O(Z-9/3) ; ee (3) 
and hence by substitution in (1) we have 
B= — ($) 8278 + 472 2 (3) 28758 4 OZ), (4) 


On the other hand, the TF approximation yields the result —($)1/8Z us 


for all Z, and is therefore strictly correct only in the limit Z—> oo as one 
would expect from the nature of the method. 
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The term 4Z? in (4) is just Scott’s boundary correction, but we see 
that there is also now a term of order Z*/3, and this is not negligible to 
the approximation considered by Scott. It seems clear to us from this 
derivation that to be consistent this term must be included in the binding 
energy formula. It is approximately one third of the exchange energy. 
In this way we obtain the modified formula 


E=—0-7687Z7/3 4 $2? 0-2936Z°4 O(Z4/8), . . | (5) 
We show in the table results for five atoms considered by Scott. For 


details of the way in which the ‘ observed ’ values were obtained we refer 
the reader to Scott’s discussion. 


Binding Energies (Non-relativistic) 


Z 6 ©) 13 18 23 
‘ Observed ’ 37:8 99-7 242-4 528 945 
Scott 36-7 97-6 236-9 518 933 
Eqn. (5) 38-1 100-4 242-1 527 947 


It will be seen that the results are in quite remarkable agreement with 
the observed values, a marked improvement being shown over the 
already satisfactory results obtained by Scott. In fact, the agreement 
is much better than one might expect from the nature of (5). However, 
the treatment we have given at least suggests a reason for this, for it 
will be seen from (4) that terms of order Z4/°, Z, and Z?/> are absent. 
Whilst this will certainly not be the case with the formula (5), the results 
for the case of a Coulomb field make it seem very reasonable to suppose 
that these terms will be small, and hence that binding energies calculated 
from (5) will be more accurate than would appear at first sight. 

Tn conclusion, we ought to point out that we have not of course thrown 
any light on the validity of the approximation implied in the use of a 
Coulomb field for the evaluation of the boundary correction. It does 
seem, however, that this is a less drastic approximation than one might 
imagine at first (for example, by noting that the first term of (4) is 
—1-14Z7/3, whereas that of (5) is —0-77Z7/), for the bulk of the correction 
arises from the K-shell electrons, as was demonstrated by Scott. Concern- 
ing his derivation, the omission of the term in 2°? which we have displayed 
here arises from the limit » — co which he takes, this procedure giving 


only the term of highest order in 21’. 
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XXVIII. REVIEWS OF BOOKS 


The Physics of Particle Size Analysis. British Journal of Applied Physics 
Supplement No.3. (London: Institute of Physics.) [Pp.218.] Price, bound 
in cloth, 35s., post free. 


In this volume there are three survey papers and thirty-three original contri- 
butions with the subsequent discussions which were presented at a conference on 
‘The Physics of Particle Size Analysis ’ held at the University of Nottingham 
in April 1954. 

In recent years two new techniques have become widely used in particle size 
analysis : electron microscopy for the direct examination of fine dusts, and 
automatized counting and sizing with photoelectronic machines. Two sessions 
of the conference were devoted to the latter technique dealing firstly with the 
theory of automatic counting and sizing and secondly with the various photo- 
electronic instruments. Only one paper deals with the application of the 
electron microscope to particle size analysis. 

At other sessions of the conference particle shape factors and visual counting 
and sizing of microscopic particles were discussed, whilst the physical properties 
of solid—liquid and solid—gas dispersions were fully treated in the first three ses- 
sions. It is rather surprising that not one paper deals with the properties and 
particle size analysis of solid—solid dispersions—a field of great interest to the 
petrologist and metallurgist. 

The editors and organizing committee are to be congratulated on the pub- 
lication of this volume so quickly after the conference. Ji 


Journal of Nuclear Energy. (London: Pergamon Press Ltd.) [Pp. 300 
approx. per volume. Quarterly at first.] Price £4 10s. Od. per vol. for 
Libraries, £3 10s. Od. for individuals. 


THE appearance of a new scientific journal is not always welcome for there 
are already so many journals to read. The Journal of Nuclear Energy, however, 
is an exception: the new and important branch of science and technology 
which it sets out to serve had hitherto no suitable outlet outside the United 
States. Nuclear energy projects are rapidly developing in many European 
and Commonwealth countries, and it is entirely fitting that there should be a 
European journal embracing all aspects of the work. : 

As the “ Atomic Age’ gradually approaches the technologist begins to take 
precedence over the scientist. It is therefore to be expected that the first 
number of this Journal places its main emphasis upon the more technical 
aspects of nuclear engineering. Thus there are papers on the design of mechanical 
pumps for use with liquid metal coolants, on the control and operation of 
reactors, and on the economics of a large-scale power reactor programme. 
Nevertheless the more scientific aspects are by no means neglected. There are 
papers dealing with neutron cross section measurements and with the physical 
and physiological effects of radiation. 

The international editorial board are to be congratulated upon the high 
standard which is set by their first number, the publishers upon the excellent 
quality of their print and plate. Engineers and scientists throughout the world 
will surely welcome the opportunity to collaborate within the pages of this 
Journal. Bsa 


[Lhe Editors do not hold themselves responsible for the views 
expressed by their correspondents. ] 


